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follows : 

By 

To  Obtain 
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metres 
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pounds 
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0.04788026 
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(mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  metre 

square 

feet 

0.09290304 

square  metres 

tons  (force)  per  square  foot 

95.76052 

kilopascals 
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CONSOLIDATION  OF  SOFT  LAYERS  BY  FINITE  STRAIN  ANALYSIS 


PART  I :  INTRODUCTION 

1.  The  importance  of  the  ability  to  accurately  predict  the  con¬ 
solidation  behavior  of  soft  clay  deposits  is  manifest  in  the  millions  of 
dollars  spent  annually  in  the  disposal  of  materials  dredged  from  the 
nation's  waterways  and  wastes  of  the  mining  industries  involved  in  phos¬ 
phate  and  other  mineral  ore  production.  To  adequately  design  the  catch¬ 
ments  necessary  to  hold  these  vast  quantities,  knowledge  of  the  rate  of 
settlement  of  the  clayey  material  is  required.  The  economics  of  the 
disposal  operation  dictates  that  each  specially  constructed  area  be  used 
to  its  fullest  potential.  Therefore,  estimating  the  consolidation  in 
each  area  is  a  prerequisite  to  determining  the  overall  area  needed  to 
support  a  specified  application  rate. 

2.  Methods  currently  available  for  computing  the  potential  set¬ 
tlements  of  soft  clay  deposits  as  a  function  of  time  are  based  on  both 
empirical  and  theoretical  relationships.  This  report  will  deal  prin¬ 
cipally  with  the  theoretical  aspects  of  consolidation  and  their  applica¬ 
tion  to  the  settlement  of  soft  clay  deposits  under  self-weight  loading 
(although  the  theory  and  techniques  employed  are  equally  applicable  to 
other  types  of  loading  as  will  be  shown  in  a  practical  example) .  It 
should  be  noted  here  that  the  method  to  be  presented  is  limited  to  one¬ 
dimensional  consolidation  of  saturated  clay  deposits  which  in  actuality 
is  no  limitation  when  applied  to  the  large  wet  disposal  sites  in  current 
use.  Other  limitations  will  be  discussed  as  they  apply  to  particular 
solution  techniques,  but  in  general  the  theory  will  require  only  that 
the  clay  deposits  be  homogeneous  in  material  type. 

3.  The  first  theory  enabling  the  prediction  of  one-dimensional 
consolidation  in  soils  was  published  by  Karl  Terzaghi  in  1924.  The 
simplifying  assumptions  adopted  for  this  original  theory  were  such  that 
its  applicability  was  effectively  limited  to  the  consideration  of  re¬ 
latively  stiff  thin  layers  at  large  depths.  For  example,  the  assumption 
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chat  there  is  a  constant  relationship  between  void  ratio  and  effective 
stress  and  that  permeability  does  not  change  within  the  consolidating 
material  is  valid  only  when  the  ultimate  change  in  effective  stress  is 
small  in  comparison  to  the  preconsolidation  effective  stress.  Because 
settlements  in. soft  clay  deposits  such  as  dredged  fill  where  strains 
greater  than  50  percent  are  not  uncommon,  the  assumption  of  small 
strains  negates  the  usefulness  of  Terzaghi's  theory  unless  soil  param¬ 
eters  and  layer  thickness  are  continuously  updated. 

4.  The  usual  form  of  Terzaghi's  governing  equation  (Terzaghi  and 
Peck  1967)  is 


3u  3^u 

at  '  cv  2 

3x 


(1) 


where  u  is  the  excess  pore  water  pressure  and  c^  is  the  coefficient 
of  consolidation.  The  independent  variables  are  time,  t  ,  and  the  ver¬ 
tical  space  coordinate,  x  .  Even  though  this  differential  equation  has 
limited  applicability  to  the  general  problem  of  soil  consolidation,  it 
has  remained  the  popular  choice  among  geotechnical  engineers  because  it 
is  the  simplest  equation  and  is  taught  in  all  basic  soil  mechanics 
courses.  Solution  of  the  Terzaghi  equation  is  simplified  because  it  is 
linear  and  the  same  as  the  heat  conduction  equation  for  which  analytical 
solutions  for  a  multitude  of  boundary  conditions  are  available  (Carslaw 
and  Jaeger  1959). 

5.  Many  authors  have  offered  alternatives  to  Equation  1  to  better 
simulate  the  actual  behavior  of  soils.  Schiffman  and  Gibson  (1964)  as¬ 
sumed  that  permeability  and  the  coefficient  of  volume  change  were  known 
functions  of  depth  and  derived  the  governing  equation  as 

3^u  _1  dk  clu  ^wmy(X)  3u 

2  k  dx  3x  k(x)  3t  '  ' 

dX 


where  k  is  permeability,  y  is  unit  weight  of  water,  m  is  coeffi- 

W  V 

cient  of  volume  change,  and  other  terms  are  as  defined  previously. 

Davis  and  Raymond  (1965)  produced  a  nonlinear  theory  of  consolidation  by 
assuming  a  constant  logarithmic  relationship  between  void  ratio  and 
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effective  stress.  Their  governing  equation  is 


32u  _  /1\2  in  3c^ 

.2  V o ' /  3x  3x 
3x 


1_  3o_^ 
o'  3t 


(3) 


where  o'  is  vertical  effective  stress  and  other  terms  are  as  previously 
defined.  Other  theories  or  variations  include  the  works  of  McNabb 
(1960)  and  Mikasa  (1965).  However,  all  of  these  variations  to  the 
original  Terzaghi  equation  have  their  own  unique  limitations  and  are  not 
suited  for  application  to  large  deposits  of  soft  dredged  fill  or  mine 
tailings. 

6.  While  the  equations  of  McNabb  and  Mikasa  did  allow  for  large 
strains,  the  first  completely  general  theory  of  one-dimensional  consoli¬ 
dation  in  soils  was  published  by  Gibson,  England,  and  Hussey  in  1967. 
Their  governing  equation,  which  will  be  fully  developed  in  the  next 
section,  is 


hs  _  \  d_  fk(e)  *1  3e  3_  k(e)  do '  3e  3e 

l  /  de  |_1  +  ej  3z  3z  >^(1  +  e)  de  3z  3t 


(A) 


where  yg  is  the  unit  weight  of  solids,  e  is  void  ratio,  z  is  a 
material  coordinate  to  be  explained  later,  and  other  terms  are  as  de¬ 
fined  previously.  The  consolidation  equation  in  this  form  is  particu¬ 
larly  suited  for  application  to  thick  soft  clay  deposits  because  it 
intrinsically  includes  the  effects  of  self  weight,  permeability  varying 
with  void  ratio,  and  a  nonlinear  void  ratio-effective  stress  relation¬ 
ship.  It  also  is  independent  of  the  degree  of  strain  which  is  the  key 
reason  it  is  suitable  for  thick  soft  clay  deposits  susceptible  to  large 
settlements.  Hereinafter,  Equation  A  will  be  referred  to  as  the  finite 
strain  theory  while  Equation  1  and  its  variations  will  be  referred  to  as 
the  small  strain  theory. 

7.  The  fact  that  Equation  A  is  a  completely  general  theory  of  one- 
dirnensional  consolidation  was  demonstrated  by  Schiffman  (1980)  when  he 
showed  that  the  small  strain  theory  and  its  principal  linear  and  nonlin¬ 
ear  variations  are  all  special  cases  of  the  finite  strain  theory.  Prac¬ 
tical  application  of  the  theory  and  a  comparison  of  results  with  those 


of  the  small  strain  theory  were  presented  by  Gibson,  Schiffman,  and 
Cargill  (1981).  Using  conventional  laboratory  data  for  a  soft  marine 
deposit,  they  demonstrated  that  faster  and  larger  settlements  are  pre¬ 
dicted  by  finite  strain  theory  although  predicted  dissipation  of  excess 
pore  water  pressure  may  be  slower  than  that  predicted  by  the  small 
strain  theory. 

8.  The  next  part  of  this  report  will  document  the  development  of 
the  finite  strain  theory  governing  equation  along  with  the  initial  and 
boundary  conditions  necessary  for  its  solution.  The  solution  technique 
to  be  employed  is  an  explicit  finite  difference  scheme  which  will  then 
be  illustrated  in  a  manner  suitable  for  computer  programming.  The  com¬ 
puter  program  CSLFS  (Consolidation  of  Soft  Layers,  Finite  Strain)  will 
be  used  to  solve  a  practical  dredge  fill  consolidation  problem  and  a 
soft  foundation  consolidation  problem  to  illustrate  the  capabilities  of 
the  program.  A  user’s  manual  for  CSLFS  is  included  in  Appendix  A. 
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PART  II:  FINITE  STRAIN  FORMULATION 
OF  ONE-DIMENSIONAL  CONSOLIDATION 


9.  The  basic  assumptions  necessary  for  the  development  of  the 

theory  of  one-dimensional  finite  strain  consolidation  are: 

The  soil  system  is  saturated  and  consists  of  a  compres¬ 
sible  soil  matrix  and  incompressible  pore  fluid.  While 
the  soil  matrix  is  considered  compressible,  individual 
soil  particles  are  incompressible. 

b^.  Pore  fluid  flow  velocities  are  small  and  governed  by 
Darcy's  law. 

c.  There  is  a  unique  relationship  between  soil  permeability 
and  void  ratio  such  that 


k  =  k(e) 


<i.  There  is  a  unique  relationship  between  vertical  effective 
stress  and  void  ratio  such  that 


o’(e) 


e.  The  material  is  homogeneous  as  to  type. 

These  conditions  are  only  slightly  restrictive  and  imply  monotonic  load¬ 
ing.  The  usual  assumption  made  in  the  small  strain  theory  restricting 
the  magnitude  of  strain  is  not  made  here. 

Coordinate  System 

10.  The  election  to  allow  unlimited  strain  makes  the  use  of  a 
fixed  coordinate  system  impractical  due  to  the  relatively  large  movement 
of  the  top  boundary  of  the  consolidating  layer.  To  simplify  the  re¬ 
quired  mathematics,  a  coordinate  system  which  moves  with  the  layer  is 
needed.  This  condition  is  satisfied  when  the  coordinates  are  defined  in 
terms  of  the  volume  of  solid  particles  in  the  layer,  which  happens  to  be 
a  constant  quantity.  These  material  or  reduced  coordinates  (Ortenblad 
1930)  are  uniquely  suited  for  use  in  the  time-dependent  consolidation 


problem  because  they  are  time  independent  and  independent  of  the  amount 
of  strain. 

11.  Before  material  coordinates  can  be  employed,  however,  a  rela¬ 
tionship  must  be  established  between  these  coordinates  and  the  more 
conventional  methods  of  thickness  measurement.  Consider  the  soil  element 
shown  within  the  consolidating  layer  in  Figure  1.  At  time  t  =  0  the 
initial  configuration  is  given  in  what  will  be  called  Lagrangian  coordi¬ 
nates.  This  system  is  related  to  "real"  measurements  at  t  =  0  .  For 
time,  t  ,  during  the  consolidation  process,  "real"  measurements  are 
made  in  terms  of  a  convective  coordinate  system  which  is  a  function  of 
the  Lagrangian  coordinate  and  time. 


LAGRANGIAN  convective 

COORDINATES  (a)  COORDINATES  (?) 

t  =0  t>0 

Figure  1.  Coordinate  systems 

12.  Both  Lagrangian  and  convective  coordinates  are  a  measurement 
of  the  soil  system,  which  includes  both  solid  soil  particles  and  the 
pore  fluid.  As  previously  stated,  the  material  coordinate  is  a  measure 
of  the  volume  of  solid  particles  only.  A  comparison  of  these  three 
systems  is  illustrated  in  Figure  2.  As  shown  in  the  illustration,  only 
the  Lagrangian  and  material  coordinates  are  constant  for  all  time  for 
particular  points  in  the  soil  layer.  It  is,  therefore,  convenient  to 
develop  the  governing  equation  in  terms  of  either  of  these  systems.  Th 
material  coordinates  will  be  used  here. 

13.  Since  material  coordinates  are  not  measurable  in  the  usual 
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Figure  2.  Comparison  of  coordinate  systems 


sense,  it  is  necessary  to  develop  a  method  of  conversion  from  one  coordi¬ 
nate  system  to  another  so  that  the  layer  thickness  may  be  expressed  in 
easily  understood  conventional  units  at  any  time.  Consider  the  differen¬ 
tial  elements  of  soil  shown  in  Figure  3.  If  these  elements  are  chosen 


LAGRANGI  AN  CONVECTIVE  MATERIAL 

Figure  3.  Differential  soil  elements 


such  that  they  encompass  a  unit  volume  of  solid  particles,  then 


da  =  1  +  e 

(7) 

0 

dC  =  1  +  e 

(8) 

and 
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dz  =  1 


where  e  is  the  initial  void  ratio  and  e  is  the  void  ratio  at  some 
o 

later  time  during  consolidation.  By  simple  ratios 

—  =  — - -  (1C 

da  1  +  e  ( 


41  =  l  +  e 
dz 


d£  =  l+e 
da  l  +  e 

o 

Thus  conversion  from  one  coordinate  system  to  another  can  be  accom¬ 
plished  by  simple  integration  such  that 


,  f _ 

11  +  e(a,( 


C  =  /  [1  +  e (z , t)  ]  dz 


These  relationships  will  be  used  extensively  throughout  the  remainder  of 
this  development  so  that  equilibrium  and  continuity  conditions  may  be 
expressed  in  the  most  easily  understood  manner  and  then  transformed  into 
the  material  coordinate  system  for  the  governing  equation. 

Material  Equilibrium 

14.  The  equilibrium  of  a  soil  element  having  unit  area 


perpendicular  to  the  page  and  a  unit  volume  of  solid  particles  is  illus¬ 
trated  in  Figure  4.  The  weight,  W  ,  of  the  element  is  the  sum  of  the 
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Figure  4.  Soil  element  in  equilibrium 


weights  of  pore  fluid  and  solid  particles: 


W  =  e  y  +  (1)  Y 

w  s 


(15) 


Therefore,  equilibrium  of  the  soil  mixture  is  given  by 

a  +  |f  dC  +  (e  yw  +  Ys)  -  o  =  0 


(16) 


where  a  is  the  total  stress.  By  simplifying  and  applying  Equation  8, 
an  equation  relating  the  spatial  rate  of  change  in  total  stress  to  the 
void  ratio  and  unit  weights  of  solids  and  fluid  is  obtained: 


e  y  +  Y 

w 

1  +  e 


s 


=  0 


(17) 


Multiplying  through  by  and  substituting  Equation  11  gives  the 

equilibrium  equation  in  terms  of  material  coordinates: 


e 


+  Ys 


=  0 


(18) 


15.  It  is  also  necessary  to  derive  an  expression  for  the  equilib¬ 
rium  of  the  pore  fluid  alone.  Considering  the  total  fluid  pressure  at 
any  time  to  be  composed  of  both  a  static  and  excess  pressure  gives 
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u  =  u  +  u 
w  o 


(19) 


where  uw  >  u  ,  and  u  are  total»  static,  and  excess  pressures, 
respectively.  Static  pressure  equilibrium  is  ensured  if 


3u 

if  +  "  0  (20) 


Therefore,  differentiation  of  Equation  19  yields 


8u 
_ w 

3£ 


=  0 


(21) 


or  in  terms  of  the  material  coordinate 


3u  . 

w  3u  ,  . 

— - t—  +  y  (1  +  e) 

8z  8z  'w 


=  0 


(22) 


Fluid  Continuity 


16.  To  determine  the  equation  of  continuity  for  the  fluid  phase 
of  the  differential  soil  element,  the  weight  of  fluid  inflow  minus  the 
weight  of  fluid  outflow  is  equated  to  the  time  rate  of  change  of  weight 
of  fluid  stored  in  the  element.  As  shown  in  Figure  5,  the  weight  of 
fluid  flowing  into  the  volume  is 


n  •  v  •  yw  (23) 

per  unit  area  where  n  is  the  volume  porosity  which  is  here  assumed 
also  the  area  porosity  and  v  is  the  velocity  of  flow.  Since  the  soil 
solid  particles  are  also  moving  during  consolidation. 


v 


v 

s 


(24) 


where  subscripts  f  and  s  represent  fluids  and  solids,  respectively. 
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K  vrw 

Figure  5.  Fluid  flow  through  a 
differential  element 

The  weight  of  fluid  outflow  is 

£ 

n  .  v  .  Yw  +  —  (n  .  v  .  Yw)dC  (25) 

By  specifying  the  differential  element  to  have  a  unit  volume  of  solid 
particles,  the  weight  of  fluid  contained  within  the  element  is 

e  (26) 

and  its  time  rate  of  change  is  therefore 

Tt  <■*  V  <27> 

Equating  this  time  rate  of  change  of  the  weight  of  fluid  within  an 
element  to  inflow  minus  outflow  results  in 

k  [»<vf  -  V,)]  d5  +  ff  -  0  (28) 

where  the  fluid  is  assumed  incompressible  and  thus  has  a  constant  unit 
weight  which  is  cancelled  in  the  equation. 

17-  Equation  28  is  the  equation  of  continuity  expressed  in  terms 
of  the  convective  coordinate  system.  Utilizing  the  chain  rule  for 
differentiation,  the  relationship 
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(29) 
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e(vc  -  v  )  =  -  — 
f  s  Y 


r  3u 

L JT*\ 


(1  +  e) 


19.  Now  Equations  18,  30b,  32,  and  35  can  be  united  to  produce  a 
governing  equation.  First,  combine  Equation  30b  and  35  to  eliminate  the 
velocity  terras.  Thus 

It  [-  ywitVtt  (jt  +  \  +  e  0_  +  It =  0  (36) 

Next,  use  Equation  32  to  eliminate  uy  in  Equation  36 

3  F  k  (da  3 o '  \  ,  3e  _  x  o 7  \ 

3z  Yw ( 1  +  e)  \3z  3z  ^w  e^w/J  3t 


and  then  Equation  18  to  eliminate  o  in  Equation  37 


k  (  da^_  \a.3e 

Y  (1  +  e)  \  ys  3z  yu)  3t 

m  W  '  '  _ 


(38a) 


( — - — \  +  i_  [ _ k _  lal"|  +  l£  = 

\1  +  ej  3z  |TW(1  +  e)  3z  I  3t 


(38b) 


Again,  by  the  chain  rule  of  differentiation,  the  relationship 

-il  =  ii  n<n 

3z  de  3z  K  ’ 

can  be  written  and  Equation  38b  thus  becomes  Equation  4: 

(ys  \  d  fk(e)  "1  3e  3  ("  k(e)  do’  3e"|  3e  . 

-  y  dT  [TTlj  aT  TT  [yw(i  +  e)  de  TTJ  +  IT  '  0  (4) 

which  is  the  same  as  the  previous  Equation  4  and  constitutes  the  govern¬ 
ing  equation  of  one-dimensional  consolidation  in  terms  of  the  void 
ratio,  e  ,  and  the  functions  k(e)  and  o’(e)  . 

20.  An  analytical  solution  to  Equation  4  is  not  possible,  but 
once  appropriate  boundary  conditions  are  specified,  its  numerical 
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solution  is  feasible  with  the  aid  of  a  computer.  Of  course,  the  rela¬ 
tionships  between  permeability  and  void  ratio  and  effective  stress  and 
void  ratio  must  also  be  known  or  assumed. 


Boundary  Conditions 


21.  Three  types  of  boundary  conditions  are  possible  for  a  soft 
clay  deposit  undergoing  consolidation.  These  are  shown  in  Figure  6  with 
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Figure  6.  Possible  boundary 
conditions 

possible  combinations  at  the  top  and  bottom  of  the  layer.  The  condition 
of  semipermeable  is  an  addition  to  the  usually  assumed  conditions  of 
either  permeable  or  impermeable.  The  semipermeable  condition  represents 
the  state  when  a  compressible  layer  is  in  contact  with  another  different 
compressible  layer  or  when  a  compressible  layer  is  in  contact  with  an 
incompressible  layer  which  has  neither  the  characteristics  of  a  free- 
draining  layer  nor  those  of  an  impermeable  layer,  but  something  in 
between. 

22.  For  the  case  of  a  free-draining  boundary,  there  is  no  excess 
fluid  pressure  and  the  total  fluid  pressure  is  equal  to  the  static 
pressure 


u 

w 


u 

o 


h  y 
w  w 


(40) 


where  h^  is  the  height  of  the  free  water  table  above  the  boundary. 
Since  the  total  weight  of  material  above  the  boundary  is  known,  total 
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stress  may  be  calculated,  and  by  the  effective  stress  principle,  effec¬ 
tive  stress  can  be  calculated.  The  void  ratio  is  then  deduced  from  the 
known  or  assumed  relationship  between  it  and  effective  stress. 

23.  At  an  impermeable  boundary,  there  is  no  fluid  flow  and  thus 


vf  =  vs  (41) 

Applying  this  to  Equation  35  results  in 

3u 

JT  +  V1  +  e>  =  0  (42> 

but  consideration  of  Equation  32,  the  effective  stress  equation,  gives 

If  -  If  +  V1  +  c>  ■ 0  <43) 

Now  if  Equation  18  is  used  to  replace  the  total  stress  term  and  the 
relationship  of  Equation  39  is  used  to  express  the  effective  stress  part 
in  terms  of  the  void  ratio,  Equation  43  can  be  written 


3e  Ys  ~  Yw 
3z  do ' 
de 


=  0 


(44) 


which  is  the  boundary  condition  where  the  compressible  layer  meets  an 
impermeable  layer. 

24.  The  boundary  condition  for  a  semipermeable  layer  is  based  on 
the  fact  that  the  quantity  of  fluid  flowing  out  of  one  layer  must  equal 
the  quantity  of  fluid  flowing  into  the  layer  across  their  common  bound¬ 
ary.  The  quantity  of  fluid  flowing  across  a  boundary  of  unit  area  is 


n(vf  -  vg) 


(45) 


There fo re 


[n(vf  -  vg)] 

upper 


[n(v  -  v  )  ] 

lower 


(46) 


20 


I 
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where  the  subscripts  indicate  upper  and  lower  layers.  Then  from  Equa¬ 
tion  33  and  the  relationship  of  Equations  29  and  11 


/  k  8u\  _  /  k  3u\ 
\1  +  e  9z \1  +  e  9z 


where  y  is  eliminated  because  the  same  fluid  is  in  both  lasers  and 
w 

1  and  2  indicate  upper  and  lower  layers,  respectively.  It  should 
also  be  noted  that  the  total,  static,  and  therefore  excess  fluid  pres¬ 
sures  must  be  equal  in  the  two  layers  at  their  common  boundary 


(u^  =  (u)  2 


25.  From  the  effective  stress  principle, 

¥■  -  IT*  ■  ¥-  (49) 

9z  3z  9z 

By  use  of  the  equilibrium  conditions  of  Equations  18  and  22,  Equation  49 
can  be  rewritten  as 


9z  ^s  9z 


which  can  also  be  written 


/  9u\  de 

y^w  9z )  da 1 


The  conditions  expressed  by  Equations  47,  48,  and  51  may  be  used  to  al¬ 
low  numerical  solution  to  the  problem  of  semipermeable  boundaries. 

Initial  Conditions 


26.  Initial  conditions  through  a  compressible  layer  will  vary 
according  to  the  stress  history  of  the  layer.  Since  it  is  necessary  to 
solve  the  governing  equation  by  an  approximate  numerical  technique,  any 
initial  distribution  of  void  ratios  is  permissible  so  long  as  it  is 


consistent  with  the  assumed  void  ratio  versus  effective  stress  relation¬ 
ship.  Typical  inital  void  ratio  distributions  in  qualitative  terms  are 
as  follows: 

a_.  A  dredged  fill  layer  will  have  a  high  uniform  initial 
void  ratio  distribution. 

Id.  A  layer  consolidated  under  self  weight  only  will  have 
relatively  high  initial  void  ratios  which  decrease 
considerably  with  depth  in  the  layer. 

c_.  A  layer  normally  consolidated  under  a  small  surcharge 
load  will  have  intermediate  void  ratios  which  decrease 
with  depth. 

dl.  A  layer  consolidated  under  a  large  surcharge  load  or 
overconsolidated  will  have  relatively  low  initial  void 
ratios  which  decrease  only  slightly  with  depth. 

The  value  of  these  void  ratios  and  their  exact  distribution  will  depend 
on  the  void  ratio-effective  stress  relationship  chosen  and  any  existing 
surcharge . 


PART  III:  SOLUTION  OF  THE  GOVERNING  EQUATION 

27.  An  analytical  solution  of  the  one-dimensional  finite  strain 
governing  equation  is  not  possible  because  of  the  nonlinear  nature  of 
its  coefficients.  However,  a  numerical  solution  of  the  equation  is 
feasible  if  these  coefficients  are  constantly  updated  during  the  solu¬ 
tion  to  simulate  their  nonlinearity.  An  explicit  finite  difference 
scheme  has  been  chosen  to  solve  the  equation  because  of  its  relatively 
simple  algorithm,  but  this  scheme  does  necessitate  stringent  stability 
criteria  which  will  be  discussed  in  a  later  section. 


Explicit  Finite  Difference  Scheme 


28.  The  finite  difference  procedure  is  a  method  of  representing 
a  differential  term  by  means  of  finite  differences.  Time  space  is 
broken  down  into  intervals  of  finite  length  denoted  r  .  The  time 
derivative  of  void  ratio  can  then  be  written 


3e 

at 


i,j+l 


) 


(52) 


where  the  subscripted  terms  are  as  shown  in  Figure  7 .  If  the  space 
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Figure  7.  Finite  difference  mesh 
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coordinate  is  divided  into  intervals  denoted  6  ,  the  derivative  of  void 
ratio  with  respect  to  space  is 


H  <*i-V  *  16  <ei+ltJ  -  W 


(53) 


by  the  central  difference  method,  and  the  second  derivative  of  void 
ratio  with  respect  to  space  is 


y  (vv  *7  -  2ei,j + 


(54) 


where  terms  are  also  as  shown  in  Figure  7 . 


Simulation  of  Nonlinearity 


29.  It  is  appropriate  here  to  rewrite  the  general  governing 
equation  (Equation  4)  in  the  form 


|Yce(e)  +  [a(e)]J  ff  +  <*(e)  -jj-§  +  Yw  |f  =  0 


(55) 


where 


and 


a(e) 


k(e)  do  * 
1  +  e  de 


(56) 

(57) 


(58) 


To  simulate  the  equation  nonlinearity,  the  functions  a(e)  and  8(e) 
are  recalculated  at  each  time  step  for  the  current  value  of  the  void 
ratio  at  each  point  in  the  z  space  grid. 

30.  In  the  computer  program  developed  for  this  report,  point  data 
are  input  relating  void  ratio  to  permeability  and  effective  stress 


similar  to  that  which  would  be  obtained  from  laboratory  testing.  To  en¬ 
sure  smooth  continuous  functions,  however,  additional  points  are  inserted 
between  the  laboratory  determined  points.  A  typical  trace  of  such  data 
is  shown  in  Figure  8.  Using  these  data,  tables  of  values  for  a(e)  and 
8(e)  at  various  values  of  e  can  be  constructed  by  numerical  differ¬ 
entiation.  Then,  by  a  linear  interpolation,  the  value  of  a(e)  and 
8(e)  for  any  value  of  e  can  be  obtained. 


x  ADDED  POINTS 

Figure  8.  Typical  plot  relating  void  ratio,  e  , 
to  permeability,  k  ,  and  effective  stress,  o' 

31.  The  solution  to  the  governing  equation  in  finite  differences 
can  now  be  written 
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From  Equation  59  it  is  seen  that  the  void  ratio  along  point  z^  at  a 
future  time,  tj+x  >  is  explicitly  determined  from  the  values  of  the 
void  ratio  at  that  point  and  its  nearest  neighbors  at  time  t^  and  func 
tions  of  the  void  ratio  at  these  same  points  at  the  present  time,  tj  . 
Thus,  once  initial  and  boundary  conditions  are  determined,  the  consolida 


tion  problem  is  solved. 


Solution  for  Initial  Conditions 


32.  Calculation  of  the  initial  void  ratio  distribution  in  a  com¬ 
pressible  layer  is  dependent  on  the  unit  weights  of  solids  and  fluid  in 
the  layer,  the  effective  weight  of  any  existing  surcharge,  and  the  rela¬ 
tionship  between  void  ratio  and  effective  stress  within  the  layer.  To 
illustrate  the  procedure,  assume  the  compressible  and  saturated  layer 
shown  in  Figure  9  is  fully  consolidated  under  its  own  self  weight  only 


Figure  9.  Initial  void  ratio  distribution  in 
a  compressible  layer  consolidated  under  self 
weight  only 


before  a  surcharge,  Aq  ,  is  added  which  will  cause  further  consolida¬ 
tion.  The  initial  conditions  in  the  layer  at  t  =  0+  are  then  the  same 
as  conditions  in  the  layer  at  t  =  0  assuming  the  surcharge  is  quickly 
added  at  t  =  0  .  This  is  so  because  the  fluid  in  the  layer  has  not  had 
time  to  drain,  and  therefore,  initially,  fluid  pressure  will  support  all 
the  added  surcharge.  Of  course,  as  time  goes  by  the  surcharge  load  will 
gradually  be  transferred  to  the  t;pil  particles  causing  the  solid  skele¬ 
ton  to  compress. 
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33.  To  determine  the  initial  void  ratio  distribution, 
the  equation 


da 

1  +  e  (a) 
o 


l 


eo(a>  * 


(60) 


must  be  solved  where  h  is  the  initial  layer  height  in  Lagrangian  co¬ 
ordinates  and  J l  is  the  initial  height  in  material  coordinates.  Since 
there  are  two  unknowns  in  this  equation,  it  cannot  be  solved  without  some 
additional  information.  In  a  fully  consolidated  state,  the  effective 
stress  distribution  through  a  layer  depends  only  on  the  buoyant  weight 
of  solids  and  any  existing  surcharge  such  that 

l 

a'(z,o)  =  /  (y  -  y  )dz  +  q  (61) 

f  s  w  o 

z 

When  Equations  60  and  61  are  used  in  conjunction  with  the  relationship 
between  void  ratio  and  effective  stress  such  as  that  shown  in  Figure  8, 
the  number  of  relationships  matches  the  number  of  unknowns  and  solution 
is  possible.  However,  even  if  the  relationship  between  void  ratio  and 
effective  stress  were  expressed  analytically  and  the  appropriate  sub¬ 
stitutions  made  in  Equations  60  and  61,  a  transcendental  equation  would 
result  which  would  require  an  iterative  type  solution.  Therefore,  an 
incremental  technique  will  be  used  here  which  will  approach  the  exact 
solution  from  the  lower  side. 

34.  It  is  first  necessary  to  divide  the  compressible  layer  into 
a  number  of  elemental  layers  of  length 

Aa  =  |  (62) 

where  N  is  any  positive  integer.  The  larger  the  N  ,  the  more  accu¬ 
rate  the  solution.  The  uppermost  elemental  layer  is  subject  to  an  ef¬ 
fective  stress  equal  to  the  effective  weight  of  any  existing  surcharge, 
q  .  When  this  effective  weight  is  used,  a  void  ratio  is  obtained  from 
data  such  as  Figure  8.  This  void  ratio  is  assumed  constant  for  the 
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elemental  layer.  Therefore,  for  the  first  layer 


Az(l) 


Aa 

1  +  eo(l) 


(63) 


and 


Aa'(l)  =  (y  -  Y  ) Az(l)  +  q 
s  w  o 


(64) 


When  Ao'(l)  is  used  as  the  effective  stress  acting  on  the  second  in¬ 
cremental  element,  the  void  ratio  of  the  second  element  can  be  deter¬ 
mined.  Following  this  technique  throughout  the  entire  layer  results  in 
the  initial  void  ratio  distribution  sufficiently  accurate  for  computation 
of  future  consolidation. 

35.  For  the  case  of  a  dredged  fill,  it  is  assumed  that  the  layer 
is  deposited  at  a  uniform  consistency,  and  after  initial  solids  sedi¬ 
mentation  the  compressible  layer  exists  at  a  uniform  void  ratio  with 
zero  effective  stress  throughout  the  layer.  Under  these  conditions, 
total  layer  height  in  material  coordinates  is  calculated  directly  from 


l 


h 

1  +  e 

o 


(65) 


where  h  is  the  height  of  the  compressible  layer  after  initial  sedi¬ 
mentation  but  before  any  consolidation. 


Void  Ratio  at  Boundaries 


36.  Void  ratio  calculation  at  a  free-draining  boundary  is  actually 
a  calculation  of  effective  stress  at  the  boundary.  This  calculation  is 
done  through  a  knowledge  of  the  total  weight  of  materials  above  the 
boundary  plus  any  existing  or  added  surcharge  and  the  distance  of  the 
boundary  below  the  free  water  surface.  Since  there  is  no  excess  fluid 
pressure,  the  effective  stress  is 

o’  =  a  +  Aq  -  y  h  (66) 

w  w 


28 


where  o  is  the  total  stress  due  to  any  existing  surcharge  and  material 


self  weight,  Aq  is  an  added  surcharge,  and  h  is  distance  of  the 

w 


boundary  below  the  free  water  surface.  With  this  effective  stress,  the 
persistent  void  ratio  can  then  be  determined  from  a  relationship  such  as 
shown  in  Figure  8. 

37.  The  determination  of  void  ratio  at  an  impermeable  boundary 
requires  the  use  of  a  fictitious  mesh  point  outside  the  boundary  as 
shown  in  Figure  10.  Using  the  initial  void  ratio  distribution  <'r 
distribution  at  any  time,  t^  ,  the  void  ratio  at  this  fictitious  mesh 
point  is  calculated  by  expressing  Equation  44  in  finite  difference 
terms.  Thus 


e.  .  =  e„  .  +  26 
0,J  2,j 


(fr) 


(y. 


V 


(67) 


j 


where 


de 

do' 


is  determined  for  e 


from  data  such  as  in  Figure  8. 


With  e 


0,j 


determined,  e^  is  then  found  from  Equation  59  and  the 


whole  process  repeated  at  each  time  step. 


Figure  10.  Void  ratio  calculation  at  an 
impermeable  boundary 


38.  When  a  compressible  layer  lower  boundary  is  neither  free 
draining  nor  impermeable,  void  ratio  calculation  at  the  boundary  is 
accomplished  by  writing  a  finite  difference  expression  for  Equation  51 
and  using  an  imaginary  mesh  point  as  was  done  for  the  impermeable  case. 


(68) 


e 


0,  j 


2,j 


where  the  term  -r—  is  either  calculated  from  the  previous  time  step  or 

O  Z 

assumed.  In  the  case  of  a  dredged  fill  overlying  a  compressible  layer, 
the  excess  pressure  gradient  at  the  layer  interface  is  assumed  to  be 
zero  for  the  first  time  step  and  thereafter  it  is  calculated  based  on 
the  previous  conditions  and  Equations  47  and  48.  The  procedure  is  shown 
schematically  in  Figure  11.  The  method  of  calculating  excess  pressure 


Figure  11.  Schematic  representation  of  void  ratio 
calculation  at  semipermeable  boundaries 
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from  void  ratio  and  vice  versa  is  given  in  a  later  section.  The  void 
ratio  of  the  top  point  in  the  compressible  lower  layer  is  based  on 
Equation  48  and  the  fact  that  the  change  in  excess  pore  pressure  equals 
the  negative  change  in  effective  stress.  In  the  case  of  a  compressible 
layer  overlying  a  semipermeable  incompressible  layer,  the  permeability, 
void  ratio,  and  a  typical  drainage  path  length  in  the  incompressible 
layer  must  be  either  measured  or  assumed.  The  calculation  procedure  is 
also  illustrated  in  Figure  11.  Only  a  typical  illustration  of  marching 
forward  in  time  is  shown,  but  this  holds  for  all  void  ratios  except  at 
the  imaginary  points  and  the  top  point  in  a  compressible  foundation 
layer. 


Settlement  Calculation 


39.  The  calculation  of  settlement  at  any  point  in  a  compressible 
layer  is  simply  the  subtraction  of  its  convective  coordinate  from  its 
Lagrangian,  or  initial,  coordinate.  If  settlement  at  a  point  is  denoted 
S(z,t)  ,  then 


S(z,t)  =  a(z,0)  -  £(z,t) 


(69) 


and  by  integration  of  Equations  10  and  11 


S(z ,t) 


+  e(z ,0) ]dz 


+  e(z, t) ]dz 


(70) 


Since  data  are  generated  aroun  mesh  points  in  the  finite  difference 
solution  of  the  consolidation  problem,  the  numerical  integration  of 
Equation  70  by  Simpson's  rule  is  a  simple  exercise. 

40.  A  common  method  of  expressing  the  state  of  consolidation  in 
small  strain  theories  is  by  the  percentage  of  excess  pore  pressure  dis¬ 
sipated.  In  the  finite  strain  theory,  degree  of  consolidation  is  appro¬ 
priately  defined  as  the  ratio  of  current  settlement  to  final  settlement 
in  the  entire  layer.  Thus 
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(71) 


r 


S(t,t) 

ut  s(*,») 

where  S ( il , °°)  Is  the  ultimate  settlement  of  the  layer  when  all  excess 
pore  pressure  has  dissipated. 


Calculation  of  Stresses  and  Pressures 

41.  Once  the  void  ratio  distribution  throughout  a  compressible 
layer  is  determined,  the  distribution  of  effective  stress  can  be  ob¬ 
tained  from  a  relationship  such  as  shown  in  Figure  8.  The  static  pore 
pressure  is  also  immediately  determined  for  each  mesh  point  as 


u  (z,t)  =  y  [h  -  C(z,t)] 

O  w  i. 


where  h^  is  the  height  of  the  free  water  surface  above  the  datum  plane, 
z  =  0  ,  and  £  is  the  convective  coordinate  of  the  mesh  point  at  the 
time  in  question. 

42.  The  total  stress  at  a  point  in  the  compressible  layer  is 
equal  to  the  total  weights  in  a  unit  area  of  all  materials  above  it  plus 
any  surcharge.  Thus 


a(z,t)  =  Yw  h2  +  /  e(z,t)dz  +  yl  /  dz  I  +  qQ 


where  h2  is  the  height  of  the  free  water  surface  above  the  top  (z  =  £) 
of  the  compressible  layer,  the  integrals  represent  the  volumes  of  fluid 
and  solids  in  the  compressible  layer,  respectively,  and  qQ  is  any 
surcharge. 

43.  With  total  and  effective  stresses  determined,  the  effective 
stress  principle  is  used  to  calculate  total  pore  pressure 


u  ( z , t )  =  o(z,t)  -  o'(z,t) 

W 


and  excess  pore  pressure  is  the  difference  between  total  and  static 
pressures , 


u(z,t)  =  u  (z,t)  -  u  (z,t) 
w  o 


Solution  Consistency,  Convergence,  and  Stability 

44.  Now  that  a  solution  technique  for  solving  the  finite  strain 
consolidation  problem  has  been  formulated,  some  assurance  that  this  tech¬ 
nique  gives  a  correct  answer  is  necessary.  Consistency  implies  that  the 
difference  equations  actually  do  approximate  the  differential  equation. 
Convergence  means  that  the  numerical  solution  is  a  close  approximation 
of  the  exact  solution.  Stability  implies  that  small  errors  introduced 
initially  or  at  a  boundary  remain  bounded  as  the  computations  progress. 
Keller  (1960)  has  shown  that  for  a  parabolic  partial  differential  equa¬ 
tion  of  the  form 

2 

4f  ~  a(z,t)  -  2b(z ,  t)  +  c(z  ,t)e  =  d(x,t)  (76) 

dt  „  Z  dZ 

dZ 

consistency,  convergence,  and  stability  are  assured  in  an  explicit 
finite  difference  scheme  if 


a(z,t) 
I  b(z, t) 


6 


where  6  and  t  are  the  spatial  and  time  mesh  spacings,  respectively, 
and  a  ,  b  ,  and  c  are  any  variables. 

45.  In  the  governing  Equation  55  for  finite  strain  consolidation 
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a(z,t)  =  -  aM  (79) 

Yw 

b(z,t)  =  -  ^“{yc  B(e)  +\l  [o(e)]}  (80) 

w 

c(z,t)  =  0  (81) 

d(z,t)  =  0  (82) 


where  a(e)  and  8(e)  are  as  previously  defined  in  Equations  58 
and  57,  respectively.  Therefore,  if 


6  < 


_ _ 2ot(e) _ 

7c  B(e)  +  h  Ia(e)] 


(833 


and 


x  < 


&2  X 


2a(e) 


(84) 


then  the  solution  should  be  consistent,  convergent,  and  stable.  To 
ensure  these  criteria  are  met  throughout  the  solution  process,  Equa¬ 
tions  83  and  84  should  be  periodically  checked  using  the  extreme  values 
of  a(e)  and  8(e)  to  be  expected  in  the  problem. 
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PART  IV:  SOIL  PARAMETERS  FOR  FINITE  STRAIN  CONSOLIDATION 


46.  Calculation  of  the  consolidation  of  soft  deposits  by  finite 
strain  theory  requires  the  determination  of  the  specific  gravity  of 
solids  in  the  compressible  layer,  the  relationship  between  void  ratio 
and  effective  stress,  and  the  relationship  between  void  ratio  and  perme¬ 
ability.  These  determinations  are  presently  routine  laboratory  proce¬ 
dures  for  fine-grained  soils  normally  encountered  in  earth  construction. 
The  use  of  standard  oedometer  tests  for  soft  deposits  which  may  be  under¬ 
consolidated  in  situ  involves  uncertainties;  for  instance,  a  thin  oedom¬ 
eter  sample  with  no  excess  pore  pressure  and  subjected  to  a  sudden  load 
increment  may  not  react  in  the  same  way  as  an  underconsolidated  thick 
sample  whose  excess  pore  pressure  is  slowly  decreased.  Additionally, 

the  consolidation  induced  by  the  hydraulic  gradient  of  a  permeability 
test  may  not  be  adequately  accounted  for  in  the  test  results.  The 
answers  to  these  questions  are  beyond  the  scope  of  this  report  and  need 
research  to  either  relate  soft  deposit  parameters  to  the  results  of 
conventional  tests  or  devise  new  test  methods  so  that  direct  measure¬ 
ments  can  be  made. 

47.  In  order  to  demonstrate  the  use  of  the  computer  program 
CSLFS,  the  soil  parameters  necessary  were  deduced  from  conventional 
oedometer  test  data  such  as  may  be  generated  in  any  well  equipped  soils 
testing  facility.  By  logical  extrapolation  of  these  data  generated  by 
the  oedometer  testing  over  the  full  range  of  void  ratios  that  might  be 
encountered,  reasonable  solutions  to  the  dredged  fill  consolidation  prob¬ 
lem  can  be  obtained.  Of  course,  the  test  results  on  a  thick  normally 
consolidated  or  overconsolidated  soil  under  a  surcharge  should  be  di¬ 
rectly  applicable  without  extrapolation. 

48.  Use  of  the  program  feature  enabling  the  specification  of 
boundary  conditions  that  are  neither  free  draining  nor  impermeable  re¬ 
quires  that  a  void  ratio,  permeability,  and  drainage  path  length  for 
the  incompressible  foundation  material  be  given.  While  it  is  generally 
possible  to  determine  void  ratio  and  permeability  by  laboratory  testing 
on  undisturbed  samples,  the  distance  required  for  dissipation  of  excess 


pore  pressures  in  the  incompressible  foundation  must  be  estimated  based 
on  engineering  judgment. 

Void  Ratio-Effective  Stress  Relationship 

49.  The  conventional  laboratory  oedometer  test  can  be  used  to 
establish  the  void  ratio-effective  stress  relationship  required  for  cal¬ 
culation  of  consolidation  by  finite  strain  theory  subject  to  the  uncer¬ 
tainties  previously  raised.  Principally,  the  only  difference  between 
testing  soft  deposits  and  the  stiffer  soils  usually  tested  is  in  the  size 
of  the  load  increments  used.  For  routine  tests  of  most  soils,  the  load¬ 
ing  schedule  starts  at  0.25  tsf*  and  is  doubled  for  each  succeeding  in¬ 
crement  until  a  total  load  of  16.0  tsf  is  applied.  Typical  tests  of 
soft  deposits  such  as  channel  sediments  or  dredged  fill  start  at 

0.012  tsf  and  are  incrementally  increased  to  1.0  tsf.  At  these  extremely 
low  pressures,  accurate  account  must  be  taken  of  the  weights  of  load 
transfer  hardware  and  even  the  force  exerted  by  dial  gage  springs 
(Palermo,  Montgomery,  and  Poindexter  1978). 

50.  Perhaps  the  best  method  of  graining  insight  into  the  behavior 
of  soft  clay  soils  is  to  examine  some  typical  oedometer  test  results. 

In  Figures  12  and  13  are  plotted  e-log  o'  curves  as  determined  in  the 
Soils  Testing  Facility  at  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station.  These  plots  have  been  corrected  from  the  originally  reported 
results  (Palermo,  Shields,  and  Hayes  in  press)  by  assuming  100  percent 
saturation  at  test  completion.  This  was  necessary  because  direct 
measurements  of  the  specific  gravity  of  soil  solids  were  not  made  and 
original  results  consistently  indicated  saturation  greater  than  100  per¬ 
cent  when  average  specific  gravity  values  were  assumed. 

51.  Figure  12  shows  four  samples  taken  from  the  Craney  Island 
dredged  material  disposal  site,  one  sample  of  channel  sediments  con¬ 
sidered  typical  of  what  goes  into  the  disposal  area,  and  one  sample  of 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  of  measurement  is  found  on  page  5. 
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Figure  12.  Oedometer  test  results  for  Craney  Island  samples 
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Figure  13.  Oedometer  test  results  for  other  samples 
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the  foundation  soil  beneath  the  disposal  area.  As  can  be  seen  from  the 
figure,  these  soft  deposits  generally  have  characteristics  similar  to 
other  soils  encountered  in  construction  practice  except  that  the  range 
of  void  ratios  these  deposits  undergo  during  consolidation  is  much 
greater.  The  tendency  for  initial  void  ratios  to  increase  as  in  situ 
confining  stresses  decrease  is  also  apparent  from  the  figure.  A  conven¬ 
tional  analysis  to  determine  the  preconsolidation  pressure  from  the 
e-log  o'  curves  is  probably  not  appropriate  since  there  is  no  way  to 
obtain  a  truly  undisturbed  sample  of  such  soft  soils.  However,  the  nor¬ 
mally  consolidated  portion  of  the  curves  should  be  a  valid  indication  of 
the  soil  behavior  as  indicated  by  the  fact  that  all  dredged  material 
curves  including  that  for  channel  sediments  are  approximately  parallel 
over  their  normally  consolidated  range. 

52.  Consolidation  characteristics  of  other  soft  materials  are 
shown  in  Figure  13.  Here  again,  the  extremely  wide  variation  in  void 
ratios  over  relatively  small  stress  ranges  should  be  noted.  The  unusual 
upturn  in  these  curves  at  the  low  end  of  the  stress  range  may  be  pecu¬ 
liar  to  the  particular  test  procedure  or  may  be  valid  indicators  of  the 
behavior  of  these  materials.  Definite  conclusions  cannot  be  drawn  with¬ 
out  further  testing. 

53.  To  illustrate  the  method  of  obtaining  the  necessary  void 
ratio-effective  stress  relationship  for  use  in  the  computer  program  CSLFS, 
consider  the  data  points  as  shown  in  Figure  12.  It  is  proposed  that 
those  points  defining  the  normally  consolidated  portions  of  the  e-log  o' 
curves  fully  describe  the  material  behavior  between  effective  stresses 

of  about  0.01  tsf  to  1.0  tsf.  Defining  the  curve  below  and  above  these 
values  Is  a  matter  of  judgment  in  the  absence  of  experimental  evidence 
dictating  otherwise.  The  arbitrary  extension  of  the  normally  consoli¬ 
dated  portion  in  a  straight  line  is  unreasonable  since  this  would  give 
an  infinite  void  ratio  at  zero  effective  stress  and  a  zero  void  ratio 
at  some  finite  effective  stress.  Probably  a  more  reasonable  assumption 
is  that  there  will  be  some  finite  void  ratio  at  zero  effective  stress 
and  that  the  curve  will  become  asymptotic  to  some  minimum  void  ratio 
depending  on  the  origins  of  the  soil.  It  is  therefore  further  proposed 
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that  the  void  ratio  at  zero  effective  stress  be  selected  as  somewhere 
between  the  void  ratio  at  the  intersection  of  the  normally  consolidated 
line  with  the  effective  stress  ordinate  0.001  tsf  and  the  measured  void 
ratio  before  oedometer  testing.  The  curve  at  effective  stresses  higher 
than  1.0  tsf  should  ideally  be  based  on  oedometer  testing  at  these 
higher  stresses,  but  in  the  absence  of  such  data  may  reasonably  be  an 
extension  of  the  normally  consolidated  portion  which  is  brought  asymp¬ 
totic  to  a  constant  void  ratio  value  between  0.4  and  0.7.  Figure  14 
shows  such  curves  constructed  from  the  data  of  Figure  12.  Void  ratios 
of  7.0  for  the  dredged  fill  and  3.0  for  the  foundation  soil  at  zero 
effective  stress  were  chosen  as  about  midway  between  the  previously  pro¬ 
posed  range  of  possibilities. 

54.  Before  the  final  decision  is  made  to  use  such  a  void  ratio- 
effective  stress  relationship  in  the  computer  program  CSLFS,  the  curve 
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Figure  14.  Void  ratio-effective  stress  relationships  for  soft  dredged 
fill  and  foundation  materials  at  Craney  Island 
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should  be  replotted  on  an  arithmetic  scale  to  ensure  the  curve  is  a 
smooth  continuous  function  without  extraneous  reverse  curvature  and  with 
continuous  derivatives.  Figure  15  shows  such  a  plot  for  the  dredged 
fill  material,  and  Figure  16  shows  the  plot  for  the  foundation  soil  from 
the  Craney  Island  site.  The  points  shown  on  the  plots  are  the  points 
to  be  used  as  program  input  for  a  practical  example  to  be  worked. 


Void  Ratio-Permeability  Relationship 


55.  The  determination  of  the  void  ratio-permeability  relationship 
necessary  for  calculation  of  consolidation  by  the  computer  program  CSLFS 
will  also  be  accomplished  through  use  of  oedometer  test  results.  Be¬ 
cause  conventional  oedometer  testing  involves  relatively  thin  samples 
and  relatively  small  load  increments,  analysis  of  this  testing  based  on 
the  assumptions  of  small  strain  consolidation  theory  will  probably  pro¬ 
duce  sufficiently  accurate  values  of  permeability. 

56.  By  small  strain  theory,  a  nondimens ional  time  factor  is  de¬ 
fined  by 


T  = 


(85) 


where  t  is  real  time,  H  is  the  drainage  path  length,  and  the  coeffi¬ 
cient  of  consolidation,  c^  ,  is 


c 

v 


k(l  +  e) 

Y  a 
w  v 


(86) 


where  k  is  permeability,  e  is  void  ratio,  and  y^  is  unit  weight  of 
water  as  previously  defined.  The  coefficient  of  compressibility,  , 
is  defined  as 


a 

v 


Ae 

Aa' 


(87) 


where  Ae  is  the  change  in  void  ratio  corresponding  to  the  change  in 
effective  stress,  Ao'  .  Combining  the  three  preceding  equations 
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results  in  an  expression  for  permeability 


T  y  Ae  H2 

w _ 

(1  +  e)  t  Ao' 


(88) 


which  involves  known  or  measurable  quantities  in  the  oedometer  test. 

57.  Typically,  consolidation  time  curves  for  each  load  increment 
are  used  to  determine  the  time,  t  ,  for  50  percent  consolidation  where 
analytically  T  =  0.197  for  an  initial  uniform  distribution  of  excess 
pore  water  pressure.  The  void  ratio,  e  ,  is  also  determined  at  t^^ 
from  a  knowledge  of  the  specific  gravity  of  solids,  total  weight  of 
solids,  and  current  sample  volume.  The  drainage  path  length,  H  ,  is 
estimated  as  one-half  the  sample  height  at  t^^  .  An  average  coeffi¬ 
cient  of  compressibility  is  obtained  by  dividing  the  total  void  ratio 
change  during  the  load  increment  by  the  load  increment. 

58.  Permeabilities  determined  in  this  manner  for  the  foundation 
soil  and  dredged  fill  of  the  Craney  Island  disposal  site  are  shown  in 
Figure  17.  While  the  data  at  the  higher  void  ratios  is  considerably 
scattered,  the  data  in  the  lower  void  ratios  which  is  less  scattered 
does  seem  to  give  a  good  fit  when  extended.  Here  again,  the  behavior 
of  the  void  ratio-permeability  relationship  outside  the  range  of  data 
points  is  purely  speculative  until  such  time  as  adequate  testing  is  de¬ 
vised  and  used  in  defining  the  curve  over  the  full  range  of  possible 
void  ratios.  However,  it  is  probably  reasonable  to  assume  that  permea¬ 
bility  becomes  infinitesimally  small  at  some  finite  void  ratio  and  thus 
the  curve  will  become  asymptotic  to  this  void  ratio. 

59.  Figure  18  shows  the  relationship  between  void  ratio  and  per¬ 
meability  for  the  same  other  samples  of  soft  deposits  described  pre¬ 
viously  in  Figure  13.  The  behavior  of  these  curves  at  the  higher  void 
ratios  may  be  an  idiosyncrasy  of  the  test  procedure  since  it  is  probably 
more  reasonable  to  expect  that  permeability  would  increase  more  dramat¬ 
ically  as  the  void  ratio  reached  some  maximum  limit  where  the  soil  no 
longer  forms  an  interconnected  network  of  solid  particles. 

60.  As  before,  it  is  beneficial  to  plot  the  void  ratio-permeability 
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Figure  17.  Void  ratio-permeability  relationships  for  soft  dredged  fill 
and  foundation  materials  at  Craney  Island 


relationship  on  an  arithmetic  scale  as  an  aid  in  determining  the  point 
data  for  use  in  the  program  CSLFS.  Figure  19  shows  such  a  plot  for  the 
dredged  fill  material,  and  Figure  20  is  of  the  foundation  soil  at  the 
Craney  Island  Site.  The  points  shown  on  the  figures  are  the  points  to 
be  used  as  program  input  for  a  practical  example. 


Semipermeable  Boundary  Parameters 

61.  As  previously  shown,  the  boundary  conditions  between  two 
compressible  layers  undergoing  consolidation  are  automatically  determined 
by  the  program  CSLFS  based  on  the  continuity  of  fluid  flow  and  current 
void  ratio  and  permeability  conditions  in  the  compressible  layers.  Where 
a  compressible  layer  bounds  an  incompressible  layer,  boundary  condi¬ 
tions  are  determined  by  the  program  based  on  current  conditions  in  the 
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Figure  18.  Void  ratio-permeability  relationship  for  other  samples 
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Figure  20.  Void  ratio-permeability  relationship  for  foundation  soil  to 
be  used  in  computer  program  CSLFS 

compressible  layer  and  specified  void  ratio,  permeability,  and  length 
of  drainage  path  for  the  incompressible  layer.  It  was  also  previously 
stated  that  the  void  ratio  and  permeability  for  the  incompressible  layer 
should  generally  be  determined  by  laboratory  testing  on  undisturbed 
samples  and  that  specification  of  the  drainage  path  length  is  a  matter 
of  engineering  judgment.  The  basis  for  making  such  a  judgment  is  dis¬ 
cussed  in  this  section. 

62.  The  drainage  path  length  is  defined  as  that  distance  required 
for  complete  dissipation  of  excess  pore  water  pressure  existing  at  the 
layer  boundary.  Together  with  this  pore  pressure,  it  is  used  to  deter¬ 
mine  the  excess  pressure  gradient  at  the  incompressible  layer  side  of 
the  boundary  by  the  equation 
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(89) 
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where  u  is  the  excess  pore  pressure  at  the  boundary,  x  is  the  drain¬ 
age  path  length  measured  in  the  Lagrangian  coordinate  system,  and  e  is 
the  void  ratio  of  the  incompressible  layer.  The  excess  pore  pressure 
is  calculated  as  previously  described  from  the  void  ratio  of  the  compres¬ 
sible  layer.  The  pressure  gradient  thus  obtained  is  used  in  Equation  47 
to  determine  the  excess  pore  pressure  gradient  on  the  compressible  layer 
side  of  the  boundary  as 
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where  the  subscripts  comp  and  incomp  refer  to  the  compressible  layer 
and  the  incompressible  layer,  respectively.  This  value  is  then  used  in 
Equation  68  for  computing  the  void  ratio  of  an  image  point  which  enables 
the  computation  of  the  void  ratio  at  the  first  mesh  point  in  the  compres¬ 
sible  layer  at  the  next  time  step. 

63.  An  examination  of  Equation  89  shows  that  if  the  drainage  path 

length  is  chosen  to  be  very  large,  the  effect  is  to  make  very  small 

and  in  the  limit  will  approach  zero  or  the  impermeable  boundary  condi¬ 
tion  which  makes  Equation  68  the  same  as  Equation  67.  At  the  other 
extreme,  if  the  drainage  path  length  is  chosen  to  be  very  small,  the 
effect  is  to  make  —  very  large  and  in  the  limit  will  approach  an 

oZ 

infinite  value.  The  computation  in  Equation  68  then  has  no  physical 
meaning,  but  the  effect  in  the  program  is  to  cause  the  void  ratio  at 
the  first  mesh  point  in  the  compressible  layer  to  be  set  at  its  final 
value  or  the  free-draining  boundary  condition. 

64.  Between  those  conditions  of  impermeable  and  free  draining, 
it  is  proposed  that  the  drainage  path  length  be  chosen  to  equal  the 
depth  of  the  compressible  layer  where  the  material  of  the  incompressible 
layer  is  the  same  or  essentially  the  same  as  that  of  the  compressible 
layer.  Where  the  material  properties  are  substantially  different,  it 

is  further  proposed  that  the  drainage  path  length  be  chosen  to  be  pro¬ 
portional  to  the  ratios  of  the  permeability  functions  times  the  depth 
of  the  compressible  layer.  In  equation  form,  this  means 
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where  x  is  the  drainage  path  length,  h  is  the  depth  of  compressible 
material,  and  k  and  e  are  average  permeability  and  void  ratio,  re¬ 
spectively,  in  the  respective  layers  near  the  interface  of  the  layers. 


PART  V:  CONSOLIDATION  PROBLEMS 


65.  In  this  Part,  the  capabilities  of  the  computer  program  CSLFS 
will  be  demonstrated  by  solving  some  practical  examples  involving  the 
consolidation  occurring  in  a  dredged  fill  disposal  site  subjected  to 
periodic  deposition  of  soft  channel  sediments  and  the  consolidation  of 

a  thick  soft  layer  subjected  to  an  additional  surcharge  due  to  some  con¬ 
struction  activity  above  it.  Figures  will  be  used  to  show  the  distribu¬ 
tions  of  excess  pore  pressure,  void  ratio,  layer  settlement  versus  time, 
and  percent  consolidation  versus  time.  Whenever  possible,  a  comparison 
between  the  results  computed  by  the  finite  strain  formulation  will  be 
compared  with  those  from  a  small  strain  theory  computation. 

Consolidation  of  Dredged  Fill  on  a 
Compressible  Foundation 

66.  In  this  example,  a  large  disposal  site  has  been  proposed  for 
an  area  of  a  bay  where  foundation  material  is  a  soft  marine  sediment 
currently  about  5  ft  below  mean  sea  level.  Considerations  of  the  area 
available  for  disposal  and  the  volume  and  type  of  material  to  be  dredged 
has  led  to  the  conclusion  that  the  site  must  be  capable  of  holding  mate¬ 
rial  deposited  according  to  the  following  schedule: 

Year  1  through  Year  2,  3  ft/year 
Year  3  through  Year  4,  2  ft/year 
Year  5  through  Year  8,  1  ft/year 

The  total  amount  for  each  year  will  be  deposited  during  the  first  few 
weeks  of  each  year  and  therefore  can  be  considered  to  be  dumped  instan¬ 
taneously  in  the  disposal  area  at  the  beginning  of  each  year.  Figure  21 
shows  the  schedule  graphically.  It  should  also  be  noted  that  the  yearly 
amounts  are  based  on  volumes  after  initial  sedimentation  has  taken  place. 
If  initial  sedimentation  is  not  complete  very  soon  after  each  particular 
dredging  operation,  due  consideration  of  the  nonsedimented  height  of 
each  layer  must  be  taken  into  account  when  calculating  the  necessary 
height  of  confinement  dikes. 


..A  «. 
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Figure  21.  Schedule  of  dredged  material  deposition 


67.  The  consolidation  behavior  of  these  dredged  fill  deposits  is 
required  to  be  calculated  in  conjunction  with  the  consolidation  behavior 
of  the  foundation  in  order  that  a  program  of  dike  construction  may  be 
instituted  that  is  neither  overly  conservative  nor  extravagant.  It  is 
further  required  that  an  estimate  be  made  of  the  time  required  for  90  per¬ 
cent  consolidation  of  the  disposal  area  and  ultimate  settlement  so  that 

an  evaluation  of  its  potential  future  use  may  be  made. 

68.  Before  consolidation  can  be  calculated,  laboratory  determina¬ 
tions  must  be  made  of  the  void  ratio-effective  stress  and  void  ratio- 
permeability  relationships  for  both  the  dredged  and  foundation  materials 
along  with  the  unit  weight  of  solids  in  these  materials  and  the  initial 
void  ratio  assumed  by  the  dredged  material  after  initial  sedimentation. 

For  this  example,  the  relationships  depicted  in  Figures  15,  16,  19, 

and  20  will  be  used.  The  dredged  material  is  assumed  to  have  an  initial 
void  ratio  of  7.0  and  a  specific  gravity  of  solids  of  2.75.  The  founda¬ 
tion  is  assumed  to  have  a  specific  gravity  of  solids  of  2.83  and  to  be 
normally  consolidated  under  its  own  weight. 

69.  It  will  be  further  assumed  that  field  borings  were  addition¬ 
ally  used  to  determine  that  the  compressible  foundation  is  20.0  ft  thick 
and  overlays  an  incompressible  layer  of  silty  material  having  an  average 
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void  ratio  of  0.65  and  permeability  of  3.0  *  10  ft/day.  The  void 
ratio  and  permeability  of  the  compressible  foundation  layer  at  the  inter¬ 
face  with  the  incompressible  silt  deposit  could  be  determined  either  by 
field  borings  or  by  assuming  the  layer  is  normally  consolidated  under 
its  own  weight  and  allowing  the  computer  program  to  calculate  its  ini¬ 
tial  conditions.  For  this  example  a  void  ratio  of  1.80  and  permeability 
-4 

of  1.03  x  10  ft/day  have  been  chosen  based  on  program  calculations. 
Equation  91  is  used  to  determine  the  drainage  path  length  for  this  semi- 
permeable  boundary  as  about  6.0  ft. 

70.  The  input  data  required  for  problem  solution  is  shown  in 
Appendix  C  of  this  report.  The  calculation  constants  t  and  6  are 
chosen  small  enough  so  that  problem  detail  and  accuracy  are  preserved, 
yet  large  enough  to  promote  computation  economy.  If  the  constants  are 
too  large  for  the  stability  criteria,  the  program  will  print  an  error 
message.  For  this  problem,  i  =  1.0  day  and  6  is  one-sixth  of  the 
initial  layer  height  for  the  dredged  fill  and  one-tenth  for  the  founda¬ 
tion.  These  selections  proved  sufficient  for  accuracy  and  stability. 

Also  included  in  the  appendix  is  calculated  data  for  the  end  of  the 
second  and  eighth  year  of  consolidation. 

71.  From  these  calculated  data,  a  visual  picture  of  the  consoli¬ 
dation  process  can  be  obtained.  Figure  22  shows  the  void  ratio  distri¬ 
bution  in  the  dredged  fill  at  the  end  of  year  2  after  two  layers  of  fill 
have  been  placed  but  before  the  third  layer  is  placed.  Also  shown  in 
the  figure  are  the  void  ratio  distribution  at  year  1  after  the  second 
layer  is  placed  (which  serves  as  the  initial  conditions  for  the  current 
consolidation  period)  and  the  final  void  ratio  distribution  if  no  more 
dredged  fill  layers  were  to  be  placed.  In  the  figure,  void  ratios  are 
plotted  against  the  material  coordinate,  z  ,  for  ease  in  comparing 
past,  present,  and  future  distributions.  The  conventional  layer  height, 

£  coordinate,  equivalent  to  z  can  be  found  in  the  problem  listing  in 
Appendix  C. 

72.  The  distribution  of  excess  pore  pressure  within  the  dredged 
fill  at  the  end  of  the  second  year  and  before  the  next  layer  is  deposited 
is  shown  in  Figure  23  along  with  the  distribution  at  year  1  after 
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CURRENT  SETTLEMENT:  0.75' 
FINAL  SETTLEMENT!  l.06‘ 

DEGREE  OF  CONSOLIDATION:  71% 
TIME  OF  CONSOLIDATION!  2  YR. 


Figure  22.  Void  ratio  distribution  at 
the  end  of  year  2 
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Figure  23.  Excess  pore  pressure  distribution  at 
the  end  of  year  2 


deposition  of  the  second  layer.  The  discontinuity  in  the  year  1  curve  is 
due  to  the  assumption  that  the  second  layer  is  deposited  instantaneously 
and  its  excess  pore  pressure  is  superimposed  on  the  existing  excess  pres¬ 
sure  before  the  layer  was  deposited.  At  the  end  of  consolidation  there 
is  no  excess  pore  pressure,  and  thus  a  final  curve  is  not  shown.  Curves 
of  this  type  are  useful  in  evaluating  strength  or  stability  using  an  ef¬ 
fective  stress  analysis.  Distributions  of  total  and  effective  stresses 
can  be  found  in  tabular  form  in  the  problem  listing  in  Appendix  C. 

73.  Figure  24  depicts  void  ratio  distributions  throughout  the 
dredged  fill  deposition  period  and  the  final  distribution  for  the  total 
amount  of  material  deposited.  This  figure  shows  that  even  after  100  per¬ 
cent  primary  consolidation,  very  high  void  ratios  will  exist  throughout 
the  dredged  fill  material  and  unless  some  later  load  causing  further 
consolidation  is  placed,  the  material  may  never  be  suitable  for  any 
engineering  purpose.  The  effects  of  surface  desiccation  and  secondary 
consolidation  are  not  considered  here,  even  though  these  factors  will 
have  an  impact  on  the  final  void  ratio  distribution.  The  effects  of 
these  factors  will  be  considered  in  future  extensions  of  the  theoretical 
basis  and  computer  program. 

74.  Shown  in  Figure  25  are  excess  pore  pressure  distributions  in 
the  later  years  of  consolidation.  Again,  this  type  of  figure  would  be 
useful  in  evaluating  strength  or  stability  using  an  effective  stress 
analysis.  Tabulations  for  year  14  can  also  be  found  in  Appendix  C. 

75.  Figures  26  and  27  are  plots  of  the  degrees  of  consolidation 
and  settlement,  respectively,  throughout  the  period  of  deposition  and 
for  9  years  after  deposition  ceases.  Also  shown  in  the  figures  are  the 
results  of  a  conventional  or  small  strain  analysis  of  the  same  disposal 
program  estimated  from  consolidation  charts  (Terzaghi  and  Peck  1967, 

Lambe  and  Whitman  1969).  The  difference  between  the  two  theoretical 
approaches  is  clearly  evident.  The  sudden  drops  in  the  degree  of  con¬ 
solidation  at  years  1  through  7  are  due  to  the  instant  application  of 
additional  dredged  fill  at  those  times.  As  can  be  seen,  90  percent  con¬ 
solidation  is  achieved  at  about  12.8  years  by  finite  strain  theory; 
whereas,  the  deposit  is  only  about  55  percent  consolidated  at  this  time 
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Figure  24.  Void  ratio  distributions  at  the 
end  of  each  year  during  deposition  and 
ultimately 


by  small  strain  theory.  The  predicted  ultimate  settlement  is  essentially 
the  same  in  both  calculations  since  the  original  individual  layer  heights 
were  relatively  small.  It  should  be  noted  that  the  small  strain  analysis 
was  a  hand  calculation  and  more  elaborate  computer  applications  of  the 
theory  may  reduce  somewhat  the  differences  shown,  but  results  from  the 
use  of  the  two  theories  will  never  match  due  to  the  basic  differences 
in  the  theories. 

76.  For  containment  area  design  purposes,  the  results  of  the 
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Figure  25.  Excess  pore  pressure  distribution 
at  years  7,  8,  and  14 


finite  strain  analysis  of  the  dredged  fill  and  compressible  foundation 
can  be  plotted  as  shown  in  Figure  28.  With  this  representation  of  the 
dredged  fill  surface  and  foundation  surface,  the  height  of  containment 
area  dikes  required  during  the  period  of  disposal  can  be  readily 
determined. 


Figure  26.  Degree  of  consolidation  by  finite  strain 
analysis  compared  to  a  small  strain  analysis 


Figure  27.  Settlements  by  finite  strain  analysis  compared 
to  a  small  strain  analysis 
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Figure  28.  Heights  of  the  dredged  fill  and  foundation 
surfaces  during  and  after  disposal  operations 


Consolidation  of  a  Soft  Thick.  Layer 

77.  This  example  will  illustrate  the  program's  capability  to  cal¬ 
culate  primary  consolidation  in  a  soft  thick  layer  which  is  normally 
consolidated  under  a  small  overburden  when  subjected  to  a  series  of 
added  surcharges.  The  layer  is  assumed  to  be  20  ft  thick  and  to  overlie 
a  coarse  sand  so  that  its  lower  boundary  may  be  considered  free  draining. 
The  layer's  void  ratio-effective  stress  and  void  ratio-permeability  re¬ 
lationships  are  those  shown  in  Figures  16  and  20,  and  the  layer's  speci¬ 
fic  gravity  of  solids  was  assumed  to  be  2.80. 

78.  It  is  further  assumed  that  initially  the  top  of  the  layer  is 
about  1  ft  below  the  water  table  and  some  years  ago  was  covered  with  1  ft 
of  sandy  material  so  that  it  is  fully  consolidated  under  about  75  psf  of 
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overburden.  It  is  planned  to  hydraulically  fill  the  area  with  an  addi¬ 
tional  10  ft  of  sand  over  the  next  three  years  to  prepare  it  for  con¬ 
struction  of  light  buildings.  The  sand  will  be  dredged  from  nearby 
sources  and  deposited  according  to  the  schedule  shown  in  Figure  29, 
which  also  depicts  initial  layer  conditions.  It  is  required  to  deter¬ 
mine  consolidation  behavior  of  the  compressible  layer  during  and  sub¬ 
sequent  to  surcharge  additions. 

79.  Based  on  program  calculations,  void  ratio  distributions  can 
be  plotted  for  any  time  during  the  consolidation  process.  Figure  30 
shows  such  distributions  for  the  first  three  years  of  the  example  in 
comparison  to  the  initial  and  final  void  ratios  in  the  layer.  The  dis¬ 
tributions  at  years  1  and  2  are  before  the  surcharges  for  those  years 
are  added.  As  can  be  seen  from  the  figure,  wide  variation  in  void 
ratios  occurs  throughout  the  layer  initially  and  until  it  is  finally  con 
solldated  under  the  total  added  surcharge.  Thus  again,  the  inapplicabil 
ity  of  a  small  strain  analysis  which  assumes  a  constant  distribution  of 
void  ratios  is  manifest. 

80.  The  distribution  of  excess  pore  pressures  at  various  times 
during  consolidation  is  shown  in  Figure  31.  The  principal  information 
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Figure  29.  Schedule  of  surcharges  added  to 
compressible  clay  layer 
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Figure  30.  Void  ratio  distributions 
in  the  compressible  layer 


to  be  gained  from  this  figure  is  the  fallacy  of  the  often-made  assump¬ 
tion  that  the  value  of  the  remaining  excess  pore  pressure  is  its  maximum 
amount  reduced  by  a  percentage  equal  to  the  degree  of  consolidation. 

For  instance,  at  57  percent  consolidation  the  remaining  excess  pore  pres 
sure  is  more  than  89  percent  of  its  maximum  value,  at  76  percent  consoli 
dation  it  is  58  percent,  and  at  91  percent  consolidation,  it  is  about 
25  percent  of  the  original  maximum  value. 

81.  Figures  32  and  33  compare  the  degree  of  consolidation  and 
settlements  respectively  as  predicted  by  the  finite  strain  analysis 
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Figure  31.  Excess  pore  pressure  distribution  in 
the  compressible  layer 


and  a  small  strain  analysis.  Once  again  the  difference  in  the  two  theore¬ 
tical  approaches  is  clearly  evident,  and  as  in  the  dredged  fill  example, 
consolidation  is  predicted  to  occur  at  a  faster  rate  by  the  finite 
strain  analysis.  Even  though  consolidation  occurs  faster,  the  dissipa¬ 
tion  of  excess  pore  pressure  is  predicted  to  occur  slower.  Figure  34 
shows  the  excess  pore  pressure  distribution  by  both  theories  at  year  6 
during  consolidation.  This  figure  shows  that  the  small  strain  theory  is 
uuderconservative  when  used  to  predict  pore  pressures  and  therefore  may 
lead  to  underconservative  safety  factors  when  used  in  stability  analyses. 
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and  small  strain  analyses  for  a  compressible  layer 


Figure  33.  Predicted  settlement  comparison  between  finite  strain 
and  small  strain  analyses  for  a  compressible  layer 


Figure  34.  Excess  pore  pressure  distribution  at  year  6 
in  the  compressible  layer  as  predicted  by  finite  strain 
and  small  strain  analyses 


82.  A  listing  of  problem  input  and  calculations  to  years  3  and  6 
are  included  in  Appendix  C.  The  calculation  constants  x  and  6  were 
1.0  day  and  one-tenth  of  the  layer  height,  respectively.  These  selec¬ 
tions  proved  to  be  sufficient  for  stability  and  provided  for  an  economic 
calculation. 
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PART  VI:  SUMMARY 


83.  This  report  has  developed  the  theory  of  finite  strain  consoli¬ 
dation  in  relatively  simple  and  concise  terms  and  shown  how  the  theory 
can  be  effectively  programmed  for  computer  computation  of  the  consolida¬ 
tion  behavior  of  very  soft  single  or  multiple  layers  of  fine  grained 
materials.  In  the  theory  development,  simplifying  assumptions  have  been 
held  to  a  minimum  which  effectively  makes  the  theory  the  most  general  in 
defining  one-dimensional  consolidation.  The  chief  advantages  of  finite 
strain  theory  over  small  strain  theory  are  its  independence  from  strain 
levels,  its  independence  of  any  set  relationship  between  void  ratio  and 
effective  stress,  and  its  consideration  of  the  variabilities  in  permea¬ 
bility  through  the  consolidating  layer  due  to  changes  in  void  ratio. 

84.  The  computer  program,  CSLFS,  documented  in  this  report 
represents  an  alternative  to  the  conventional  methods  of  calculating 
one- dimensional  consolidation  which  was  previously  unavailable.  The 
program  was  purposely  written  to  require  only  the  most  basic  soil  prop¬ 
erty  data,  i.e.,  point  data  from  laboratory  testing  relating  effective 
stress  and  permeability  to  the  void  ratio.  It  also  provides  for  the  very 
real  case  of  a  semipermeable  boundary.  Although  the  program  was  inten¬ 
tionally  structured  to  facilitate  the  calculation  of  consolidation  in 
multiple  dredged  fill  layers  deposited  on  a  compressible  foundation, 

it  is  equally  suitable  for  making  one-dimensional  consolidation  predic¬ 
tions  in  a  clay  layer  subjected  to  more  traditional  foundation  type 
loads. 

85.  As  shown  by  the  example  problems  worked  in  the  report,  this 
method  of  consolidation  prediction  is  not  merely  a  more  detailed  analy¬ 
sis  which  leads  to  essentially  the  same  results  obtained  through  a 
simpler  small  strain  analysis.  There  is  a  real  and  substantial  differ¬ 
ence  in  the  results  and  indications  are  that  the  finite  strain  method 
is  more  accurate  because  of  consistent  underprediction  of  settlements 
in  designs  using  small  strain  theories.  Therefore,  the  program  should 
prove  to  be  a  valua>le  aid  in  future  designs  requiring  a  prediction  of 
one- dimensional  consolidation  as  a  function  of  time. 


62 


REFERENCES 


Carslaw,  H.  S.  and  Jaeger,  J.  C.  1959.  Conduction  of  Heat  In  Solids, 

2d  ed..  Clarendon  Press,  Oxford. 

Davis,  E.  H.  and  Raymond,  G.  P.  1965.  "A  Non-Linear  Theory  of  Consoli¬ 
dation,"  Geotechnique,  Vol  15,  No.  2,  pp  161-173. 

Gibson,  R.  E.,  England,  G.  L. ,  and  Hussey,  M.  J.  L.  1967.  "The  Theory 
of  One-Dimensional  Consolidation  of  Saturated  Clays.  I.  Finite  Non- 
Linear  Consolidation  of  Thin  Homogeneous  Layers,"  Geotechnique ,  Vol  17, 

No.  3,  pp  261-273. 

Gibson,  R.  E.,  Schiffman,  R.  L. ,  and  Cargill,  K.  W.  1981.  "The  Theory 
of  One-Dimensional  Consolidation  of  Saturated  Clays.  II.  Finite  Non- 
Linear  Consolidation  of  Thick  Homogeneous  Layers,"  Canadian  Geotechnical 
Journal,  Vol  18,  No.  2,  pp  280-293. 

Keller,  H.  B.  1960.  "The  Numerical  Solution  of  Parabolic  Partial  Dif¬ 
ferential  Equations,"  Mathematical  Methods  for  Digital  Computers, 

A.  Ralston  and  H.  S.  Wilf,  Ed.,  Wiley,  New  York,  pp  135-143. 

Lambe,  T.  W.  and  Whitman,  R.  V.  1969.  Soil  Mechanics,  Wiley,  New  York, 
pp  406-421. 

McNabb,  A.  1960.  "A  Mathematical  Treatment  of  One-Dimensional  Soil 
Consolidation,"  Quarterly  of  Applied  Mathematics,  Vol  17,  No.  4, 
pp  337-347.  . . 

Mikasa,  M.  1965.  "The  Consolidation  of  Soft  Clay,  A  New  Consolidation 
Theory  and  Its  Application,"  Reprint  from  Civil  Engineering  in  Japan, 

Japan  Society  of  Civil  Engineers,  Tokyo. 

Ortenblad,  A.  1930.  "Mathematical  Theory  of  the  Process  of  Consolida¬ 
tion  of  Mud  Deposits,"  Journal  of  Mathematics  and  Physics,  Vol  9,  No.  2, 
pp  73-149. 

Palermo,  M.  R. ,  Montgomery,  R.  L.,  and  Poindexter,  M.  E.  1978.  "Guide¬ 
lines  for  Designing,  Operating,  and  Managing  Dredged  Material  Contain¬ 
ment  Areas,"  Technical  Report  DS-78-10,  U.  S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  Miss. 

Palermo,  M.  R. ,  Shields,  F.  D. ,  and  Hayes,  D.  L.  In  press.  "Development 
of  a  Management  Plan  for  Craney  Island  Disposal  Area,"  Technical  Report, 

U.  S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 

Schiffman,  R.  L.  1980.  "Finite  and  Infinitesimal  Strain  Consolidation," 
Journal  of  the  Geotechnical  Engineering  Division,  American  Society  of 
Civil  Engineers,  Vol  106,  No.  GT2,  pp  115-119. 

Schiffman,  R.  L.  and  Gibson,  R.  E.  1964.  "Consolidation  of  Nonhomogen- 
eous  Clay  Layers,"  Journal  of  the  Soil  Mechanics  and  Foundation  Division, 
American  Society  of  Civil  Engineers,  Vol  90,  No.  SM5,  pp  1-30. 

Terzaghi,  K.  1924.  "Die  Theorie  der  Hydrodynamischen  Spanungserscheinun- 
gen  und  ihr  Erdbautechnisches  Answendungsgebeit,"  Proceedings,  First 


63 


APPENDIX  A:  USER’S  MANUAL  FOR  CSLFS 

1.  This  appendix  will  provide  information  useful  to  users  of  the 
computer  program  CSLFS  to  include  a  general  description  of  the  program 
processing  sequence,  definitions  of  principal  variables,  and  format  re¬ 
quirements  for  problem  input.  The  program  was  originally  written  for 
use  on  the  WES  Time-sharing  System  but  could  be  readily  adapted  to  batch 
processing  through  a  card  reader  and  high-speed  line  printer.  Some  out¬ 
put  format  changes  would  be  desirable  if  the  program  were  used  in  batch 
processing  to  improve  efficiency. 

2.  The  program  is  written  in  FORTRAN  IV  computer  language  with 
eight-digit  line  numbers.  However,  characters  9  through  80  are  formatted 
to  conform  to  the  standard  FORTRAN  statement  when  reproduced  in  spaces 

1  through  72  of  a  computer  card.  Program  input  is  through  a  quick  access 
type  file  previously  built  by  the  user.  Output  is  either  to  the  time¬ 
sharing  terminal  or  to  a  file  which  must  be  saved  by  the  user  at  the  end 
of  a  run.  Program  options  will  be  fully  described  in  the  remainder  of 
this  appendix. 

3.  A  listing  of  the  program  is  provided  in  Appendix  B,  and  typi¬ 
cal  solution  output  is  contained  in  Appendix  C. 

Program  Description  and  Components 

4.  CSLFS  is  composed  of  the  main  program  and  ten  subroutines.  It 
is  broken  down  into  subprograms  to  make  modification  and  understanding 
easier.  The  program  is  also  well  documented  throughout  with  comments, 
so  a  detailed  description  will  not  be  given.  However,  an  overview  of 
the  program  structure  is  shown  in  Figure  Al,  and  a  brief  statement  about 
each  part  follows: 

Main  Program.  In  this  part,  input  data  are  read  according  to 

the  option  specified  and  the  various  subroutines 
are  called  to  print  initial  data,  calculate  con¬ 
solidation  and  stresses,  and  print  solution 
output . 

Subroutine  INTRO.  This  subprogram  causes  a  heading  to  be 

printed,  prints  soil  and  calculation  data,  and 


Al 


from  input  relationships  between  void  ratio,  ef¬ 
fective  stress,  and  permeability. 

Subroutine  RESET.  In  this  subroutine  initial  conditions  are 
modified  each  time  a  new  dredged  fill  layer  or 
surcharge  is  added  to  the  consolidating  layers. 
The  subprogram  also  calculates  new  final  settle¬ 
ments  and  resets  the  bottom  boundary  pressure 
gradient . 

Subroutine  FDIFEQ.  This  is  where  consolidation  is  actually 
calculated.  A  finite  difference  equation  is 
solved  for  each  nodal  point  in  the  consolidating 
layers  at  each  time  step  between  specified  out¬ 
put  times.  Void  ratio  functions  and  pore  pres¬ 
sure  gradients  at  layer  boundaries  are  also  re¬ 
calculated  at  each  time  step.  Just  before  each 
output  time,  consistency  and  stability  criteria 
are  checked. 

Subroutine  VRFUNC.  The  functions  a(e)  and  8(e)  required  at 
each  time  step  in  FDIFEQ  are  calculated  in  this 
subprogram. 

Subroutine  STRESS.  Here,  the  current  convective  coordinates, 
soil  stresses,  and  pore  pressures  are  calculated 
for  each  output  time. 

Subroutine  INTGRL.  This  subroutine  evaluates  the  void  ratio 
integral  used  in  determining  convective  coordi¬ 
nates,  settlements,  and  soil  stresses.  The 
procedure  is  by  Simpson's  rule  for  odd  or  even 
numbered  meshes. 

Subroutine  DATOUT.  DATOUT  prints  the  results  of  consolidation 
calculations  and  initial  conditions  in  tabular 
form.  Examples  are  shown  in  Appendix  C. 

Subroutine  DATAIN.  This  subprogram  reads  the  data  from  a  pre¬ 
vious  program  run  so  that  future  consolidation 
can  be  calculated  without  having  to  recalculate 
previous  consolidation. 

Subroutine  SAVDAT.  The  data  from  the  current  program  run  is 
written  to  a  file  in  the  format  required  to  be 
read  by  DATAIN. 


Variables 


5.  The  following  is  a  list  of  the  principal  variables  and  vari¬ 
able  arrays  that  are  used  in  the  computer  program  CSLFS.  The  meaning  of 
each  variable  is  also  given  along  with  other  pertinent  information  about 
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it.  If  the  variable  name  is  followed  by  a  number  in  parentheses,  it  is 
an  array,  and  the  number  denotes  the  current  array  dimensions.  If  these 
dimensions  are  not  sufficient  for  the  problem  to  be  run,  they  must  be 


increased  throughout  the  program. 


A(101) 

Al(ll) 

AF(101) 

AFl(ll) 

AHDF(IO) 

ALPHA(51) 

ALPHA1 (51) 
BETA(51) 

BETA1(51) 

BF(101) 

BFl(ll) 

DA 

DSDE(51) 


the  Lagrangian  coordinate  of  each  space  mesh 
point  in  the  dredged  fill  layers. 

the  Lagrangian  coordinate  of  each  space  mesh 
in  the  compressible  foundation  or  layer. 

the  function  a(e)  corresponding  to  the  cur¬ 
rent  void  ratios  at  each  space  mesh  point  in 
the  dredged  fill  layers. 

the  function  a(e)  corresponding  to  the  cur¬ 
rent  void  ratios  at  each  space  mesh  point  in 
the  compressible  foundation  or  layer. 

the  initial  height  of  added  dredged  fill  layers 
in  Lagrangian  coordinates  or  the  amount  of 
added  surcharge  on  a  compressible  layer. 

the  function  a(e)  corresponding  to  the  void 
ratios  input  when  describing  the  void  ratio- 
effective  stress  and  permeability  relationships 
for  the  dredged  fill. 

the  function  a(e)  as  above  except  for  the  com¬ 
pressible  foundation  or  layer. 

the  function  B(e)  corresponding  to  the  void 
ratios  input  when  describing  the  void  ratio- 
effective  stress  and  permeability  relationships 
for  the  dredged  fill. 

the  function  B(e)  as  above  except  for  the  com¬ 
pressible  foundation  or  layer. 

the  function  B(e)  corresponding  to  the  cur¬ 
rent  void  ratios  at  each  space  mesh  point  in 
the  dredged  fill  layers. 

the  function  B(e)  corresponding  to  the  cur¬ 
rent  void  ratios  at  each  space  mesh  point  in 
the  compressible  foundation  or  layer. 

the  difference  between  the  Lagrangian  coordi¬ 
nates  of  space  mesh  points  in  the  dredged  fill 
layer . 

the  calculated  value  of  corresponding  to 

de  t 

the  void  ratios  input  when  describing  the  void 
ratio-effective  stress  relationship  for  the 
dredged  fill. 
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DSDE1(51) 

DU0 


DUDZ10 

DUDZ11 

DUDZ21 

DZ 

DZ1 

DQ 

E(101) 

E0 

E00 

El(lOl) 

Ell(ll) 

EFFSTR(lOl) 

EFIN(lOl) 

EFINl(ll) 


the  calculated  value  of  as  above  except 

for  the  compressible  foundation  or  layer. 

the  drainage  path  length  in  an  incompressible 
boundary  layer  used  for  computing  the  semi- 
permeable  boundary  condition.  This  value  is 
originally  input  in  Lagrangian  coordinates 
but  is  changed  to  material  coordinates  by  the 
program. 

the  excess  pore  pressure  gradient  in  an  incom¬ 
pressible  foundation  at  its  boundary  with  the 
compressible  layer. 

the  excess  pore  pressure  gradient  in  the  com¬ 
pressible  foundation  or  layer  at  its  boundary 
with  an  incompressible  foundation. 

the  excess  pore  pressure  gradient  in  the  dredged 
fill  layer  at  its  boundary  with  a  compressible 
foundation  or  incompressible  foundation. 

the  difference  between  the  material  or  reduced 
coordinates  of  space  mesh  points  in  the  dredged 
fill. 

the  difference  between  the  material  or  reduced 
coordinates  of  space  mesh  points  in  the  com¬ 
pressible  foundation  or  layer. 

the  initial  additional  surcharge  placed  on  a 
compressible  layer. 

the  current  void  ratios  at  each  space  mesh 
point  in  the  dredged  fill. 

the  void  ratio  in  the  incompressible  foundation 
at  its  boundary  with  the  compressible  layer. 

the  initial  void  ratio  assumed  by  the  dredged 
fill  after  initial  sedimentation  and  before 
consolidation. 

the  initial  void  ratios  at  each  space  mesh 
point  in  the  dredged  fill. 

the  initial  void  ratios  at  each  space  mesh 
point  in  the  compressible  foundation  or  layer. 

the  effective  stress  at  each  space  mesh  point 
in  the  dredged  fill. 

the  final  (100  percent  primary  consolidation) 
void  ratios  at  each  space  mesh  point  in  the 
dredged  fill. 

the  final  (100  percent  primary  consolidation) 
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void  ratios  at  each  space  mesh  point  in  the  com¬ 
pressible  foundation  or  layer. 


EFSTRl(ll) 

the  effective  stress  at  each  space  mesh  point 
in  the  compressible  foundation  or  layer. 

ELL 

the  total  depth  of  the  dredged  fill  in  material 
or  reduced  coordinates. 

ELL1 

the  depth  of  the  compressible  foundation  or 
layer  in  material  or  reduced  coordinates. 

ER(ll) 

the  current  void  ratios  at  each  space  mesh 
point  in  the  compressible  foundation  or  layer. 

ES(51) 

the  void  ratios  input  when  describing  the  void 
ratio-effective  stress  and  permeability  rela¬ 
tionships  in  the  dredged  fill. 

ESI  (.51) 

the  void  ratios  input  when  describing  the  void 
ratio-effective  stress  and  permeability  rela¬ 
tionships  in  the  compressible  foundation  or 
layer. 

F(101) 

the  void  ratios  at  each  space  mesh  point  of  the 
previous  time  step  in  the  dredged  fill. 

Fl(ll) 

the  void  ratios  at  each  space  mesh  point  of  the 
previous  time  step  in  the  compressible  founda¬ 
tion  or  layer. 

FINT(lOl) 

the  void  ratio  integrals  evaluated  from  the 
bottom  to  the  subscripted  space  mesh  point  in 
the  dredged  fill. 

FINTl(ll) 

the  void  ratio  integrals  evaluated  from  the 
bottom  to  the  subscripted  space  mesh  point  in 
the  compressible  foundation  or  layer. 

GC 

the  buoyant  unit  weight  of  the  dredged  fill 
soil  solids. 

GC1 

the  buoyant  unit  weight  of  the  soil  solids  of 
the  compressible  foundation  or  layer. 

GS 

the  unit  weight  of  the  dredged  fill  soil  solids. 

GS1 

the  unit  weight  of  the  soil  solids  of  the 
compressible  foundation  or  layer. 

GSBL 

the  specific  gravity  of  the  soil  solids  of  the 
compressible  foundation  or  layer. 

GSDF 

the  specific  gravity  of  the  dredged  fill  soil 
solids . 

GW 

the  unit  weight  of  water. 

HBL 

the  initial  height  of  the  compressible  founda¬ 
tion  or  layer  in  Lagrangian  coordinates. 
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HDF 


the  initial  height  of  the  first  dredged  fill 
layer  in  Lagrangian  coordinates. 

HDF1  the  initial  height  of  later  dredged  fill  layers 
in  Lagrangian  coordinates. 

IN  an  integer  denoting  the  input  mode  or  device 
for  initial  problem  data  which  has  the  value 
"10"  in  the  present  program. 

INS  an  integer  denoting  the  input  mode  or  device 
for  problem  data  from  a  previous  computer  run 
which  has  the  value  "12"  in  the  present  program. 

I0UT  an  integer  denoting  the  output  mode  or  device 
for  recording  the  results  of  program  computa¬ 
tions  in  a  user's  format  which  has  the  value 
"11"  in  the  present  program. 

IOUTS  an  integer  denoting  the  output  mode  or  device 
for  recording  the  results  of  program  computa¬ 
tions  in  a  format  for  continuing  the  computa¬ 
tions  in  a  later  run  which  has  the  value  "13" 
in  the  present  program. 

LBL  the  number  of  data  points  used  in  describing 

the  void  ratio-effective  stress  and  permeability 
relationships  in  the  compressible  foundation  or 
layer. 

LDF  the  number  of  data  points  as  above  except  for 
the  dredged  fill. 

MTIME  the  number  of  additional  output  times  when  con¬ 
tinuing  a  previous  computer  run. 

NBDIV  the  number  of  parts  the  initial  dredged  fill 

layer  is  divided  into  for  computation  purposes. 

NBDIV1  the  number  of  parts  the  compressible  foundation 
or  layer  is  divided  into  for  computation 
purposes . 

NBL  an  integer  denoting  the  following  options: 

1  =  consolidation  calculated  for  dredged 

fill  layers  and  a  compressible 

foundation. 

2  =  consolidation  calculated  for  dredged 

fill  layers  only. 

3  =  consolidation  calculated  for  a  single 

compressible  layer  only. 

ND  the  total  number  of  space  mesh  points  in  the 
dredged  fill  layers. 


NDATA1 


NDATA2 

NDIV 

NDIV1 

NFLAG 

NM 

NND 

NNN 

NPROB 

NPT 


NST 


an  integer  denoting  the  following  options: 

1  =  this  is  a  new  problem  and  data  will  be 

read  from  file  "10". 

2  =  this  is  a  continuation  of  a  previous 

computer  run  and  data  will  be  read  from 
file  "12". 

an  integer  denoting  the  following  options: 

1  =  do  not  save  data  for  later  computer  run. 

2  =  save  data  on  file  "13”  so  that  calcula¬ 

tions  can  be  continued  in  a  later 
computer  run. 

the  number  of  space  mesh  points  in  the  initial 
dredged  fill  layer. 

the  total  number  of  space  mesh  points  in  the 
compressible  foundation  or  layer. 

an  integer  denoting  the  following: 

0  =  print  current  conditions  heading. 

1  =  print  initial  conditions  heading. 

an  integer  counter  which  is  used  in  tracking 
the  output  times  for  each  computer  run. 

an  integer  used  to  denote  the  total  number  of 
parts  into  which  the  dredged  fill  layers  are 
divided  for  computation  purposes. 

an  integer  counter  which  is  used  in  tracking 
the  total  number  of  time  steps  through  which 
consolidation  has  proceeded. 

an  integer  used  as  a  label  for  the  current 
consolidation  problem. 

an  integer  denoting  the  following  options: 

1  =  make  a  complete  computer  run,  printing 

soil  data,  initial  cnditions,  and  cur¬ 
rent  conditions  for  all  specified 
print  times. 

2  =  make  a  complete  computer  run  but  do  not 

print  soil  data  and  initial  conditions. 

3  =  terminate  computer  run  after  printing 

soil  data  and  initial  conditions. 

an  integer  line  number  used  on  each  line  of 
data  input  and  on  data  lines  output  for  use  in 
a  later  computer  run. 


NTIME  the  number  of  output  times  during  the  initial 
computer  run  of  a  consolidation  problem. 

PK(51)  the  function  j-  —  -  corresponding  to  the 

void  ratios  input  when  describing  the  void 
ratio-permeability  relationship  in  the  dredged 
fill. 

PK0  the  function  -  for  the  incompressible 

foundation  layer. 

k 

PK1(51)  the  function  j —  —  corresponding  to  the  void 

ratios  input  when  describing  the  void  ratio- 
permeability  relationship  in  the  compressible 
foundation  or  layer. 

PRINT(25)  the  real  times  at  which  current  conditions  in 
the  consolidating  layers  will  be  output. 

Q0  the  initial  overburden  on  a  compressible  layer. 

Q1  the  current  total  surcharge  including  overburden 

on  a  compressible  layer. 

RK(51)  the  permeabilities  input  when  describing  the 
void  ratio-permeability  relationship  in  the 
dredged  fill. 

RK1(51)  the  permeabilities  input  as  above  except  for 
the  compressible  foundation  or  layer. 

RS(51)  the  effective  stresses  input  when  describing 
the  void  ratio-effective  stress  relationship 
in  the  dredged  fill. 

RS1(51)  the  effective  stresses  input  as  above  except 
for  the  compressible  foundation  or  layer. 

RWL(IO)  the  new  height  of  free  water  surface  above  the 
bottom  of  the  compressible  foundation  or  layer 
after  a  new  dredged  fill  layer  or  surcharge 
has  been  added. 

SETT  the  current  settlement  in  the  dredged  fill. 

SETT1  the  current  settlement  in  the  compressible 

foundation  or  layer. 

SFIN  the  final  settlement  in  the  dredged  fill  layer 
presently  existing. 

SFIN1  the  final  settlement  in  the  compressible  founda¬ 
tion  or  layer  under  present  loading  conditions. 

TAU  the  value  of  the  time  step  in  the  finite  dif¬ 
ference  calculations. 


TIME  the  real  time  value  after  each  time  step. 

TPRINT  the  real  time  value  of  the  next  output  point. 

TOSTRl(ll)  the  current  total  stress  at  each  space  mesh 

point  in  the  compressible  foundation  or  layer. 

TOTSTR(lOl)  the  current  total  stress  at  each  space  mesh 
point  in  the  dredged  fill. 

U(101)  the  current  excess  pore  pressure  at  each  space 
mesh  point  in  the  dredged  fill. 

U0(1O1)  the  current  static  pore  pressure  at  each  space 
mesh  point  in  the  dredged  fill. 

1101(11)  the  current  static  pore  pressure  at  each  space 
mesh  point  in  the  compressible  foundation  or 
layer. 

Ul(ll)  the  current  excess  pore  pressure  at  each  space 
mesh  point  in  the  compressible  foundation  or 
layer. 

UCON  the  current  degree  of  consolidation  in  the 
dredged  fill. 

UC0N1  the  current  degree  of  consolidation  in  the  com¬ 
pressible  foundation  or  layer. 

UW(101)  the  current  total  pore  pressure  at  each  space 
mesh  point  in  the  dredged  fill. 

UWl(ll)  the  current  total  pore  pressure  at  each  space 
mesh  point  in  the  compressible  foundation  or 
layer. 

VRI1  the  initial  total  void  ratio  integral  for  the 
compressible  foundation  or  layer. 

WL  the  initial  height  of  free  water  surface  above 
the  bottom  of  the  first  dredged  fill  layer. 

WL1  the  initial  height  of  free  water  surface  above 
the  bottom  of  the  compressible  foundation  or 
layer. 

XI(101)  the  current  convective  coordinate  of  each  space 
mesh  point  in  the  dredged  fill. 

XIl(ll)  the  current  convective  coordinate  of  each  space 
mesh  point  in  the  compressible  foundation  or 
layer . 

Z(101)  the  material  or  reduced  coordinate  of  each 
space  mesh  point  in  the  dredged  fill. 

Zl(ll)  the  material  or  reduced  coordinate  of  each 

space  mesh  point  in  the  compressible  foundation 
or  layer. 
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ZK0  the  permeability  in  the  incompressible  founda¬ 
tion  at  its  boundary  with  the  compressible 
layer. 


Problem  Data  Input 


6.  The  method  of  inputting  problem  data  in  CSLFS  is  by  a  free 
field  data  file  containing  line  numbers.  The  line  number  must  be  eight 
characters  or  less  for  ease  in  file  editing  and  must  be  followed  by  a 
blank  space.  The  remaining  items  of  data  on  each  line  must  be  separated 
by  a  comma  or  blank  space.  Real  data  may  be  either  written  in  exponen¬ 
tial  or  fixed  decimal  formats,  but  integer  data  must  be  written  without 
a  decimal. 

7.  For  an  initial  problem  run  (i.e.,  NDATA1  =  1),  the  data  file 
should  be  sequenced  in  the  following  manner: 

a.  NST,  NPROB,  NDATA1 ,  NDATA2 

b.  NST,  NPT,  NBL 

c.  NST,  GSBL,  HBL,  WL1,  LBL,  Q0,  DQ 

d.  NST,  ES1(I) ,  RS1(I) ,  RK1(I) 

e.  NST,  GSDF,  HDF,  WL,  LDF,  E00,  GW 

£.  NST,  ES(I) ,  RS( I) ,  RK(I) 

£.  NST,  E0,  ZK0,  DU0 

h.  NST,  NBDIV,  NBDIV1,  TAU,  NTIME 

i_.  NST,  PRINT(I) ,  AHDF(I),  RWL(I) 

It  should  be  pointed  out  here  that  NST  may  be  any  positive  integer  but 
must  increase  throughout  the  file  so  that  it  will  be  read  in  the  correct 
sequence  in  the  time-sharing  system. 

8.  The  following  exceptions  and  explanations  should  also  be  noted 
for  particular  line  types: 

Line  type  c:  Q0  and  DQ  have  nonzero  values  only  if  NBL  =  3. 

If  NBL  -  2,  all  data  values  are  set  to  zero 
except  NST. 

Line  type  d:  There  are  LBL  of  these  lines  unless  NBL  =  2,  and 
then  there  will  be  one  line  with  all  values  set 
to  zero  except  NST. 
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Line  type  e:  If  NBL  =  3,  all  values  on  this  line  are  set  to 
zero  except  NST  and  GW. 

Line  type  f:  Ther  are  LDF  of  these  linei  unless  NBL  =  3,  and 
then  there  will  be  one  line  with  all  values  set 
to  zero  except  NST. 

Line  type  i:  There  are  NTIME  of  these  lines. 

9.  For  the  continuation  of  a  previous  problem  run  (i.e., 

NDATA1  =  2),  the  input  data  file  should  be  input  in  the  following 
sequence: 

Line  type  aa.  NST,  NPROB,  NDATA1,  NDATA2 

Lae  type  bb.  NST,  MTIME 

Line  type  cc.  NCT,  AHDF (NTIME),  RWL (NTIME) 

Line  type  dd.  NST,  PRINT(I),  AHDF ( I ) ,  RWL(I) 

10.  The  following  explanations  should  be  noted  for  particular 
line  types: 

Line  type  cc:  AHDF  and  RWL  are  the  values  from  the  last  line 
of  the  previous  computer  run. 

Line  type  dd:  There  are  MTIME  of  the  lines. 

11.  All  input  data  having  particular  units  must  be  consistent  with 
all  other  data.  For  example,  if  layer  thickness  is  in  feet  and  time  is 
in  days,  then  permeability  must  be  in  feet  per  day.  If  stresses  are  in 
pounds  per  square  foot,  then  unit  weights  must  be  in  pounds  per  cubic 
foot.  Any  system  of  units  is  permissible  so  long  as  consistency  is 
maintained. 


APPENDIX  B:  CSLFS  PROGRAM  LISTING 


1.  The  following  is  a  complete  listing  of  CSLFS  as  written  for 
the  WES  time-sharing  system. 
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noni  CiufifH  CSLFS  CONSOLIDATION  OF  -.OFT  LAYERS  BY  FINITE  STRAIN 


n  ooi  om  oc 

0 in'll  On? or 

0  n  o  i  o  ii  ■  or 

♦  4 

•10  01  ft  ft  4  I'lC 

♦ 

♦ 

0  o  o  i  o  ns  or 

♦ 

CSLFS 

♦ 

on  oi  one  or 

♦ 

♦ 

0  0  f  1  0  0  r  or 

♦ 

ONE-DIMENSIONAL  FINITE  STRAIN 

CONSOLIDATION 

♦ 

n  o  o  1 0  ns  nr 

♦ 

♦ 

0  o  o i on roc 

♦ 

OF 

♦ 

nooi m  nor 

♦ 

♦ 

0  0  0 1  0 1  IOC 

♦ 

HOMOGENEOU:  SOFT  CLAY 

LAYERS 

♦ 

0  0  0  1  0 1  ?  oc 

♦ 

♦ 

0  o  o  i  o  i  or 

♦  4 

0  0  o  1014  or 

ii  no  i  n  1 5  or 

i'i  i'i  0  101k  I'm" 

i,  ri  i'i  1  ri  1  7  fiC 
rn'i ill  01  HOC 
nor  i  of?  nr 
0  n  fi  1 1  !=■  0  nr 
i'i  no  i  >'•?  i  nr 
0  om  oss  or 
on ii i  ri i  -■  o r 

0  0  01  0.?  4  01- 

no  0 1  0?  5  or 
0  o  0 1  it?  6  or 
nooi  n??n>; 

I'lOn  1  I'i,?:-:  iV 

0  I'I  01  I'i?"?  I'lC 
0001  0?  I'  t'L 

o  ri  o  i  o  ?.  i  or 

o  0  I'I  i  o  3?  0 

no  oi  o?  '30 

0  0  0 1  0  •;  4  0 
00  01  OS'50 
00  01  0 Si  ft 
011  01  0370 
on 01 03? o 
on oi o? ao 
0001 0400 
000 l 04  I  0 

Ol' 01  I'i 4?  0 

i'i  0  01  04?  0 
0  0  01  044(1 
I'i  001  04?  0 
0  0  0 1  044  nr 


♦  CSLFS  COMPUTES  THE  VOID  PATIOS.  TOTAL  ANP  EFFECTIVE  ♦ 

♦  STRESSES*  F'DPE  WATER  PRESSURES*  .SETTLEMENTS*  AND  ♦ 

♦  DEGREES  OF  CONSOLIDATION  FOP  HOMOGENEOUS  IOFT  CLAY  ♦ 

♦  L  Ft  YEP S  OF  DREDGED  FILL  DEPOSITED  ON  A  COMPRESS  I FLE  ♦ 

♦  OP  INCOMPRESSIBLE  LAYER  BY  FINITE  STRAIN  CONSOLIDATION 

♦  THEORY.  LCMFP  BOUNDARY  OF  THE  BOTTOM  COMPRESSIBLE  ♦ 

♦  LAYER  MAY  BE  COMPLETELY  FREE  DRAINING.  IMPERMEABLE.  ♦ 

«  OP  NEITHER.  THE  VOID  RATIO-EFFECTIVE  STRESS  AND  ♦ 

♦  VOID  PATIO-PERMEABILITY  RELATIONSHIPS  ARE  INPUT  AS  ♦ 

♦  POINT  VALUES  AND  THUS  MAY  ASSUME  ANY  FORM.  ♦ 


COMMON  DA. DUO. DUDZ1 0.DUDZ11 . DUDSRl . DZ >  IC 1 . DC  >  E  0. E  0  0. ELL  .ELL! . 
GC . GC 1 . GS . GS 1 . GSBl . G IDF? GM . HBL .HDF.HDF1. IN. IN? • I DUT . 

I  Ol  ITS .  LBL .  LDF .  MT I  ME .  NE  Ii  I V .  NBD I V 1  .  NBL  .  ND .  HP  I V .  ND I  V 1 . 
NFLAG.NM.NPPOE.NPT.NNIi.NNN.NTIME.PI':  0.00.  Cl «  SETT .  SETT  1  • 
5F  IN.SFIN1.TAU.TIME.  TF'P  I  NT .  UCON  •  UC  ON  1 .  VP  1 1 .  ML  .  ML  1 «  Zt  0 , 

ao  on »  ai  an  .  af  a  on «  afi  <  in  .alpha  <5  n .  alpha  1  rsi  > . 

BETA  -5 1>  .  BETA  1  ''5 1  >  .  BF  ■  i  (.  i  >  .  BF 1  •  1 1  > .  DS  PE  <5 1  >  .  DSDE 1  •  5 1  >  . 
€  •:  1  0 1  •  .El  a  0  n  •  E 1 1  <11  •  •  EF  I N  a  0 1>  .  EF I N 1  an  »  EP  •' 1  1'  . 

ES  « 1  > .  ES 1  >'5 1  > .  EFFSTP  <  1  0 1 ')  » EFSTP 1  •  1 1  •  -  F  <1  0 1 '»  .  F 1  •  1 1  ■  • 

FiNTaon .pint i  an .Fiasn  .pta *5n  >pi  <sn .wa  <'5i> . 
ps -5i '  .psi  '5i -  .totstp  (ion  .TQSTPi  an  *u> ion .ui  an. 
u o ■  ion  *uoi  at  •  »um •  1  on  *uwi  ■  in  *xi  <ion  .xn  a  n . 
zaoi  •  .zi  a  i  ■ 

D I  MEN  S I  ON  AHDF  <  1  0  >  .  PRINT  •:  35Y  *  PJ.iL  ■  1  O' 
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0001 04700 

0001 048  OC  ...SET  IHPIJT  AND  OUTPUT  NODES 

00010490  IN  =  10 

00010500  I OUT  =  11 

00010510  INS  =  12 

00010520  I OUTS  =  13 

00  01 053 Of  ...READ  PROBLEM  INPUT  FROM  FREE  FIELD  DATA  FILE 

oOOl 054  00  . CONTAINING  LINE  NUMBERS 

00010550  100  FORMAT  <\>  > 

000105600  . PROBLEM  NUMBER.  DATA  OPTIONS.  INTRO  OPTION.  FDT  OPTION 

00010570  READ  (IN.  100.)  NST.  NPROB.  NDATA1 .  NDATA2 

0001 0590  IF  i'NDATAI  .EO.  2>  GOTO  4 

00010590  READ  > IN. 100*  NST.NPT.NBL 

OOrtlOGfiOC  . SOIL  DATA  FOP  FOUNDATION  LAYER  OP  SOFT  LAYER 

00 0 1 06 1 0  READ  f IN. 1 n 0>  NST. GSBL » HBL . ML 1 >  LBL . 0 0 . DO 

00010620  DO  1  1=1. LBL 

0001 0630  READ -  IN. 100)  NST. ESI  CD  > RSI  *'P  .PK  1  '  I) 

00010640  1  CONTINUE 

0001 065 OC  . SOIL  DATA  FOR  DREDGED  FILL 

00010660  READ  ‘'IN.  1 00)  NST. GSDF. HBF. ML. LDF. E00. GW 

00010670  f  DO  2  1=1. LDF 

00010660  READ 0 IN. 1 00)  NST. ES <I> .RS <1 > . RK 01) 

0001 0690  2  CONTINUE 

OOOl 070 Of  . CONSOLIDATION  CALCULATION  DATA 

0 1)010710  READ  •  I N .  1 0  0)  NST .  E  0 ,  Zfc  0 .  DU  0 

U0010720  READ <  IN. 100)  NST. NBD1 V. NBDI VI . TAU. NT IME 

00010730  DO  3  T= 1 . NTIME 

OOOl 074  0  READ i IN. 1 00)  NST.  PRINT  Cl )  « ANDF <  I ) . RWL  <I> 

00010750  3  CONTINUE 

0001 076 OC 

OOOl 077 OC  ...SET  INITIAL  VARIABLES 


n  n o 1076  o 

ELL1  =  0.0  ? 

DZ1  =  0.0 

nooto79o 

TIME  =  0.0 

OOOl 0300 

ucon  =  n. o  ; 

UCON  1  =  0.0 

OOOl 031 0 

SETT  =  0. 0  ! 

SETT1  =  0.0 

0001 0320 

spin  =  o.o  ; 

SFIN1  =  0.0  !  VRI 1  = 

OOOl 0930 

NNN  =  1  .  NM 

=  1 

0 0  0 10340 

DA  =  0.0  ?  D2 

=1.05  HDF1  =  0. 0 

OOOl 035 n 

DUDZ1 1  =0,0 

5  DUD221  =0.0 

00010860 

01  =  00  +  DO 

0001037  OC 
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0001  088 nr  .  .  .PRINT  INPUT  PRTR  RNl>  MfilE  INITIAL  CALCULATION" 

0001  0*9 1'  .CALL  I  NT  PCI 

00010900  IF  .HPT  .EO.  3>  STOP 

00010*10  GOTO  6 

0001 09c 00 

oooio9:-or  . . , nfw  consolidation  time:  and  dpt r 

00010940  4  REftflCIN. 100>  NST.MTIME 

0001 0900  CRLL  DATA  IN 

0001  0960  PEAPtIN,  100>  NST.RHDFfMM-l>,PHLoHM-t> 

00010-70  DO  T  I =HM. HTIMF 

00010990  READ f IN,  1  On  .  NST, PRINT  ''!•  .  RHPF  .  I  :.  , Pf,<L  ‘  I > 

00010991'  5  COMTIHUF 

1)0  01 1  00  Of 

0001  1  01 OC  _ PERFORM  CALCULATIONS  TO  EACH  PPINT  TIME  PHD  OUTPUT  RESULTS 

OOOllOE'O  6  DD  9  K=HM,HTIME 

0 0  o  11050  Tpc  I  NT  =  PRINT  oo 

0001  1040  IF  •:>:  .EO.  1>  GOTO  7 

00 01 1 050  HpPl  -  AH  OF ■k-l. 

00  o  1  1  06  0  I.1L  1  =  F  I.JL  <K  - 1  > 

0001 1070  CRLL.  RESET 

000)1 080  7  CRLL  FDIFEO 

00011090  CRLL  “TRESS 

0001 11 00  CRLL  DR TOUT 

0 OO11110  3  CONTI NOE 

00011  1 5  OC 

0001 1 130  IF  CNPA TR5  .£0.  2>  CRLL  SAVPRT 

0  0  0 1  1  14  Of. 

00011  ISO  STOP 

Coon  160  END 

0001 1 1 7  OC 
OOOJ 1 130C 
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n  n  n?  n  o  r.  o 
no  os  ofii  iif 
mi  n? nn? nr 
nun?  on  ■?  nc 
o  iiii?  iii  14  or 

i  lill'V-  Of:  5  ff 
n  o  ii?  o  ilk  lie 
mm?  iin  ?n 
i'i  il  i'i?  il  f.  9  o 

o  on?  on-?  i' 

n  fin?  m  iin 
non? ni  l  ii 
000201? i ' 

n nil?  nt  ;n 
nnn? n  1 40 
non?  ni  5  ii 
Il  I'I  n  ?  o  1  ?.  o 

ii  0  n?  o ;  7  ii 

i  i  il  0  ?  ft  1 9  n 

0  o n?  o  1 9  o 
n  0  il?  0?  0  oc 
n  n  o?  o?  i  of 

Oil  i'i?  il??  i'i 
il  0  il?  1.1?.?  il 
0  ii  ii?  n?4  ii 
il  0  f1?:  025  f. 
Ilfill?fl?kli 

ii  o  n?  o?  7  o 
noil?  0?  i?  n 
no  o?  Oc"?  nc 

n  0  0?  03  0  il 

r.  oo?  ns  io 
non?o.3?0 
00020330 
il  nil?  09411 

non? 0350 

il  0  f.  ?  il  ?  ?  il 

non? 03 70 

n  0  i'i?  0  ?  ft  o 
0  0  0?  039  il 
il  0  03  04  il  il 
0003041  111? 
00030430 
n  n  n?  04  3  0 
00030440 
0003045 f 

o  no?  04  ?.  o 
00030470 
00030490 
0  Oil 5  04 90 
noo?n5.nn 
non? 051 0 


SUBROUTINE  INTRO 


♦  INTRO  PRINT?  INCUT  E«r9T«=*  AND  RESULTS  OF  INITIAL  ♦ 

♦  CALCULATIONS  IN  TABULAR  FORM.  ♦ 


COMMON  DA • DU 0 *  DUDS 1 0 , DUDS 11 * DUDS? 1 ,  DZ ,  DZ 1  *  DU .  E 0 * E  0 n , ELL  *  ELL  1 . 
ft  GC * GC 1  * GS* GS 1  * GSBL * GS DF *  GW • HBL • HDF * H OF 1  *  IN*  INS  *  IOUT* 

ft  TOUTS  *  LBL .  LDF*  MT IME,  NED  I V*  HDD  IV 1  *  NBL  *  ND.  NDTV.  NDIVl . 

ft  NFLAG*  NM*  NPPOB  *  NPT  <  NMD  >  NNN  *  1ST  I  ME*  PKO*  C>0*  01  *  SETT*  SETT  1 , 

ft  SFIH*  SFINl.TAU,  TIME*  TPPINT*  UCGN*  UC0H1  >  VP  I  1 ,  WL.*  WL1  *  ZKO* 

?<  A  •  1  0 1 »  •  A 1  •  1 1  >  *  AF  *'  1  0 1  >  •  AF 1  <  1 1  '  *  ALPHA  <5 1  >  *  ALPHA  1  ■  5 1  >  * 

BETA <5 11  *  BETA1  <51  ■  ,BF<101 *  ,  BC1  ■'  1 1  ■  •  DSDE «  5 1  >  *  DS  DE1  <511  > 
ft  E<  1011  *  El  <101  >  .Ell  >1 1  >  *  ERIN  <1 01  • « EFIN1  -  111  ,EP<1 1  >  * 

ft  ES <51>  * ES 1 *511 .EFFSTR <1011 .EFSTP1 <111  *F' 1011  *  FI <11 1 . 

ft  FINT<101) • F I NT 1 <111  *  PL  <51 > *PK1 <51  > *PK<511  *  PI  1 <511  . 

PS  <511  .PS)  <51  >  *  TOTS  TP  <  101*  * T0STF1  <11  *  *IJ<1  Oil  ,  U1  <1  1  >  . 

ft  u  0  <  1011*  u  0  lvin.  uw  <  1 0 1 1  *  uw  1*111*  x  1  <  1 0 1  > . : :  1 1  •  1 1 1 . 

ft  Z<1  01  >  >Z1  *  1  1 1 

...PRINT  PROBLEM  NUMBER  AND  HEADING 

WRITE • IOUT* 1 001 

WRITE ‘IOUT* ion 

WRITE < I OUT. 1 0? 1 

WP I TE  <  I  OUT  *  1  03'*  NP  ROB 

CALL.  SETUP 

IF  <UPT  .  EC .  31  PETIJFH 
IF  <NBL  . EO.  31  GOTO  3 

...PRINT  SOIL  DATA  FDR  COMPRESSIBLE  FOUNDATION 
WRITE  < IOUT  ? 1041 
WRITE < IOUT. 105) 

WRITE < I OUT, 1 0G> 

WR I TE • IOUT , 1 071  HBL , GSBL *  WL 1 , 0 0 
WRITE < IOUT* 1031 
WRITE  < IOUT , 1091 
DO  1  1=1* LBL 

WRITE  *  IOUT,  1 1  01  I.  ESI  <11  *  RSI  <H  ,  RK  1  <11  »PK1  <11  ,  B.ETA1  <11  * 
ft  DSDE 1 <11, ALPHA 1 < I 1 

1  CONTINUE 

IF  -NBL  .EC.  31  GOTO  4 
...PRINT  SOIL  DATA  FOP  DREDGED  FILL 
3  WRITE* IOUT* 11 11 
WRITE  < IOUT  *  1121 
WRITE <IOUT, 1131 

WR  I  TE  <  I  OUT  *114:*  HPF  *  G  SDF  *  WL  *  E  0  0 ,  GW 
WRITE < I OUT* 1031 
WRITE  < I OUT  *  1 091 
DO  3  1=1* LDF 

WR  I  TE  <  I  OUT  *  1 1  O'*  I  *  ES  <  I  >  ,  PS  <  I )  *  PI*  <  1 1  *  PK  •"  1 1  *  BETA  ■:  I  •  , 
ft  DSDE <11 , ALPHA* I* 

•3  CONTINUE 
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00 02 052 DC  ...PRINT  CALCULATION  DATA 

00020530  4  MPITFt IDUT, 115) 

Cl 00$ 05-4  0  WRITE  <  IDUT ,  116) 

0 O 02 055 0  MR ITE  , IDUT > 1  ID 

il  0  it?  056  0  UP  I TE  <1 OUT » 1 1 8)  TAU , E  0 , 2K  0 » DU  0 

0 0 02 057 OC  ...PRINT  TABLES  QF  INITIAL  CONDITIONS 

00020580  NFL  AG  =  1 

ooo?  056.0  CALL  DATOUT 

00  0?  06  00  NFLAG  =  0 

00 OP 061 OC 

O0O?06?OP  .  .  .  FORMATS 

0 0 0? 06  3 0  100  FORMAT '1  HI  xxxxx9X,60  ilH*.  ) 

0 00? 064  0  101  FORMAT  C9X *  49HCDNSQL I PAT ION  OF  SOFT  LAYERS  BY  FINITE  STRAIN  —  . 

000? 0650  g  12HPREDGF.D  FILL) 

0 0 02 066 0  1  0?  FORMAT  ,9X»  60  < 1  H*>  > 

0  00? 0670  103  FORMAT  >  .- 9X,  14HPP0BLEM  NUMBER.  14) 

0  0  0? 0630  104  FORMAT  , vxxxxis* ih*> • 37HS0IL  DATA  FDR  COMPRESSIBLE  FOUNDATION. 

0  0  02  063  0  X  1 7  >'  1 H  ♦ '  > 

0002  070  0  105  FORMAT <x'6X,5HLAYEP, 6X« 16HSPECIFIC  GRAVITY, 4X> 1 1HWATEP  LEVEL, 

0  0 0?  0?  1  0  S-  ?X,  7HINITI  Al.) 

000?07?0  106  FORMAT  .'4X,  9HTHICKNESS  •  AX,  9H0F  SDL  IDS,  7X.  1  1HFRDM  BOTTOM,  8X, 

n  00?  07?  0  g.  ?H  SURCHARGE) 

000? 074  0  1 07  FORMAT <  '4X« R8. 3, 7X.F8. 3, 2  <10X»F8.3) ) 

0  0  0?  075  0  1  0?  FORMAT  OxAX « 4H  VO  I D ,  2X ,  9HEFFECT I VE  •  3X  <  5HPEFM- ,  5X ,  5HK/ 1  +E  > 

0  p  02  076  o  l  n?  FORMAT  ■:  4X,  SH I  RAT  1 0 , 4X ,  6HSTRESS ,  3X ,  8HEAE I  LI  TY ,  4X ,  £HPK ,  7X «  4HBETA  • 
00020770  *,  6X ,  4HDSDE » 5X ,  5HALPHA) 

00020780  110  FORMAT  <£X,  IS-.  IX,  F6. 3, 6E1 0.  3) 

00020730  111  FORMAT <1HV>  ,£6HS0IL  DATA  FDP  DREDGED  FILL,  23 <1H*> , 

00 0? 03 0 0  112  FORMAT <x x5X, 5HLAYER, 5X, 16HSREC IF IC  GRAVITY, 3X,  1 1HWATEP  LEVEL- 

0 0 02 OS 1 0  %  5X , 7HINIT I AL , 4X . 11 HUN I T  WEIGHT) 

00020*20  113  FORMAT  <■  3X •  9HTH I CKNES S ,  7X .  3H0F  SOL  IDS ,  6X,  1 1HFR0M  BOTTOM, 

00020A30.  *.  3X,  1  OH  VO  ID  RAT  ID,  5X,  8HDF  WATFP) 

00020840  114  FORMAT (.  2X,F3. 3, 8X-F8. 3, 9X, F8. 3, 5X.F8.3, 7X, F6. 2) 

00020850  115  FORMAT  <Vxxxx?8C1H4>, 16MCALCULATI0N  DATA,  £8  «'1H+)  ) 

00020860  116  FORMAT • xxSX, 3HTAU, 1  OX. 11HL0WEP  LAYER, 7X, 1 1HL0WEP  LAYER, 7X, 

no 020870  *,  13HDRA INAGE  PATH) 

00 0? 0880  117  FORMAT 02 IX, 10HVOID  RATIO, 8X, 12HRERMEABILI TV ,9X . 6HLENGTH • 

00020830  118  FORMAT  ox4X, El  1 . 5, 8X. F8. 3, 9X* E 1 1 . 5, 7X, 3HZ  »,F8.3» 

00020900C 
0  0  02  03 1  OC 

00020320 

0 00203  3  0 
0 0 02  094  OC 
0  0  02  095  OC 


RETURN 

END 


SUBROUTINE  SETUP 


0  0  0  5  ft  c>  0  0 
000 30 i'll  OC 
n ou3  or 3  nr 

on OSOOOijC  ♦  SETUP  MALES  INITIAL  CALCULATIONS  AND  MANIPULATIONS  ♦ 

0003004OC  ♦  OF  INPUT  DATA  FQP  LATER  USE.  ♦ 

on 03 005 oc  ♦♦♦♦♦♦• 

0  0  0  3  0  Ob  OC 

no  030*170  COMMON 

000  3009  0  ‘i- 

Ii00?0090  x 

0  0  0  3010  0  & 

0003011  0  3. 

0  0  0  3  0 1  2  0  3. 

0 1 1 0  3  0 1  90  i. 

00030140  5 

000.90150  :i 

00  03  0 160  4= 

0  0  03  01 70  f. 

00030190  &, 

00030190  & 

0  009  0?  0  Of 

0003021 OC  ...SET  CONSTANTS 

O0030220  MB  I V.  =  NBDIV  +  1 

00030230  MB  =  Nil  IV 

00030240  62  =  GSDF  ♦  GW 

00030250  GC  =  5S  -  GW 

00030260  GS1  =  6SBL  ♦  GW 

0003027*0  GC  1  =  GS  1  —  GW 

00030230  NDJV1  =  NBPIV1  +  1 

00 03 0290  PKO  =  ZKO  •  <1.0+E0* 

00030295  DUO  =  DUO  • 1. 0+E0> 

00030300  IF  <NBL  .EC.  2>  SOTO  10 

0003031 OC 

00 09 092 OC  . . . CALCULATE  ELL  POP  COMPRESSIBLE  FOUNDATION  LAVER 

00030330  PZZ  =  0.  0 

00090340  NBD  =  10  ♦  NBDIV1 

000:30350  DABL  «  HEL  /  FLOAT  <NBP> 

00 03 0360  EFS  =  QO 

00030370  DO  4  1=1, NED 

00030330  DO  1  N=2>  LEL 

00030390  SI  =  EFS  -  PS! <N> 

00 03 04 00  IF  -SI  .LE.  0. 0>  GOTO  2 

00030410  1  CONTINUE 

00030420  V  =  ESI  'LEL*1  5  GOTO  3 

00030430  2  NN  =  N-l 

0 0 030440  V  =  ES  1  < N >  +  <* S I ♦  < ES 1  •  NN>  -ESI  <N>>/* <RS  1  <NN>  -PS  1  < N>  >  > 

00030450  3  TDZ  =  DAFL  <1.0+V> 

00030460  EFS  =  EFS  ♦  GC1*TDZ 

00030470  DZZ  *  DZZ  +  TDZ 

00030430  4  CONTINUE 

00030490  ELL1  =  DZZ 

00030500  DZ1  =  ELL  1  ✓  FLOAT  <NBDIVO 

0003051 OC 


DA . DU 0 *  PUDZ 1 0 , DUPZ 1 1 » DUDZ2 1 , DZ > DZ 1 > DO , E  0 ,  E  0 0 ,  ELL , ELL  1 
GC  *  GC  1  1  bS ,  6S  1 , GsBl. ,  bS DF , GW ,  HBL ,  HI'F ,  HDF 1 ,  IN,  I  NS  ,  I DUT , 

I OUTS >  LEL »  LPF , MT I ME , NED I V , NED I V 1 , NBL , ND , ND I V , ND I V 1 , 
NFLAG,  NM,  NPRC1B,  NPT ,  NND,  NNN,  NT  IME ,  PKO,  00,  01 »  SETT ,  SETT1 
SF I N , SF I N 1 , TAU , T I ME - TPP I NT , UCON ,  UC ON 1 • VP 1 1 , ML » WL 1 , ZK 0 
A  <100  »A1  <  1 1  •  » AF<1  01  >  ,  AF1  <  1 1  >  •  ALPHA  '5 1  >  ,  ALPHA  1  '50  , 
BETA  <5 1  > ,  EE T A 1  <5 1  >  ,  BF  • 1  0 1 »  ,  BF 1  <  1 1  >  ,  DSDE  <50,  DSDE 1  <5 1 1 
E<101>  ,E1  <100  »E1 1  <11<  >EFIN<100  .EFIN1  <1  O  *EP  <1 1>  , 

ES  <51> , ESI <51 > , EFFSTP  <101) »EFSTP1 ' 1 1 > ,F' 1 0 1 > ,F1 <10 , 
FINT  <1  OO  •  FINT1  <1  O  » PK  <50  ,  PF1  <51  >  ,  RK  <51  >  ,  PK1  <50  , 

PS  <50  »RS1  <5t>  ,  TOTSTP  <100  ,T0STR1  <10  , U  <  1 OO  ,  U1  <  1 1>  , 

uoaoo  «uoi  <1  o  »uw<ion  ,u«i  <11  >  »xi  <1  oi>  ,xi  1  <1  o  , 

Z  <  101  ,Z1  <10 


0  II 0  3  052  OC 

.  .  .  CALCULATE  INITIRL  COOPDINRTES  AND  VO 1 0  PRTID7 

oon?o5?or 

...POP  COMPRESSIBLE  FOUNDATION  LAYER 

no 03 H540 

21 <i>  =o. o  ;  Ri‘iiso.0  ;  :<i  1  a  >=o.  o 

0  0  0 3  055 C 

FFS  =  SCI  ♦  ELL1  +  00 

00030560 

DO  8  1  =  1 «  NBIV1 

00030570 

DO  5  N=c'«  LBl 

00 03 053 0 

SI  =  EPS  -  PS1 ONl 

00  0-:  0590 

IF  <S  1  .LE.  0.01  SOTO  6 

000  -  0600 

s 

com  I  HUE 

00ft 30*  1 

Ell  <11  =  ES1CLBL.1  !  SOTO  7 

000306c 0 

h 

MH  =  N-l 

00030630 

Ell  '  ll  =  ES1CN1  +  <S1*<ES1  'NN1-ES1  <N>  >  ✓  <RS1  <NN>-RS1  <1 

00030640 

? 

F1<I>  =  Ell <11 

00030650 

EPCIl  =  E11CI1 

0  ft 03  ftb*  0 

EPS  =  EFS  -  SC14D21 

00  03  06  70 

s 

CONTINUE 

00 03 0650 

CALL  I NTGPL  v£R  >  DZ 1 >  NB I V 1 >  F I NT 1 > 

0  0  03  Oe.9  0 

DO  9  I=?>  NDIV1 

00030700 

21  <11  =  Z1 CI-l)  +  D21 

00030710 

fll<I)  =  21<I>  +  FINTl  '11 

00  0":  0730 

xii  <■  i  >  =  Riai 

00030730 

00 030740C 

9 

CONTINUE 

00 03 075 OC 

...CRLCULRTE  ELL  FOR  FIRST  PRFBGED  FILL  LAYER 

00030760 

1  0 

ELL  =  HPF  ✓  a.O+EOO.) 

00030770 

0 00307300 

IF  <NBL  .EO.  3>  GOTO  15 

00  03  079  Of. 

...CRLCULRTE  INITIRL  C00RD1NRTES  RND  SET  VOID  RRTIOS 

00 03 03 00 

DZ  =  ELL  x  FLORT  UIPPIV. 

00 03 03! 0 

2<l)»o.o  «  R <  1)  =0.  0  ;  xr <11*0.0 

00030330 

€1 < 1 1 =E  00  !  F  < 1 1  =E  0  0  »  E  < 1 1  =E  0  0 

00030970 

DR  =  HBF  FLORT  <NBDI VI 

0003094 0 

DO  11  I=2.NDIV 

00 03 OS 50 

II  =  1-1 

00030960 

Z  *'I.i  =  ZCIM  +  DZ 

00030S70 

R .  1 1  =  flan  +  DR 

00 0303 SO 

xi  <i>  =  r <: 1 1 

0  003  OS'?  0 

El  ai  =  EOO 

00030900 

F  <  I )  =  EOO 

0003091 0 

E<I1  =  EOO 

00 03 0930 

00 03 093 DC 

11 

CONTINUE 

00 03 094 OC 

. ..CRLCULRTE  FINRL  VOID  RRTIOS  FOR  DPEPGEP  FILL 

00030950 

DO  14  1=1 <NBDIV 

00030960 

SI  =  GC ♦ <ELL-Z  <I> 1 

000:30970 

IF  <S1  .LT.  0.0.1  SI..*  0.0 

00030950 

DO  1£  N=2. LDF 

00030990 

S£  =  SI  -  RS  <N> 

00031 000 

IF  <S2  .LE.  0.0)  GOTO  13 

00031 01 0 

1  s 

CONTINUE 

on 03 J 030 

EFIN  <11  =  ES  <LDF1  5  GOTO  14 

on 031 030 

13 

NN  =  N-l 

00031 040 

EFIN  <11  =  ES  <Ni  +  <’S£*  '.ES  <'NN)  -ES  <N>  )  /  <RS  <NN1  -RS  <H>  >  ) 

00031050 

14 

CDNT I HUE 

00  0:31  060 

00031 07 nc 

EFIN(NDIV)  =  EOO 

00031080C  . .  .CALCULATE  FINAL  VOIP  RATIOS  FOR  FOUNDATION 

00031090  IF  WEI  .  EQ.  2>  SOTO  20 

00031  100  15  Cl  =  ELL  1  ♦SC  1  5  C2  *  ELL»GC  ♦  01 

00031110  SI  =  Cl  +  C2 

00031120  DO  18  1=1 >NPIV1 

00031  130  S2  =  SI  -  Z 1  •' I  >  ♦•3C 1 

00031140  00  16  N=2.LBL 

00031150  S3  =  "2  -  RSI <N> 

00031160  IF  (S3  .LE.  0. o>  GOTO  1 7 

00031170  16  CONTINUE 

00031130  FFINlfn  =  ESI <LBL>  !  GOTO  18 

00031 190  17  NH  =  U-l 

00031800  EFIN1 <I>  =  ES1<M>  +  <S3* <ES1 (NN> -ES l <N, > ✓ <RS1 <NN ■ -RS 1 • N < > ) 

0 0 0 1210  IS  CON T I NlJE 

00031220  IF  •  NBL  .EO.  3>  ER<NMV1>  =  EFIN1  •  NPI  Vl"> 

00031 23 OC 

00031240C  ...CALCULATE  INITIAL  STRESSES  AND  FORE  PRESSURES 

00 03 l 25 uC  . FOP  FOUNDATION  LAYER 

00031860  DO  19  I=1,NDIV1 

00031270  IJ01  OI'*  =  GW  ♦  <WL  1  -X II  < I > > 

00031280  l.tl<I>  =  C2  -  00 

00031290  UM1  <I>  =  U01  Cl)  *  uia» 

00031  300  ffstpi  <  I  >  =  Cl  -  SC  1  *Z  1  •' I +  00 

00031310  T0STP1 CI>  =  EFSTR1 <l>  +  UW1 <l> 

00 03 1320  19  CONTINUE 

O0031330C  . ULTIMATE  SETTLEMENT  FOR  COMPRESSIBLE  FOUNDATION 

00031340  VRI 1  =  F INTI <NPI Vt ) 

00031350  CALL  INTGPL  <EFIN1 *  DZ1 » NDIV1 » FINT1T 

00031360  SF IN!  =  VRI1  -  FJHTI <NBIV1> 

00031370  IF  -NBL  .EC).  3)  GDTD  25 

00 031 38 OC 

000 31 39 OC  . FOR  DREDGED  FILL  LAYER 

00031400  20  DO  21  I=1»NDIV 

00031410  uoa>  =  GW  ♦  OWL -XI  d  >> 

00  031420  ua>  =  GC  ♦  (  ELL  -2 <  I  >  ' 

00031430  UW(I>  *  UOd>  *  U«D 

00031440  EFFSTR(I)  =  0.0 

00031450  TOTSTP  •:  IT  =  UW<I> 

00031460  21  CONTINUE 

00031470C  . ULTIMATE  SETTLEMENT  FOR  DREPGFD  FILL 

00031480  CALL  INTGPL  (EFIN> PZ* NDI V» FINT » 

00031490  SFIN  =  EOO*ELL  -  FHtT(NDIV> 

00031 500C 

0003151 OC  ...CALCULATE  FUNCTIONS  FDP  DREDGED  FILL 

0003152 OC  . PERMEABILITY  FUNCTION 

00031530  DO  22  I  =  1  > LDF 

00031540  PK  I >  =  PK  ( I  >  /  >"1 .  O+ES  <I>  > 

00031550  22  CONTINUF 


B9 


'Vi 


%  - 


000315600 
00 031 57 OC 
00031580 
000?  159  f* 
00051600 
0003161 0 

. SLOPE  OF  PERMEABILITY  FUNCTION  —  BETA 

. AMO  3 LOPE  OF  EFF  STRESS-VOID  RATIO  CURVE  —  DC BE 

CD  =  ES ■  2>  -  ES  C 1 > 

BETAC1>  =  <PK<£  >-PkCl  V>  ✓  CD 

DSDE  c  1  >  =  'PS  C2)  -PS  c  1 )  )  /  CD 

L  =  LDF  -  1 

< 

00031680 

00 03 163  0 
00031640 
0005 1 650 

DO  33  1=3. L 

11=1-1  5  IJ=I+1 

CD  =  tS'IJ>  -  ESCII) 

BETA<  I)  =  CPk  ■ -.1  J)-PK(in)  /  CD 

00031660 

00031670 

33 

DSDE* I)  =  CRS<IJ>-PS<II>>  ✓  CD 

CONTINUE 

00  03 

CD  =  ES  <LDF)  -  ES  CL) 

00031690 

BETA  <LDF'i  =  (Pk  CLDF>  -PK  CL)  >  CD 

) 

00031700 
0003171 OC 

DSDE  <LDF)  =  CRS<LDF)-RS<L)>  ✓  CD 
. PERMEABILITY  FUNCTION  TIMES  DSDE  —  ALPHA 

l 

00031730 

00031730 

DO  24  1  =  1,  LDF 

ALPHA (I)  =  PKCI)  ♦  DSDE <1 > 

i 

00031740 

34 

CONTINUE 

00031 750 

00  03 17600- 
no  03 177  Of: 

IF  CNBL  .EQ.  2)  GOTO  39 

...CALCULATE  FUNCTIONS  FOP  COMPRESSIBLE  FOUNDATION 

0003173 OC 
00031790 

?5 

. PERMEABILITY  FUNCTION 

DO  36  1=1, LBL 

j 

00031800 

00031610 

36 

PKia>  =  pki <i>  ✓  <1. 0+esi -:d) 

CONTINUE 

i 

00 031 83 OC 
000 31 83 OC 
00031840 

on 031 850 

. SLOPE  OF  PERMEABILITY  FUNCTION  —  BETA1 

. AND  -SLOPE  OF  EFF  STRESS-VOID  PAT  ID  CURVE  —  DSDE1 

CD  =  ES1<2>  -  ESI  Cl) 

BETA1  <1  >  =  CPK1  C£)-PK1  <1>)  /"CD 

i 

0003186 0 
00031870 
00031880 

0 Of 3 1890 

DSDE 1 1 1 >  =  CPS  1 C3) -RS 1  C 1  >  )  /  CD 

L  =  LBL  -  1 

DO  87  1=2, L 

11=1-1  ;  ij=i+i 

i 

00031900 

00031910 

0003)930 

00031930 

37 

CD  =  ESI  el J)  -  ESI  Cl  I) 

BETA 1 C I >  =  CPK 1 C I J > -PK 1 C 1 1 >  >  /  CD 

DSDE  1  •' I  >  =  CPS  1  <  I  J>  -PS  1  C 1 1 )  >  CD 

CONTINUE 

00031940 
00031950 
00031960 
00031 97 OC 

CD  =  ES1CLBL)  -  ESI  CL) 

BETA! (LED  =  CPK1 CLBL) -Pk 1  CL) >  ✓  CD 

DSDE1CLBL)  =  CRS1 CLBL) -RSI  CL) >  CD 

. PERMEABILITY  FUNCTION  TIMES  DSDE  —  ALPHA1 

i 

00031980 

00031990 

DO  38  1=1, LBL 

ALPHA 1  Cl)  =  PklCI)  ♦  DSDF1CI) 

00033000 

33 

CONTINUE 

0003301 OC 
000330? OC 
00035 040 
00033 050 

39 

. . . CALCULATE  BOTTOM  -BOUNDARY  DUD? 

DUD? 10  =  U1C1>  ✓  DUO 

IF  CNBL  .EG.  2)  DUDZ10  *  U<1)  ✓  DUO 

4 

0003806 OC 
00 033 07 OC 
00033080 
00038090C 

0  00331  one 
00033110 
00038180 
00033130C 
000331 4 OC 

...COMPUTE  VOID  RATIO  FUNCTION  FDR  INITIAL  VALUES 

CALL  VRFUNC 

RETURN 

END 

i 

BIO 


00114011110 

0004001  or 

O0G4OO£0C 
00  04  00  3  DC 
00040040C 
0  004  0  of,  OC 
000400*01' 
0004 0070 
00040030 

00 04 0040 
000401 00 
0004011 0 
00040120 
00 04 01 30 
000401 40 
0 0  0401 5  0 
000401*0 
00040170 
oo 04 01 SO 
00040190 
O0i14O?00C 
i'i  0  04  02  1  0 
0004 0?20 
0  0 114  0 2  3  0 
00040240 
0004 025 OC 

0  i'i  04  Of*  0 
0 0  04  02*5 
000402** 

0  0 04  0270 

00040220 

00040290 

00040300 

0004031 0 

00  04  0.32  OC 

00040320 

00040340 

00040350 

000403*0 

00040370 

00 04 03S0 

00040390 

0  0  04  04  0  0 

0  0 04  04 1  0 


SUBROUTINE  RESET 

♦  RESET  CREATES  PREVIOUS  CALCULATIONS  TO  HANDLE  ♦ 

♦  ADDITIONAL  DEPOSITIONS  OF  DREDGED  FILL.  ♦ 


COMMON  DA< DUG* DUDZ1 0. DUDZ1 1 , DUDZ21 . DZ.  DZ 1 . DO. EG- EOC > ELL* ELL  1 . 

S,  GC*  GC 1 • GS*  6S 1  * GSBL  > GSDF * Gl.i •  HBL  * HDF* HDF1  *  IN.  IMS.  TOUT* 

*.  I  OUTS .  L  BL  ■  LDF .  MT I  ME » H6D I V .  MED  I V I .  NBL .  ND  *  ND I V  •  ND I V 1 . 

$,  NFL  AG .  NM .  NPROE  >  NPT .  NND .  NNN  *  NT  I  ME .  PK  0 . 0  0 . 0 1  *  S  ETT  >  SETT  I  . 

SF I N . SF 1 N 1 . TRU  *  T I ME . TPR I NT , UC  DM . UCON 1  *  VP 1 1 .  ML .  ML 1 » 2K  0 . 

s  a  <1 oi)  * ai  a  n .af.ioii  *afi  an  .alpha. 511  .rlfhai  <511  . 

i:  BETA  <511  *  BETA1  <511  » BF  <  l  01  1  *  BF1  <  1 1 1  *  I'SDE  <5 1  1  *  DSPE 1  <511  * 

*■  E<1011  *E1  <1011  *E11  <111  *EFIN<1011  -EFINl  <1 11  *EP<111  * 

S:  ES>511  .ESI  <511  * EFFSTP <1  011  .EFSTP1  ‘  1 1 1  * F <1  011  . FI  <1 11  * 

&  FINT  .101)  .  F I  NT  1  <111*  PK  .511  *  PK  1  <511 .  PK  <5 1  >  .  RK1  <511  * 

t:  PS  <511  ?PS1  <511 ,  TOTSTP.  1011  >  TOSTP1  <1 11  .  U  <1  01  >  .  U1  <1 11 » 

s,  uo  <ioi>  *uoi  an  ,uwa  on  *u«i  <1  n  »xi  <1011  *xn  <1 11 . 

4-  z  on  .21  <111 

IF  .NBL  . EO .  £1  ML  =  ML  1 
IF  <HDF1  .LE.  0.01  RETURN 
IF  .NBL  .EG.  3>  01  =  HDF1  +  Q1 
IF  <NBL  .EO.  31  GOTO  5 

...CALCULATE  ELL  FOP  NEXT  DREDGED  FILL  LAYER  AND  RESET  CONSTANT 
EL  =  HDF1  .'1 .  O+EOO. 

IF  <NBL  .EQ.  21  U<11  =  U<11  +  EL*GC 

Ul<11  =  Ut<11  +  EL*6C 

NDZ  =  I  FIX  (EL-'DZI 

ELL  =  ELL  ■  +  DZ»FLOAT  <NDZ1 

NV  =  ND  +  1 

ND  =  ND  +  ND? 

NB  =  ND  -  1 

...CALCULATE  ADDITIONAL  COORDINATES  AND  SET  VOID  RATIOS 
DD  1  I =NV  *  ND 
II  =  1-1 

Z<I>  =  Z< III  +  DZ 
fl.I1  =  Adi)  +  DA 
XI <11  =  XI <111  +  DA 
El  <H  =  EOO 
F<11  =  EOO 
E  < 1 1  =  EOO 
1  CONTINUE 


Bll 


0  0  ft  4  04.?  00 
ft 004 043 ft 
00040440 
00 04 0450 
ft  0 04 04b 0 
00040470 
00  04  043  ft 
0004 0490 
00040500 
00040510 
00  04  05c  0 
0004 0550 
00040540 
0004055 OC 
0004 05b OC 
00 04 0570 
ft 004 055  ft 

00040590 
ftO  04  ftb  ft  ft 
ft  0  04061  0 
00040650 
00 04  Ob  50 
II II 04  ftb  4  ft 
0004065ft 
0004 0660 
ft 004 ftb 70 
00040650 
ft  ft  ft4  069  0 
000407005 
0004071 0 
ft  0  04 1175  0 
ft 004 07 5 OC 
00 04 074 OC 
01104  0745 
000411750 
00040755 
0004076  0 
0004077OC 
HO  04  075  OC 
00040790 
000405 00 
0004051 OC 
00 04 055 OC 
00040530 
00040840 
00040550 
00040860 
00040570 
mi  04  083  OC 
0  0  04  089 OC 
0 004 08 00 
0004091 0 
0004 095 OC 
00040950C 


.  .  . CAL  CULATE  FINAL  VOID  PATIOS  FOP  DREDGED  FILL 
00  4  I  =  1 <  MB 

51  =  GC*  vELL-Z <1)1 

IF  Cl  .LT.  0.01  S1  =  0.  0 
DO  8  N=8>  L  OF 

52  =  SI  -  PS <M) 

IF  ca  .LE.  0.0>  GOTO  3 
a  CONTINUE 

EF I N 'I<  =  ES'LDF)  ;  GOTO  4 

3  NN  =  H-l 

EF I N  <  J  >  =  ES<N>  +  C3*  fES  ''NN)  -E  S  iN)  s  fRS  iNN>  -PS  <N)  >  > 

4  CONTINUE 
EFINi.NO'  =  FOO 

...CALCULATE  FINAL  VOIP  PATIOS  fop  FOUNDATION 
IF  iNBL  .EG.  a>  GDTQ  9 

5  Cl  =  ELLl^GCl  ;  C£  =  ELL*GC  +  01 

si  =  ci  +  ca 

DO  8  I  =  l«  ND I V 1 

sa  =  si  -  21<i>*gci 

DO  6  N=Z»  LBL 

53  =  sa  -  PS1 <N> 

IF  C5  .LE.  O.o.)  GOTO  7 

6  CONTINUE 

EFIMl-I)  =  ESI  il.BL)  ?  GOTO  3 

7  NN  =  N-l 

EFIN1  <1  •  =  ESI  if-o  +  C 3*  vESl  >  NN)  -ES  1  «W  )  ✓  iRSl  iNN»  -RSI  <N)  >  > 

8  CONTINUE 

. ULTIMATE  SETTLEMENT  FOP  COMPRESSIBLE  FOUNDATION 

CALL  INTGRL <EFIN1 »  DZ1 >  NDIV1 »  FINT1) 

SF INI  =  Vp II  -  FINT1 iNDIVl) 

. . . PESET  BOTTOM  BOUNDARY  DUD2 
IF  iNBL  .EQ.  3)  Uia>  =  UU1)  +  HDF 1 
DUDSIO  =  U 1  •  1 )  DUO 

IF  iNBL  .  Et>.  3)  RETURN 

9  IF  iNBL  .  EG.  3>  IUJDZ10  =  Uil)  y  DUO 

. ultimate  setlement  fop  tdtal  dredged  fill 

CALL  INTGPL  iEFIN*  DZ»  ND»  FINT> 

SPIN  =  EOO*ELL  -  PINT 'NID 

...SET  VOID  PATIO  FUNCTIONS  FOP  RESET  VALUES 
N  =  NV-1 
DO  10  I *NV *  HD 
AF  1 1 )  =  AFiN) 

BF' I)  =  BFiN) 

10  CONTINUE 


RETURN 

END 


v 


SUBROUTINE  FPIFEG 


00050000 
0005001 ne 

0005002 nr 

Of.  1.500  7:  OC  «  FPIFEO  CALCULATES  NEW  VOIP  PATIOS  ft?  CONSOL  I  PAT  ION  PROCEEDS  ♦ 

0i.05004f,:  •  F,  AN  EXPLICIT  FINITE  DIFFERENCE  SCHEME  BASED  ON  PREVIOUS  ♦ 

00 050 05 OC  ♦  VOIP  PATIOS.  SOIL  PARAMETER  FUNCTIONS  ARE  CONSTANTLY  ♦ 

III. 050050,:  ♦  UPDATED  TO  CORRESPOND  WITH  CURRENT  VOID  RATIO.  ♦ 

00050070C  ♦♦♦♦♦♦ 

00  05  005  0.: 

00 050090  COMMON 

00  05  01  00  li 

00050110  4 

00  05  01 2  0 
0  0  05  0 1?0  :< 

00  0501  40 
00050150  * 

000501  A 0  t 

00 0501 70  % 

00 050 ISO  8 

00050190  & 

00  05  Or  00  :<• 

00050210  Sc 

0005022 OC 

00050236C  ...SET  CONSTANTS 

00050240  CF  =  TAUx<GW*P2> 

00050250  DZg  =  DZ-*2.  0 

000502*0  NNP  =  ND  -  1 

00050270  IF  (NFL  .EO.  2>  GOTO  5 

00 05  029  0  DZ12  =  DZ 1*2.0 

00 05 0290  CF1  =  TAU  <GW»DZ1 > 

0  005  0795  IF  <NBL  . EG.  3'<  ER<NDIV1>  =  EFIN1  CNDIVl > 

00050300C 

0005091 OC  . LOOP  THROUGH  FINITE  DIFFERENCE  EQUATIONS  UNTIL  PRINT  TIME 

0005  032 OC 

O0050330C  ...CALCULATE  VOID  RATIO  OF  IMAGE  POINT  AND  FIRST  PEAL  PD  I  NT 

00 05 034 OC  . FOP  COMPRESSIBLE  LAYER 

no 05 0350  1  DO  £  I =2*  LFL 

000503F.fi  Cl  =  ER  <  1  •  -  ESI  <I> 

010050370  IF  ‘O'  1  .  GE.  0.  0>  GOTO  3 

00 05 0330  2  CONTINUE 

00050390  DSED  =  PSDE1 <LBL)  ?  GOTO  4 

00050400  3  II  =  1-1 

0005041  0  DSED  =  DSDE1-I)  ♦  <C1* • DSDE1 ' I > -DSDE 1 < 1 1 > ) ' <ES 1 <I> -ES 1 ■  1 I) ) > 

00050420  4  F10  =  FI  <£>  +  DZ12*  <'GC  1 -t-PUDZl  1  >  xDSED 

00050430  DF  =  <F1<2>-F10>  *  2. 0 

00 05 0440  PF2DZ  =  <F1 <£>-£. 0*F1 <1>+F10>  x  DZ1 

00050450  AC  -  <AF1 <2>-AFl <1>>  x  DZ1 

000504*0  ER'1>  =  F 1  <  1  >  -  CF1*.:PF*<GC14BF1  a>+AC>+DF2DZ*AFl  -  l>i 

00050470  IF  -EP- 1>  .LT.  EFIN1<1>>  EP  <1>'  =  EFIN1  ■ 1> 

00050430  IF  <ER  <  1 P  .  GT.  E11<1>>  EP<1>  =  EU<1> 

00050490  IF  -MBl  .EG.  3)  GOTO  24 


DA * DU 0 . DUDZ  t  0 *  DUDZ 1 1 »  DUDZ21 *  DZ, TZ1 ,  DO*  EO*  EOO*  ELL*  ELL1 * 
GC >  GC  1 » GS * GS 1 , GSBL*  GSPF * GW?  HFL .  HPF » HDF 1. IN* INS* 1 01* T * 

I OUTS  * LBL i LDF  > MT I ME . NFD I V  > NBD IV 1  * NBL *  ND* ND I V *  ND IV 1 . 
NFLRG . NM* NPPOB *  NPT *  NND* NUN  *  NT I ME  *  PK 0  *  0  0  *  0 1  *  S  ETT . S  ETT 1  * 
SPIN*  SFIN1  *  TAU  *  TIME,  TPPINT.UCON.UCONl ,  VP  1 1  *  WL  >  *>IL  1  *  Zl'O. 
A  <1  01  >  *A1  (11)  ,AF<101>  ,AF1  <11>  » ALPHA  *5 1  >  *  ALPHA  1  -  51  ,  * 
BETA  -  51) , BETA1 <51) , BF • 1 01>  » BF1 <1 1 >  * DSDE <51 ■ . DSDE1 <51>  * 
EClOlJ.El  <101),  Ell  <11>  .EFIN'101) .EFIN1  < 1 1 > , EP  <1 1- • 
ES<51  ■  .ESI  <5 1  >  ,  tFFSTR  >101 >  -EFSTPl  ■  1 1  >  ,  F  .;  1  01  '<  *  F 1  -  1 1  •  * 
FINT-  lOn  * FINT 1  •.  ll.<» PK  <5  0  •  PK  1  <51>  ,RK<5P  *  PI  1  <51  > , 
PS<51> *  RS 1 <51 > *  TOTSTR 1  1  OP . T0STP1 ■ 1 1 > * U ■ 1 OD * U1 • 1 1 • * 
uo< too  .uoi  aic  lUwaoi)  .uwi  <n>  ,xi  aoo  ,xi  p  1  p  , 
zaoo  »zi  <u ' 


B13 


no  05  05  one 
no  05  os  1 0 

00050520 

o  o  0*5  or- }  o 
rocfor-40 

00050^50 
0005050-0 
00050570 
00050590 
00050590 
n 005 Oft 00 
00  05  Oft  1  0 
0 0050690 
00  05  Oft  3  i'i 
00  05  Oft  4  i'i  C 
000506-500 
00 05  Oft ft  0 
00 0 5 Oft 70 
00050630 
00050690 
00050700 
0005071 0 
00 05 0720 
00050730 
00050740 
00 05 0750 
Hu  05  07ft  0 

00 05 07 70 
0  0  05  07?-  0 
no  05 07? n 
0  0  05  Oft  0  0 
00050? 1 0 
00 05 Oft 20 
nVi  05  0.9?  o 
00 05 09400 
i'i  0  M  5  Oft  5  I'i  C 

00  05  Oftt  o 
00050870 
0  Oi  05  0930 
00  05  0:3?  0 
00 05 0900 
O0050?l 0 
00050920 
0005 0930 
00050941' 

0  0  05  035 n 
00 050aft0 
0005 0970 
00050990 
00050990 
00 051 000 
00051010 

0005-1  02  0 
00051  o?.n 
00051 040 
no 051 050 
00051 OftO 
00051 070 
00051 090 
00  05-1  09  0 
00051 i one 


. F DC  DREDGED  FILL 

5  DO  6-  I  =2  .  LDF 
Cl  =  E  *  1 >  -  ESd) 

IF  ('L  l  .  GE .  0.  0<  GOTO  7 
6-  CONTINl-F 

DC ED  =  DSDElLDF-  !  GOTO  3 

7  II  =  1-1 

DSED  =  DCPE'I)  *  .  Cl  ♦  <  DC  DE  (I  >  -DSDE  <  1 1  >  >  ■ F7  <  I  >  -EC  <  1 1  ■  ■  > 

8  FO  =  F  '-ft)  ♦  D22*  •  GC+DUD221>  DSED 

DF  =  < F  ■'  ft)  -F 0>  2.0 

DF2DZ  =  >'F  ■  2>  -ft .  0*F  <  1  ">  +F 0  >  DZ 

8C  =  (HF(?)-hF(1>)  V  IC 

E«l;i  =  F  ( 1 )  -  CF*1  DF*  «GC*BF  <  19  +60  +DF2D2*l4F  '  1  >  * 

IF  •  E  (  1  •  .LT.  EF  IN<  1  •  ">  E(l>  =  EFINd) 

...C8LCUL6TE  VOID  P6TI0  DF  TOP  POINT  IN  COMPRESSIBLE  L6YEP 
IF  CNBL  .EC.  2)  GOTO  27 
DO  9  1=2* LDF 

Cl  =  £(l>  -  EC  <I> 

IF  (C  1  .GE.  0.0>  GOTO  10 

?  continue 

EST  =  PC (LDF)  5  GOTO  11 

10  II  =  1-1 

e't  =  rs 1 1 :■  +  'Ci*|Pv,:  p-psai>>6(FCi'i>-E?  ai):u 

11  DECT  =  EST  -  EFFSTP  <1 > 

LIT  =  Ud>  -  DECT 

EFS1  =  EFSTPl 6NDIV1 •  +  DECT 
DO  12  1=2. LEL 

Cl  =  EF 3  1  -  PC  Id) 

IF  'Cl  .LE.  0.0)  GOTO  13 

12  CONTINUE 

EP  (ND I V 1 >  =  EC1 (LEL1  i  GOTO  14 

13  II  =  1-1 

EP(NDIVl)  =  ESKI)  +  (Cl* (ESI  d I) -ESI  <I)  .)  6  (RSI  (ID  -RSI  <I>  )  ) 

. ..RESET  P0UND8RY  DUTC  FOF  DREDGED  FILL 

14  DO  15  1=2. LEL 

Cl  =  EP'NBDIVO  -  ESI  (I) 

IF  (Cl  .GE.  0.0)  GOTO  16 

15  CONTINUE 

EST1 =  PCI  (LEL 1  5  GOTO  17 

16  II  =  1-1 

EST  1  =  PS1 <I>  +  (Cl ♦(PCI (f)-RSl (II>)y<E51 « I) -ESI  1  II) >  > 

17  UT 1  =  Ul(NEDIVl)  -  EST1  +  EFSTP1 (NBDIV1 > 

PUD2 1 ft  =  (UT  -  UT 1 )  D21 

DO  18  I =2. LEL 

Cl  =  EP(NDIVl)  -  EC  1 < I > 

IF  (Cl  .GE.  0.0*  GOTO  1? 

18  CONTINUE 

REFER  =  PF 1 'LEL 1  ?  GOTO  20 
1?  II  =  1-1 

REFER  =  PF 1  •  I  •  ♦  <C  1  ♦  ( PF 1  ■  I )  -EF  1  (  II  "•  )  /  (ES 1  <  I .)  -ES 1  1  II)  >  > 

20  DO  ftl  1=2. LDF 

Cl  =  E(l>  -  ESd ) 

IF  (Cl  .GE.  0.0)  GDTD  22 

21  CONTINUE 

EKE  =  PF  (LDF)  ?  GDTD  23 

22  II  =  I-t 

PFE  =  PF  ( I  >  +  (CM  (PF  •'  I)  -FF  ( 1 1>  )  '  (ES  (  I  (  -ES  (II))) 

2?  TIUP721  =  DUD 2 1 2  ♦  PPV  ER  6  PKE 


B14 


* 


no 05 1 1 1  oc 

. .  .CRLCULRTE  MEM  VOIP  RRTIOS  FOR  REMRINDEP  OF  MRTEP1RL 

0  0  05 1  l  8  or 

. in  COMPRESS I BLE  FOUUDRTlON 

00 051 1 ?o 

£4 

on  85  I =8 *  NBDI V 1 

00051 140 

II  =  1-1  :  IJ  =  1+1 

00051 150 

OF  =  • F 1  >  I  J  ■  -FI  'll  .  »  £ .  0 

00051160 

DF8D2  =  <F  1  1 1  J>  -F  t  *  I  •  ♦? .  0+F 1  ■ :  1 1>  >  02 1 

00051  1  70 

ftC  =  ‘RF1  ■  I J  >  -FtF  1  <  1 1  o  x  Iiz  1 8 

00051 . 50 

EP  ■  I  (  =  Fia>  -  CF1*‘DF*<GC1*BF1  •  I>+RC.:'+0FED2*F|F1  (I>> 

00051 1 9 0 

5C( 

CONTINUE 

0005 1 60  or 

. RESET  FDR  NEXT  LOOP 

00051 8 1 0 

DD  86  1=1 • MDIV1 

0  0  05 1 32 0 

F 1  •  I  >  =  EP  ’'I  j 

00051c 30 

4?  6 

CDNT I  HUE 

00051 840 

IF  <NBL  .EO.  3>  GOTO  30 

ft  n  05 1  c  *5  OC 

00051?6nr 

...MEM  VOID  RRTIOS  IM  OREOGEO  FILL 

00051 27 M 

£7 

00  8"?  1=8*  NMD 

0  0  05  1  8  6  o 

II  =  1-1  5  IJ  =  1+1 

00  051  c'Sil 

OF  =  '  F  ■  I  J)  -F  •:  H  >  )  /  3.0 

0  0  05 1  3  0  0 

OF 8 02  =  (F(I  J>-F<IJ*3.  0+FtII>>  D2 

0005 1? 1 0 

RC  =  'RF ■: I J> -RF •II')  '  028 

00051  ?*?  0 

E>.  I  *  =  Ftl>  -  CF*  <  OF*  <:GC*BF  ( I  >  +RC  >  +0F3B2*RF  <  JJ 

00  05)  ";3Q 

29 

continue 

0  0  051  ;.40lT 

. RESET  FOR  NEXT  LDDF 

00051350 

00  89  1  =  1*  UNO 

0  0  05 1 36  0 

FiT)  =  E • I > 

00 05 1370 

39 

CONTINUE 

00  051 3'rfOC 

0005 1  "  =•  ilC 

.  . . RESET  BOTTOM  BOUMORRY  RW2  FOP  COMPRESSIBLE  LAYER 

0005 1400 

IF  (NEL  . EO.  8>  GOTO  34 

0005141 0 

30 

DC  31  I =8  *  LBL 

0005143  0 

Cl  =  ERd)  -  ESI  d  > 

00051430 

IF  'Cl  .GE.  0. 0>  GOTD  33 

0005144H 

31 

CONTINUE 

no 05 1450 

PF'KER  =  PCI  'LBL- 

0  0  05 1460 

E3T1  =  RSI- LBL-  ?  GOTO  33 

00051470 

32 

II  =  1-1 

00051480 

CS  =  Cl  <ESl  d>-ESl  < 1 1 >  > 

0005 l 480 

PPKER  =  PK 1  < I >  ♦  CS*  < F'K  1  < I  >  -PB 1  ■  1I>1 

00051500 

EST  1  =  PS  1  d  >  +  €£♦  <RS  1  <  I J  -RS  t  d I  ■'  > 

00051510 

33 

01*02 1 1  =  DUO?  10  ♦  PKO  RPKEP 

00051580 

UTl  =  Him  -  EST  1  +  EFSTRim 

00051530 

DUD? 10  =  UTl-  •  DUO 

00 05 1540 

GOTO  38 

0  0  05 1 55  OC 

I 

! 


f 


I 


00051560C  ...RESET  BOTTOM  BOUNDARY  DUDZ  FOP  DREDGFH  FILL 

00051570  34  DO  35  1=3. LBF 

00^51530  Cl  =  E  >1 ■  -  ESd) 

00051530  IF  <C 1  . GE.  0. 0>  GOTO  36 

0  0  05 1  6  0  0  35  CONT  I  MlJF 

00051610  PKE  =  PFiLDF) 

00051630  EST  =  PS>LDF9  5  GOTO  37 

0  0  05 1 6  3  0  36  II  =  1-1 

00051640  C3  =  Cl  -  fEJiO-Eidl)) 

0  0  051  650  PKE  =  PK  <  I  9  *  C£»<FK  •  I )  -PK  LI  I  *  9 

00051**0  EST  =  PS  09  +  C3»  < PS  >  I  9  -PS  <  1 1  >  9 

00051670  37  BUDZ81  =  DUDZ  10  ♦  PKO  ■'  PK E 

00 05 1 6. SO  UT  =  UO)  -  EST  +  EFFSTP  0  9 

00051690  DUBZ10  =  UT  DUO 

>>0051  700C 

00051 71 OC  ...CALCULATE  ALPHA  FIND  BET  ft  FOP  CURRENT  VOID  RATIOS 

00051 73 0  33  CALL  VPFUNC 

0005] 73 OC 

00051 7407  ...CALCULATE  CUPPENT  TIME  AND  CHECK.  AGAINST  PRINT  TIME 

00051750  TIME  =  TAU  ♦  FLOAT -'NNN  > 


00051 760 

NNN 

=  NNN  +  1 

00051770 

IF 

(TIMF  .LT. 

TPPINT  .AND.  NBL  .EC. 

GOTO 

1 

0  0  05 1 73  0 

IF 

•TIME  .LT. 

TPPINT  .AND.  NBL  .EC. 

GOTO 

er 

00051790 
0005 1 SOOC 

IF 

(TIME  .LT. 

TPPINT  .AND.  NBL.  .EC. 

3> 

GOTO 

l 

0  0  05 131 OC 

,  ,  ,  1 

CHECK  STABILITY  AND  CONSISTENCY 

00051330 

IF 

•NFL  .EO. 

3>  GOTO  39 

00051330 

STAB  =  ABS  O  B 

Z 1  ♦♦3*6W9  '  ‘c.  0*AF  1  (1  >  ;> 

00051340 

IF 

(STAB  .LT. 

TAIJ9  UP ITE  (I OUT.  1  0 09 

NPPOB 

0  ft 05 1 95  0 

CON 

S  =  ABS  < (3 

.  0*AF1  Cl 9  »  <GC  1  ♦BF  1  (19 

•  'I 

0  0  05 1 36 0 

IF 

(CONS  .LE. 

DZ 1 9  WRITE* I OUT.  101) 

NPPOB 

00051370 

IF 

(NBL  . EO. 

39  RETURN 

000516:3:0  39  STAB  =  ABS  ■  • BZ»»3*6W9  -  <8.  0*AF  <1)  »  ) 

00051390  IF  (STAB  .LT.  TAU9  WRITE <IOUT» 1089  NPPOB 

00051900  CONS  =  ABS0  3.  0*AFO9>/(GC*BFO9  9  > 

00051910  IF  (CONS  .LE.  DZ>  WRITE (I OUT. 1 039  NPPOB 

00 05 198 PC 

0 0  05 1.9  3 OC  .  .  .  FORMATS 

00051940  100  FORMAT  L s//s '38MSTABILI TV  EPPDP  — FOUNDATION  — PROBLEM* 159 
00051  950  101  FORMAT  <✓  v  4 OHCONS I STENCV  ERROR  — FOUNDATION  — PROBLEM*  159 
00051960  10?  FORMAT ■ -  "  ^4 OHSTAB I L I TY  EPPDP  —DREDGED  FILL  — PPDBLEM- 1^9 

00  05 1970  1  03  FORMAT  (,-v- -48HC0NSI STENCY  EPPDP  —DREDGED  FILL  —PROBLEM.  15 

00 051 93 OC 
0  o  05 199  OC 

00053000  RETURN 

00053  0 1 0  END 

0  0  052  03  OC 
000580'*  OC 
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AO-All*  112  ARMY  ENGINEER  WATERWAYS  EXPERIMENT  STATION  VICKSBURO— ETC  F/0  8/13 
CONSOLIDATION  Of  SOFT  LAYERS  BY  FINITE  STRAIN  ANALYSIS. <U> 

MAR  82  K  0  CARGILL 

UNCLASSIFIED  WES/MP/GL-82-3  ML 


00 06 00 00 
0006001  Of 
001**002013 
00  06008013 

000600400 
000600500 
00  Or.  006  Of 
00  03  0  070 
00 06 0080 
0  0  «6  i'  090 
000801 00 
000801 1 0 
00080180 
00 118  01 3  0 

00080140 
no 08 0 i 5 0 
00080160 
00 06 01 70 
00060180 
000801 90 
0006  0800C 
00 08 08 1  0 
0006 088 OC 
00  06  0?30 
00080840 
00060850 
00060560 
00  06  Or.  70 
00060890 
00060:90 
00060300 
00 06 031 0 
0 006 0330 

00 06 03 30 

00080840 
0006035  0 
00 06 086 OC 
0006 03700 
000603:30 
noo6  0.390 
00 06 04 00 
0006041 0 
no 06 0420 
00 06 04 30 
00060440 
00060450 
00060460 
00060470 
00 06 0480 
no 06 0490 
0006 0500 
0006051 0 
0  0  06  052  or 
00 06 05 30C 
00 06 0540 
no  06  055.0 
00 06 056 OC 
0006057013 


SUBS OUT IMF  VRFIJNC 


♦  VPFUNC  CALCULATES  ALFhA  AND  BETA  FUNCTIONS  • 

♦  FOP  CURRENT  VP I  0  PAT  IQS.  « 


COMMON  DA «  DU  0.  OUOZ 1  0  •  OUTS  1 1 . 00022 1,02-  02 1 »  DO  -  E  0 .  F  0  0  -  ELl  *  ELL  1 
A  GC»  6C1 < 6S  *  SSI  * GSBL<  GSOF-  Gi«<  HPL<  HOF  *  HOF  1  -  IN-  IMS  •  IDUT. 

*  I  OUTS  -  LBL  -  L  OF  f  AT  I  ME  -  NBD  !  V .  M  80 1  1 «  MBL  -  NO  -  NO  I V  •  NO  I V 1 , 
NR.  AS  -  NM ,  HPRQB  -  NPT  >  NND- N*‘l  U  NT  I  ME  *  PK  0  •  0  0  -  0 1  »  S  ET  T  -  3  ETT 1 
SFIN. SFIN1 . TAU<  TIME  - TPPINT<UC0N<LIC0N1<  VP! 1 » ML  •  ML  1 « 2* 0 

S  A  -  1  OP  <  A1  -'1  1>  <  AF  <1  01)  -  AF1  ‘ID  <  ALPHA  (51 )  <  ALPHA  1  <51*  < 

*-  BETA <51  >  .BETA1  <51  >  •  PF  <  1  OD.BFl  <1  1)  .OSOE  <51'  ,05 PEI  <51* 

•i  E < 1 01*  »E1  <1  On  <  El  1  <1 1>  »FFIN<:  101'  < EF INI  <11)  ,  EP  <  1 1  > . 

'i  ES  <51 )  « ESI  <51)  -  EFFSTP  <101)  .EFSTP1  <1 1)  » F  <1 01)  <  Fl<  1 1 '  . 

*  FI  NT  <1  01>  <FINTi  an  >PK<51)  »PK1  <5n  ,PK<51*  -PL  1  <51  '  * 

S:  PS  <51>  < RSI  <51  >  .  TQTSTPClon  »  TQSTF1  <  !  1)  < LI <101)  <U1  <1 1  >  < 

uo<ion  »uoi  <1 1  >  *  l*.<»  <101  < .  uw  1  <in ,  >: i  a  on  <vi  1  <  1 1  •  < 

•»-  zaon  *zi  <n:> 

IF  <NBL  .EQ.  8)  SOTO  4 

...FDP  COMPRESS I BLE  FOUNDATION 

DO  3  1  =  1  *  N  0 1 V 1 

OC  1  N=8<  LBL 

Cl  =  EP<n  -  ESI  CN) 

IF  <Ct  .  GE.  0.0'  GOTO  2 

1  CONTINUE 

AF1 <I>  =  ALPHA1<LIL> 

PF 1 < I '  =  BETAS  <LBL>  ?  GOTO  3 

2  NN  =  N-l 

CM  =  Cl  <ES1  CN'-ESl  <NM.» ) 

AF  1  <  I )  =  ALPHA  1  <N'<  +  CM*  <ALPHA  1  <N) -ALPHA  1<NN  •  > 

BF 1 < I >  =  BETA1 <N>  +  CM*  <EETA1 <N) -BETA1 <NN> ) 

3  CONTINUE 

IF  <NBl  . EQ.  3'  PETUPN 

...FOP  DREDGED  fill 

4  DC  7  I » 1. NND 

00  5  H=2<  LDF 

Cl  =  E<I)  -  ES <N' 

IF  <C1  .GE.  0.0)  GDTD  6 

5  CONTINUE 

AF  < I )  =  ALPHA < LDF) 

BF  <  I )  =  BETA-:  LDF''  !  GOTO  7 

6  UN  =  N-l 

CM  «  Ct  <E3  'N>  -ES  <MN  ' ) 

AF  < I )  =  ALPHA  <N)  *  CM*  <Al.PHA  <N) -ALPHA  <NN)  ) 

PF  < I *  =  BETA  <N>  +  CM* > BETA <N) -BETA <NN>) 

7  CONTINUE 

AF  <M0>  =  Al  PHh < 1 > 

BF  <,ND>  =  BETA  <  1) 


PETUPN 

END 


SUBROUTINE  STRESS 


00070000 

i'i  0  fi  7  0  i!  1  fiC 

nnnrnnEnr 

iifii-i7n>i-.r,r 

000700400 

I'i  I'I  II 7  Mil'1;  In: 

no  07004.  iy; 

i  •  0  <■ 7  00  7  nc 
OP  0700:50 

ii  0  07  0  09  0 

no 07 oi on 

000701 1 o 
0 u 07012 0 
r  on 7 oi so 

0  0  0 7014  0 
fin  07  01  SO 

000701 4 0 
00 07 01 70 
000701 80 
00 07 01 90 

no 07 or on 
0  0  07  08 1  in' 
no 07 088  or 

no  07093 0 

00070840 

0  007  OS'S  i'i 

on 07094  0 
nori7n?70f 
000702808 
00 07 089 OC 
00070300 
00070*  1  0 
00070820 
00070  -.30 
00070340 
00 07 08 SO 
00070380 
00  07  087  0 
00070380 
00070390 
00070400 
000704 1 0 
0007048  0 
00 07 04 30 
00070440 
00  07  04 S 0 
00070460 
00070470 
00070480 
00 07 049 0 
00 07 0500C 


♦  STRESS  CALCULATES  EFFECTIVE  STRESSES.  TOTAL  ^TRESSES.  ♦ 

♦  AND  POPE  WATER  PRESSURES  BA 'ED  ON  CURRENT  VOID  RATIO  * 

♦  AND  VOID  PATIO  INTEGRAL.  ♦ 


COMMON  DA .  DU 0 ,  D'JDZ  1 0 .  DUDZ 1 1 ,  PUDZE  1 .  DZ.  DZ 1 .  DO.  E0»  E 00 ,  ELL .  ELL  1 » 
A  GC ? GC 1 « GS. GS 1 ?  GS  BL • G SDF ,  GW, HBL  ,  HPF  >  HPF 1  *  IN  > I N  S , I  OUT • 

4.  I  OUTS  1 LBL  «LDF. MTIME?  NBDI V.  NBD I VI .  UBL , ND , ND I V .  NP I V 1 , 

S-  NFL  AG ,  NM .  NPROB  *  NPT .  NND » NHM , NT  I  ME  *  FT  0,00.01.  -SETT .  S  ETT  1 , 

A  SPIN,  SFIHt .  TALI.  TIME.  TPRINT. UC ON » UC  ON  1 .  VP  1 1 ,  WL  » WL  1 ,  ZK 0, 

4.  A  •'  1  0 1  •  ,  A 1  ( 1 1  '  ,  AF  <101)  ?  AF 1  <  1 1  > .ALPHA <5 1>  ,  ALPHA  1  <51  >  , 

$,  BE TA  <5 1  >  ,  BETA  1  <5 1 > • BF • 1 0  D . EF 1 < 1 1  >  •  DS DE  <51>  •  PSDE1  <51>  • 

A  e  <10D  .El  <1  OH  ,  Ell  <1 1.)  .EFIN<1  01>  .EFIN1  <11>  »EP<11>  . 

A  ES  <51>  .ESI  <5D  « EFFSTP  <1  01  >  ,  EFSTR 1  <11'.  .F.  101  ,F1  <11-  . 

A  F INT  <1  01  >  , FINT 1  <11  >  ,  PK  <51  >  ,PK1  <51  •  »Pk  <51  ■  .  FT  1  <5D  - 

i,  PS  <51>  »RS1  <51  > ,  TOTSTP  <1 01  >  ,  TQSTR1  <11  >  ,U<101>  ,  U 1  <  ID  , 

g.  uociod  »uoi  <1 1:>  ,ui.i<i  oi:» ,  lm.ii  1  <n>  .xi  <101 » ,xi  1  <  1 1> , 

Z<101>  ,Z1  <1D 

...CALCULATE  VOID  PATIO  INTEGRAL 
IF  <N Bl.  .EO.  3>  GOTO  1 
CALL  INTGPL  <E.  DZ,  ND.  R'lNT) 

IF  <NBl.  .EO.  2>  GOTO  7 
1  CALL  I NT GRL  <ER • DZI , ND I V 1 » F I NT1 > 

...FOR  COMPRESSIBLE  FDUNDAT I ON 

. CALCULATE  XI  COORDINATES  AND  STRESSES 

DO  8  I*l»Nfim 

XII  <I>  =  Z1  <I  >  +  FINT1CD 

8  CONTINUE 

WL  =  1.1.I L 1  -  NI1<NDIV1> 

SI  =  ELL  •  GC  +  01 
Ml  =  FIMT1  <NDIV1>  +  WL 
DO  6  I=1,NPIV1 
DO  3  N=2« LBL 

C1  =  EP • D  -  ESI <N> 

IF  <C1  .GE.  0.0)  GOTO  4 

9  CONTINUE 

EFSTP1CD  =  RS  I  <1  BL>  !  GOTO  5 

4  NN  =  N— 1 

EFSTR 1 <D  —  RSI <N>  +  <C1* <RS1 <N> -RSI <NN> ' ✓ <ES1 <N> -ES 1 • NN> ) > 

5  U  0 1 <  D  =  GW  ♦  <WL 1 -X 1 1 •  I  > > 

TO S  TP  1  <  I ')  =  GW*  I.  lit  1  -F I  NT  1  <  I  >  >  +  GS  1  ♦  CELL  1  -Z1  <  I )  3  +  G 1 
UW1<I>  =  TO? TP!  CD  -  EFSTR 1<I> 

Ul  <I>  =  UW1  <I>  -  U01  <D 

6  CONTINUE 

IF  <NPL  .EO.  3>  GOTO  13 
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<1007051  OC  ...FOP  DREDGED  FILL 

0007052 OC  . CALCULATE  XI  COORDINATES  AND  STRESSES 

00070530  7  DO  3  I=1.ND 

00070540  XI  a?  =  2<I>  +  FINT(I) 

00070550  8  CONTINUE 

00 0? 0580  WLL  *  ML  -  XI <ND> 

00070570  Ml  =  FI  NT  »ND>  WLL 

00070580  DO  1?  I  =  1 » ND 

000705*0  DO  9  N=2»LDF 

00070600  Cl  =  E < I >  -ES <N> 

00070610  IF  fCl  .GE.  0.0)  GOTO  10 

00070620  9  CONTINUE 

00070630  EFFSTR ■' I>  =  PSCLDFi  !  GOTO  11 

00070640  10  NN  =  N-l 

00070650  EFFSTR  ‘'I>  =  RS  <N>  +  >C1*  <RS  <N1  -PS  <NN1 )  s  (ES  <N  »  -ES  *'NN1 

00070660  11  U0-'I>  =  GW  ♦  <WL-XI«>> 

00070670  TOTSTP  <l>  =  GW* XW1-FINT  f  II  >  +  GS* < ELL-2 <I>> 

00070630  UW<I>  =  TDTSTR <I>  -  EFFSTR <I> 

00070690  Ua>  =  UW«I>  -  U0<I> 

00070700  1 ?  CONT I NUE 

0007071 OC 

OO070720C  ...CALCULATE  SETTLEMENT  AND  DEGREE  DF  CONSOLIDATION 
00070730  IF  tNBL  .£<?.  2>  GOTO  14 

00070740  13  SETT  1  =  At  t'NDIVl)  -  XI1<TNDIV1> 

00070750  UCDN1  =  SETT1  '  SFIN1 

00070760  IF  <NBL  .  EO.  3>  RETURN 

00070770  14  SETT  =  A <ND  •  -  XI <ND> 

00070780  UCON  =  SETT  ✓  SFIN 

00070790C 

0D070S00C 

00070810  RETURN 

00070320  END 

00070830C 
00 07 084 OC 


SUBROUTINE  INTGRL 'E » DZ» N« F> 


0008 00 00 
ados  oi'u  oc 
00 080 03 OC 

n003003nc  ♦  INTGPL  EVALUATES  THE  VOID  PATIO  INTEGRAL  TO  ♦ 

000@On40C  *  EACH  MESH  POINT  IN  THE  MATERIAL.  ♦ 

r,  nos 005  m: 

000300300 

00030070  DIMENSION  E < 1 01) j F <1 01.) 

nnnsonsnr  ...BY  SIMPSONS  PULE  FOP  ALL  ODD  NUMBERED  MESH  POINTS 

no os  0030  F<1>‘  =  0.0 

non? 01  no  DO  1  1=3. N. 3 

ftilfic 01 1  ft  p<I>  =  Ff  l-g>  +  D2*<£*I-3>+4.  0*E<I-1>+E^I>>^3.  0 

non-- ni  jo  1  continue 

00080!  JC.C  ...BY  SIMPSONS  3^8  PULE  POP  EVEN  NUMBERED  MESH  POINTS 

00080140  DO  3  I =4. N«  3 

0.0030150  Fei>  -  F<I-3>  +  DZ*  <E*  1-3)  +3.  04<ECI-3>  +E  C I  —  I  >  >  +E  •  I)  >  ♦  1  3.  0'8.  O' 

000301 80  3  CONTINUE 

000301 70C  ...BY  DIFFERENCES  FDR  FIRST  INTERVAL 

00 03 0 1 3 0  F£  =  DZ*  '-E  <£)  +4 .  0>E  f  3>  +E  <4>  )  .'3 .  0 

00030130  Fes')  =  Fe4>  -P£ 

0005C3 OOC 
000:3031  OC 

00050330  RETURN 

00030330  END 

0003054 OC 
0003035 OC 
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...... 


run’ll  on  on  SURPCU  aNF  DATOUT 

000900) Of 

'>0090030:7  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦•♦♦♦♦♦♦♦♦ 

no 0900? or  ♦  BATDUT  PRINTS  RESULTS  DF  CQNSOl  I  PAT I ON  CALCULATIONS  AND  ♦ 

0009004  0.;  ♦  BASE  I.ATA  in  tabular  fdpm.  ♦ 

00090050C  ♦♦♦♦«♦ 

0 009 00b 00 

00090070  COWON 

00090030  $. 

00  09  0090  i, 

o 0090  j  no  ?. 

000*011 0  4 

00090130  & 

00090130  i. 

000901 40 
00090150 
00090160  t- 

000*0170  t 

00090)30  i. 

00090190  i- 

00 0*030 or 

00090c  IOC  ...PRINT  CONDITIONS  IN  COMPRESSIBLE  FOUNDATION 

0  00*033 0  IF  <NEL  .EC!.  3>  GOTO  3 

00090? 30  IF  < NFL AC  .EO.  1>  WRITE 'I OUT. 100* 

000*034  0  IF  ''NFL  AG  .EC.  O'.  WRITE  <IOUT»  1  OS) 

0 0 09 035 0  WRITE*: IOUT.  101> 

00090360  WRITE  <IOUT» 1 03> 

00 09 0370  DO  1  1=1 . NLIV1 

000*0380  J  =  NB IV 1+1 -I 

000*0390  WRITE  a  OUT . 1 03’  Al < J> • XI 1 < J>  » 21 <J> .Ell <J> >  EP  <J> . EFIN1 <J> 

00090300  1  CONTINUE 

00090310  WRITE -  IOUT.  104> 

0  0  09  033 0  MR I TE  *  I OUT .105) 

000903 SO  DO  3  1  =  t . NDI VI 

00090340  J  =  NPIV1+1-I 

00090350  WRITE  aOUT.  1 C<3>  XII  <J>  >  T05TP1  <J> .EFSTPl  <J>  »UW1  ..J..IJOI  <J>  ,U1  <J> 

00 09 01 60  3  CONTINUE 

00090970  WPITE  a  OUT . 1 07>  TIMF.UCDN1 

00  0*90380  WPITE  CIOUT .  1 1  0>  SETTI.SFIN1 

00090390  WRITEaOUT.Hl>  BUPZ11 

00090400  WPITE  <  IOUT.  U2>  WLt 

00090410  IF  <NEL  .EC.  3)  RETUPN 

00 03 043 OC 


DA • DU 0 . DUD2 1 0 • BUDZ 1 1 . DUDZ3 1 . D2 . DZ 1 . DC . E  0 • E  0  0 « ELL . ELL 1 . 
GC  *  GC 1  *  GS. 6S1 .  GSBi. .  6S  I*F .  GW.  HBL .  HDF .  HPF 1.  IN.  INS.  I  OUT . 
IOUTS.  LBL.  LDF >  MTINE.  NBDI'-. NBDIV1 » NBL«  ND»  NBI V*  MM  VI » 

NFL AG . NM . NPPDB. NPT . NND . NNN. NTIME.PKO. 00. 01 »  JETT. JETT1. 
SFIN.SFINl . TAU .TINE. TPPINT • UCDN? UCDN 1 » VP 1 1 . ML . WL 1 . ZK  0 » 
A  < 1 0 1 >  »  A 1 < 1 1 >  »  AF  < 1 0 1 > • AF 1 <1 1> . ALPHA <51> . ALPHA 1  <51> . 
BETA  *51  >  ,  BETA1  <51  >  .  BF  <  1 01* .  BF1 . 1 1>  ,  DSI-E  <51  >  .  P7DF1  <51>  . 
£  <  1  01  >  .El  <101  >  .Ell  *  1 1  > .  EF I N  *  1  0 1  >  .  EF IN  1  *1  \>  « EP  < 1 1 >  > 

ES  <5 1 > . E  SI <5 1 > . EFFSTP  < 1 01 > . FFSTP 1 < 1 1 > . F  < 1 0 1 > . F 1 < 1 1 > . 
PINT <101 > .FINT 1 <11> .PK  <51> .PK1 <51 » .PK*51 > . Rk 1 <51 > . 

PS  <51  >  . PS1  <51  >  .  TDTSTP  <101...  TOSTR 1  <  1 1  >  .  U  <  1 0 1  >  .  U 1  <  1 1  >  . 

uo<ioi>  .uoi  <1  r>  »uw<ioi>  »uwi  • 1 1>  »xi  a  on .  x 1 1  a  r> , 

2  <  1  0 1  >  »  Z 1  <  1 1  > 


000904300  .  .  . PRINT  CONDITION?  IN  DREDGED  PILL 

1.0090440  3  IF  • NFLAG  .  EQ.  1>  WRITE  ‘'I  OUT .  106) 

00090450  IF  'NFLAS  .EO.  0)  WRITE  <1  OUT.  1 09) 

in;. 03 <1460  WRITE<  IOUT-  1  01  > 

ft  f*  03  04 7  0  MR  I TE  <  I DJT .  1 03  > 

00030430  DO  4  1=1 » HD 

00 03 0430  j  =  ND+l-I 

00 03 05 00  WRITE  .'I  OUT .  1  03)  A.J>  ,XI  >'J>  , Z  C.u  ,E1  <  J>  ,E<J<  .EFIN.  J) 

00030510  4  CONTINUE 

00030580  WPITEdOUT.  104> 

00  03  053 0  MR  I TE  < I OUT . 1 050 

•  10030540  DO  5  1=1 .ND 

00 03 0550  J  =  ND+l-I 

00030560  WRITE  <1  OUT.  103)  XI  <J> .  TOTSTR  .J) . EFFSTP < J <  .  UW  <  .1  •  ,  U i>  <  .i .  .  IJ  <  J) 

O0i-i305r0  5  CONTINUE 

00  030530  WRITE 'TOUT. 10?)  TIME.UCDN 

00030530  WRITE  .''l OUT.  11  0>  SETT. SPIN 

00 03 OF 00  WRITE ‘I OUT .111)  PUD251 

0 0 03 06 10  WRITE1 I OUT. 11 3>  ML 

00 03065 DC 

0009063GC  . • . FORMATS 

00  03  064  0  100  FORMAT  <  1  HI  •'^.•'•^■'1 4  ■  1H4  .•  .34HINITIAL  CONDITION*  IN  COMPRESSIBLE* 

00 03 0650  '<•  11H  FOUNDATION.  13 a HV.O 

00030660  101  FORMAT. '/3X. 5 <1H*) . 13H  COORDINATES  .  5  <1H*'  *  1 :3X,  5  C1H*  >  » 

00 03 06 7  0  $,  13H  VOID  RATIOS  *5<1H*>> 

0003 0680  105  FORMAT  <y?X» 1HA. 10X.2HXI. 1 IX. 1HZ. 7X. 8HEIN ITI AL. AX. 1HE.8X. 

00030630  i-,  6HEFIHAL) 

O008070P  1 03  FORMAT • ?X« 5 'FI  0. 4. £X) ,F1 0. 4) 

0003071  0  104  FORMAT  < ✓✓15X* 5 <1HV>  .  1  OH  STRESSES  .  5  ■ 1H»>  >  7X.  5  HH* )  ? 

00030750  S  16H  PORE  PRESSURES  »5'1H*V> 

00030730  1  05  FORMAT  •'-  6X.  5HXI .  3X»  5H TOTAL  ,  5X>  INEFFECTIVE*  5X.  5HTDTAL  .  6X* 

00030740  !{•  6HSTAT IC.  6X*  6HEXCESS) 

00030750  106  ROPMAT  •:  1 H 1  /s/s  1 9  <1 H*  >  » 34H IN I T I  AL  CONDITIONS  IN  DREDGED  FILL. 

00030760  «  1 3  <1 H*) > 

00030770  107  FORMAT-: -vi  OX*  7HT I  ME  =  .  El  0.4. 5X»  56HPE6PEE  OF  CONSOLIDATION  = 
00030780  S'  FI  0. 6' 

00030730  108  FORMAT  •:  1  HI •  •" 1 4  •:  1 H*>  .  34HCURPENT  CONDITIONS  IN  COMPRESSIBLE. 

00090800  11H  FOUNDATION.  13.  IH*V> 

00090810  109  FORMAT  y\W\/ ■'//,'  19  TIH*)  ,  34HCURPEN  T  CONDITIONS  IN  DREDGED  FILL. 
00090330  %■  1 9 •:  1  h*:>  ) 

00090330  110  FORMAT *V1 OX. 13HSETTLEMENT  =  , FI 0. 4. 5X> 19WFINAL  SETTLEMENT  =  , 

00030840  ?v  FI  0.4) 

0OO90850  Hi  FORMAT.' '1  OX. 57HR0TT0M  BOUNDARY  GRADIENT  =  .F10.4) 

00  09  0860  115  FORMAT.' '1  OX.  37HWATEP  LEVEL  ABOVE  BOTTOM  =  »F10.4> 

00090870C 

00090880C 

non* 0890  RETURN 

on 03 0900  END 
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Si.lPRni.iT  I  ME  DATA IN 


0  it  1  0  0  i'i  0  0 
full  00 ill  i.iC 
001 00020C 

O010O030C  ♦  DATAIN  READS  THE  DATA  FROM  A  PREVIOUS  PROGRAM  RUM  FROM  * 

001 000400  ♦  FILE  SO  THAT  FUTURE  CONSOLIDATION  CAM  PE  CALCULATED  ♦ 

001 000* or  ♦  WITHOUT  REDOING  ALL  PREVIOUS.  ♦ 

ooi oo of.  or 

Ail  1  0  00701" 

0010 0 OS  0  COMMON  DA  .DUu *  PUT'S  1 0.  PUDZ 1 1  *  BUPZ2 1 >  D? .  DZ 1 .  DO ,  E  0 , E  0 0  «  ELL ,  ELL  1 

0  0 1  0  I'i  0  9 "  r  GC  ,  GC 1 » GS»  G  S 1 ,  GSPL  ,  GSPF  *  GW,  HBL  «  HDF » HDF 1 » IN «  I  MS  >  I  OUT  , 

0  0 1  0  0 1  0  0  a.  I QUTS ,  L  EL » LDF  ,  MT I  ME .  NPD I V , MPD I V 1 ,  NPL , ND , HD  I V  *  ND I V 1 , 

OOI  001 1  0  •;  NFLAG?  NM«NPROB ,  HPT ,  NND"  NNN,  NTIME«  PKO«  l?0«  Cl  ,  SETT,  SETT  1 

001  0  0  ISO  SFIN- SF INI ,TAU»  TIME,  TPR I  NT.  UCON.  IJC0N1  .  VP  1 1 ,  Wl  .  ML  1 .  ZKO 

001001  30  A<1  01>  ,Rt  <1  n  .  AFC  101)  .  AF1  *1 1  )  ,Ai  PHA<51)  .  ALPHA  1  <51 )  . 

ooi  0014  0  *.  BETA 'SI)  >  PETA1  .'51>  ,BF<101)  • BE I'll),  DSDE 1 5 I ) »  DSDE 1 '51) 

001  00150  *•  E  <  1  0 1 )  »  E 1  <  1  0 1 )  ?  E 1  1  r  1 1  >  , EF I N  < 1 0 1>  .  EFIN1  <  1 1)  .  ER  '"1 1  )  « 

001  0  0  ISO  ft,  ES  <51) .ESI 05D . EFFSTR 0101) .EFSTP1 oil), F '101). FI  Oil.. 

001  00170  S  FI  NT  01  01  )  .FINT1  Cl  1)  ,PK  <51 )  .PK1  -'51  > ,RK <51 )  ,RK1  <51)  . 

iiO  1  00 1  SO  %,  PS  <5 1 )  ,PS1  <51 )  .  TOTSTP  <101  1  , T0STP1  Oil)  ,U<101)  ,U1  <11)  , 

00100190  ?.  UO'  lOl)  ,U01  <11  :• . UW  0101)  ,|JW1  <1 1)  .  XI  <1 01 > , XI 1 <  1 1 ) , 

0 i'll  00200  i.  Z  <  1  01)  ,21  <  1 1  > 

OOI 009 1 oc 

00100220  READ' INS, 100)  NST, IN. INS. I PUT. IQUTS, LBL. LDF 

i'i-OI  00.-30  READ  <  IMS  •  100)  NST,  NEDI V,  NEDI VI .  NDI V,  NBI V 1 ,  NEL 

001  OO.940  PEAD<INS,  1  00)  NST ,  ND*  NFL  AG,  NM.  NND,  NNN,  NT  IME 

00 1 0 Oc 5 0  READ  0 INS, SO 0)  NST , DA , DUDZ 1 1 , BUBZ2 1 , DZ , DZ 1 

00 1  0  OPE  (.1  READ  <  I  NS  >  £00)  NST,  EOO,  ELL ,  ELL  1 ,  GC ,  GC  1 

001  00270  READ  'IMS  « £Oii)  NST,  GSi,  GS 1 ,  GSBL ,  GS DF*  GW 

001 00380  READ < INS, 200  NST, HBL, HDF. HDF 1 . SETT, SETT  1 

001  00 S90  READ <  INS  ,200>  NST,  SFIN.  SFIN1 ,  TAU,  TF'PINT 

001  00300  REAP  <  INS.  200)  NST. UCDN, UCC1N.1 ,  VPI 1 ,  ML ,  WL  1 

001 00  31  0  READ  < INS, 200)  NST. DUO. DUDZ 10, DC, EO 

001 0  0  3  c  0  READ  <  I  NS'  ?  2  0 1.1)  NST ,  ZK  0 .  F‘k  0 , '"'  0 , 19 1 

001 00  33  or. 

001 00340  IF  'NEL  .EQ.  3)  GOTO  £ 

001 00350  DO  1  1=1, ND 

001  00360  READ < INS, £00)  NST, A  CD , AF  < I ) ,  BF  <I)  ,  E  CD  ,  El <I> 

00 1  ■>:>  OS  7  0  PE  AD  <  I  NS  ,  £  0  0)  NST,  EF  IN  ■  I )  •  EFFSTR  <  I )  ,  E  <  I )  .  F I  NT  Cl),  TOTSTP  <  I  > 

00)  00380  READ  <IMS  « £00)  NST.U'I)  •  IJ 0  <  I )  ,UW*  I)  «XI  Cl)  ,Z<I> 

001 00390  1  CONTINUE 

00100400  IF  'NFL  .EC.  2)  GOTO  4 


nni  0!i4i  or 

full  0O430  2  PD  3  1  =  1*  NDIV 1 

mil  00430  REAP (INS* 200)  NST.Al  <I)  .  AF  1  *  I)  *  BF 1  ( I)  ,  ER  ( I)  » El  1  (  I) 

00100440  REAP ('INS*  2C0>  NST, EF INI (I)  ,  EFSTP1 ( I) *  FI <I>.FINT1 (I) *  TQSTP1 <I> 

OO1O045O  REAP ( INS. 200>  NST.'Jl  <I> ,U01  < I >  . U(*>  1  < I>  . XI 1  ("I)  *  Z1  < I) 

00100460  3  CONTINUE 

00100470  IF  <NBL  .  EG.  3)  'SOTO  6 

00100420 C 

00100440  4  DO  5  I  =  1 , L  PF 

0 0 1 0 05 00  READ (INS* 20 0>  NST* ALPHA f 1 1 . BETA ( I ) . DSDE  < I > . ES  <1 >  * PK ( I > 

001O051O  REAP* I NS *2 DO)  NST, PK 'I > • PS <1 > 

00100520  5  CONTINUE 

001 00530  IF  <NBL  .EO.  2)  AOTO  S 

ooi 0054 or 

00100550  6  PD  7  I  »1  •  LBL' 

00100560  RE AD  a  NS  *  2  0  0)  NST  *  ALPHA 1 < I >  *  BETA 1 < I ) *  DSDE 1 < I )  *  ES 1 < I ) • FK 1  (I > 

O01OO57O  READ  (INS  *  200)  NST, RK1 <I> , RSI <I) 

00100580  7  CONTINUE 

001 00530C 

001 0060 DC  . . . RESET  TIME  CONTROL 

00100610  8  NM  =  NT  I  ME  -*•  1 

00100620  NTIME  =  NT IMS  +  MTIME 

00100630  WRITE flOUT, 300)  NPROB 

ooi  0 f»64 Of 

OO 1 006 5 of  . . . FORMATS 

0  ft 1 o  066  0  1  ft 0  FORMAT  *15,718) 

00100670  200  FORMAT *15, 5E13. 6) 

Oft! 00620  300  FORMAT (V9X, 30HC0NT INUATION  OF  PROBLEM  NUMBER* I4> 

OOl  00690C- 

00100700  RETURN 

00100710  END 

OOI 00720C 
OOI 0073 OC 


oril  i  norm 
nm  i  not  oc 
0,'U  1  002  OC 
ooi  i  o  03  or 

not  l  o ri4 os 
nm  10050C 
oon  on-,  nc 
ft 01 1 007 OC 
nm  1 0030 
nm  i  oo9o 
nni i  moo 
non  oi  10 
nm  l  oi?n 

001  1 01  sO 

nm  1 0140 
OOl 1 0150 
003  1  i'll  30 

nm  l  oi?o 
OOl 1 0130 
not i 0190 
001 1 0300 
OOl  i 031  DC 

no i l o££o 

001 1 0230 

001  1 0340 

nm  1 0950 
OOH  0230 
001 1 0270 
001 1 0280 
001 1 0290 
001 1 0300 
00110310 
001 1 0320 
001 1 0320 
001 1 0340 
00110350 
OOl 1 0380 

non  0370 

001 1 0320 
001 1 0390 
001 1 0400 
001 1 04 1 0 
0011 0420 
OOl  1 0490 
001 1  044 OC 


SUBROUTINE  SPVD3T 


♦  23V DPT  S3''E7  THF  PPT9  FROM  8  PREVIOUS  PPOGRPM  PUN  ON  ♦ 

♦  FILE  70  THPT  FUTURE  EXTENSION"  TO  THE  RUN  M3Y  BE  MPPE  * 

♦  WITHOUT  »EC3LCUL3TTNG  PREVIOUS  CON "Dl.  I  OPTION.  ♦ 


DP .  DU O.DUPZIO,  DUD 711.  DUDZ2 1 .  BZ«  D21  ■  DC'.  £0*  E00«  ELL . ELL  1 
GC , GC1 , GS, GS 1 . GSBL . G3PF , GW . HBL . HDF . HDF 1 . IN, INS . I OUT . 

1  OUTS ,  LBL.LDF  •  MTI ME, NBPI V.  NBDIV1 ,  NBL » ND,  NMV.  NDI VI , 
NFLPG,  NM.  NPROB »  NRTpNND.NNN.  NT IME.Pt1  0,00.  i*U  ,  SFTT.  SErTl 
SCIN. "F I M 1 . TPU .TIME. TPR I NT . UC  CN . UCON 1 , VR 1 1 • WL , Wl 1 . 210 
3' 1  01>  ,P1  <  n  .  ,3F  <1  OI',  .  HF  1  <  1 1  >  «  3LPH3  .51  >  .PLPH31  >51'. 
BET3 .5 1  >  ,  BET3 1  <5 1 , .  BF  *  1 0 1 ,  BF 1  <  1 1  > ,  DSPE  * 5 1 » ,  T,  SDF X  <5 1 » 
E<101> ,E1 <10r>,£ll <11>,EFIN- 101> ,EFIN1 *11 . EP  <1 1  •  , 

ES  *51''  .ESI  <5 1  >  •  EFFSTR  <  1  01>  .EFSTPl  *1  1>  ,  F  '  Of'  ,F1  <11  ■  , 
PINT <1 01 > .PINT  1 • 1 1> ,P>  <51> . Pf 1 <51, , RK  <5 1 , , PX 1 <51 , . 

RS  <51 >  »RS1 <51> . TOTSTP  *  1 01> , TOSTPt •  1 1 >  ,  U <1 01>  « U1 <11  • , 
HO*  1  01>  .U01  <11  >  . UI.J  <  1 0 1  >  »UW1  <1 1>  .XI  <101 ,  ,xi  1  *1 1  >  , 

2  <  1  0 1  > ,  2 1  <  1 1  > 


N'T  =  1 

WRITE <IOUTS» 100> 
NST  =  N'T  +  1 
WRITE  <  I  OUTS,  100', 
NST  =  NST  +  1 
WRITE < I OUTS. 1 00> 
NST  =  NST  +  1 
WRITE <IOUTS. 200- 
NST  ■  NST  +  l 
WRITE <IOUTS, 200, 
NST  =  NST  +  1 
WRITE <IOUTS,200> 
NST  =  NST  +  1 
MR  I TE  <  I  OUTS ,  8  0  CO 
NST  =  NST  +  1 
WRT  TE  < TOUTS. £0(0 
NST  =  NST  +  1 
WR I TE  < I OUTS « £  0  0> 
NST  =  NST  +  1 
WRITE < I OUTS. 200) 
NST  =  NST  +  1 
WRITE  '  I  OUTS  >200 


NST, IN, INS, I OUT, I OUTS, LBL.LDF 


NS  T , NB D I V, NED IV 1 , ND I V , N P I V 1 . NBL 


NS  T , N  P «  NFL3G • NM • NNp» NNN, NT I ME 


NST, P3, DUP21 1 . PUP221 . 02, D21 


NST,  E  0  (• .  ELL ,  ELL  1 ,  GC .  GC  1 
NST ,  GS>,  GS  1  ?  GSBL  ?  GS IiF,  GW 
NST, HBL , HDF , HDF 1 , SETT, SETT l 
NST «  SF IN, SF INI , T3U . T IME , TPR I NT 
NS  T . UC ON • UCON 1 » VP 1 1 • WL , WL 1 
NS T , DU  0 , DUDZ 1 0  >  DO , b  0 


NS  T,  2k 0, PK 0  »  0  0 . 0 1 


001  1045'.  IP  (NFL  .EG.  3>  GOTO  2 

0GH0460  do  1  i=i.nd 

Of.  110470  NST  a  NST  1 

001  1  04 SO  WRITE  <  IDUTS.200'*  NS  T » ft  *  I  >  »  ftF  ‘  I  > .  FF  <I>  tE  <  I>  »E1  <  I> 

001104*0  NST  =  NST  +  1 

001  1  0500  WRITE  fl  CUTS.  200>  NST.  EE  IN ..  I  >  ,  EFFSTP  .'I >  ,F  ■  I  "•  .PINT  a  >  »TDTS  TP  f  I  > 

oo  11  05 10  N''T  =  NST  +  1 

001  10520  WRITE- IDUTS.  £00)  NST.Ud  ■  .  U  0  •  I  ■'  .UW<I>  .XI  •  I  >  .  "  (I  • 

001 105 SO  1  CONTINUE 

00110540  IF  .'NFL  .EO.  2)  GOTO  4 

0011 055 OC 

00110560  2  DD  3  1=1 • NPIV1 

00 it  0570  NST  =  NST  +  1 

0011  0530  WRITE*  IDUTS.200'.  NST.  ft  1  <l>  .  ftFl  • I> .  BF1  d>  .  ER  <  I  >  .  El  1  (  P 

00110590  NST  =  NST  +  1 

001  1  OR  0  0  WRITE (IDUTS. 200>  NST.  EF  INI  (I> .EFSTFl  (I  1  .FI  •.  I  >  .FI  NT  1  ■  I  ■  .  TS:  T=  1  •  I  ■ 

001  !  ORl  0  NST  .=  NST  +  1 

001  1  06£  0  WRITE  (IDUTS.  200*  NST.  U1  <  I  • .  U01  .'I).  iJWl  <  I  >  , :  1 1  ■  I  >  .  £'  i  1 1  ■ 

001106-0  5  CONTINUE 

001  1064  0  IF  .'NFL  .EO.  3>  GDTO  6 

1*11*1 1 1  Gro  OC 

ii.il  1  0**0  4  DP  5  I-1«LDF 

nmioero  NST  =  H?r  ♦  i 

0011  0630  write  (TOUTS  .200)  NST.ftt_PH*a>  •  BETftd>  »  DSI'Ed  >  .ES  ' I>  .PR  (I  ' 

00110690  NST  =  NST  +  1 

nonoron  WRITE  dOUTS.  200>  NST. PL  ( I .)  . P S  ( I O 

00110710  5  CONTINUE 

00110720  IF  (NFL  .EO.  2>  RETURN 

001  1  07:3  Of 

00 110740  6  PD  7  I = 1 . LEL 

00110750  NST  =  "ST  +  1 

00110760  WPITF  (IDUTS.  20CO  NST.ftLPHftl  (I>  •  FETftl  (P.PSPE1  (I.)  .ESI  (I)  .PK'l  <l> 

00110770  NST  =  NST  +  1 

001 1 0730  WRITE • IOUTS»200>  NST.PR1 d>  .PS1 • I.) 

00110790  7  CONTINUE 

ooii os one 

not  1 031 OC  . . .  FORNftTS 

0 0 1 1  082  0  1  0  0  FCFMPT  ( 1 5 .  7 1 9~> 

00110330  200  FDPMftT (I5.5E13.6> 

0  011034  or 

00! 1 0350  RETURN 

00110360  END 


APPENDIX  C:  SAMPLE  PROBLEM  LISTINGS 

1.  The  following  pages  contain  sample  data  input  and  calculation 
results  from  the  two  practical  applications  previously  discussed. 


Cl 


2.  This  page  and  the  next  contain  the  input  data  file  used 
the  dredged  fill  with  compressible  foundation  example. 


1  no 

1  M  1 

£ 

1  Ml 

2  1 

2  n  m 

2.33 

20.0  2 

5 . 0  36  0 

20 1 

3.  MM 

0.  0 

1.21 0E-03 

;?0c 

2.35 

4.2 

1.1 12E-03 

?03 

2.30 

ft  ft 

1. 03GE-G3 

204 

2 . 35 

14.0 

9. 494E-04 

2.80 

13.3 

3. 854E-04 

?  0£ 

tct 
i>-  ■  i  J 

25.4 

3.234E-04 

?  ft  ? 

2.  70 

?«?■  ft 

7.616E-04 

.-'MS 

2.65 

33 .  0 

7. 000E-M4 

SMS 

3 . 6  0 

48.  0 

6. 382E-04 

S  t  M 

£.  55 

58 .  0 

5. 788E-04 

SI  1 

2. 5ft 

70.  0 

5.227E-04 

SIS 

2.45 

86 .  0 

4 . 68 0E- 04 

SI  3 

2.40 

l  04 . 0 

4. 234E-04 

S 1 4 

2.35 

128.  0 

3.  830E-04 

St  5 

2.  30 

1  54 .  0 

456E ~  04 

2  lb 

2.25 

190.  n 

3.  036-E-04 

•51  ^ 
fc-  A,  « 

2 . 2  0 

238.  0 

2. 736E-04 

Sis 

2.  15 

233.  0 

2.448E-04 

SI  9 

2.10 

344.  0 

2. 160E-04 

SSM 

2 .  05 

420.  0 

1 . 344E-04 

SS 1 

c*.  Oft 

510.  0 

1. 714E-04 

•TiT’  C 

1.35 

640.  0 

1 .512E-M4 

c  3 

1.30 

780.  0 

1 . 325E-04 

SS4 

1 . 35 

350.  0 

1 . 17 OF- 04 

cr, 

1.80 

1 1 6  0 .  0 

1. 034E-04 

sss 

1.75 

1 4  0  0 .  0 

3.  0O0E-05 

ft  2  7 

1.70 

1700.  0 

7. 720E-05 

2  ft 

1.65 

2 04  0 . 0 

6. 624E-05 

229 

1.60 

2540. 0 

5 . 832E- 05 

S3  0 

1.55 

31 00. 0 

5.  1 12E-05 

231 

1.50 

3750. 0 

4. 332E-05 

?3  2 

1.45 

4600.  0 

3.  773E-05 

S3  3 

1 .40 

5540.  0 

3.  137E-05 

S3  4 

1 . 35 

6800.  0 

2.  736E-05 

235 

1.30 

84  00.  0 

2. 333E-05 

236 

1 . 25 

1  0400.  0 

1.337E-05 

.300 

3.75 

3.  0  5. 

0  31  7.0 

301 

7.  00 

0.  0 

8.568E-03 

30£ 

6.95 

0.3 

8. 30SF-03 

303 

6.90 

1.0 

7.S4SE-03 

304 

6.80 

3.3 

7. 300E-03 

305 

6.  60 

5.4 

6.  091F-03 

3  0(3 

6.40 

8.  8 

5.  098E-03 

307 

6.30 

13.3 

4.  176E-03 

.3  0  3 

6.  00 

18.3 

3.  442E-03 

309 

5.80 

34.8 

3.833E-03 

3t  U 

5.60 

33.3 

2.  318E-03 

311 

5.40 

44.  0 

1 . 386E-03 

31  £ 

5.30 

57.  0 

1 . 570E-03 

313 

5.  0(- 

73.  0 

1 .267E-03 

314 

4.80 

96.  0 

1. 037E-03 

315 

4.60 

135.  0 

8.  358E-04 

316 

4.40 

163.  0 

6.  768E-04 

.31 7 

4.30 

31  0.  0 

5.429E-04 

313 

4.  00 

374.0 

4. 378E-04 

319 

3.80 

358 .  0 

3.49-9E-04 

33  0 

3 . 6  0 

468 .  0 

8.  794E-04 

331 

3.40 

600.  0 

3.  218E-04 

■'u'  iZ 

3.30 

790.  0 

1 . 735E-04 

333 

3.  00 

1 030. 0 

1.354E-04 

334 

3.80 

1330. 0 

1. 082E-04 

335 

3.60 

1740. 0 

7.  488E-05 

336 

3.40 

3340. 0 

5.  323E-05 

337 

3.30 

3000. 0 

3.  696E-05 

333 

3.  00 

4000. 0 

2. 506E-05 

•339 

1.80 

5480. 0 

1.656E-05 

93  0 

1.60 

7500. 0 

1. 094E-05 

331 

1.50 

9000. 0 

3.  7S4E-06 

40  0 

0.  65 

3. OF- 04 

6.  0 

40 1 

6  10 

1.0  4 

403 

365 

3.  0  36 

.  0 

403. 

730 

3.  0  38 

.  0 

404 

1  095 

3. 0  39 

.  0 

405 

1460 

1.0  30 

.  0 

C3 


62 


F 

t 

r. 


3.  Below  are  the  calculation  results  after  2  years.  A  total  of 
6.0  ft  of  dredged  material  has  been  deposited.  Results  for  the  compres¬ 
sible  foundation  are  not  shown. 


►CURRENT  COMP  IT  ION?  IN  DREDGED  FILL* 


1 


♦♦♦♦♦  chord inates 


A 

X! 

6. 0000 

5.851  0 

5.5000 

4 . 7628 

5. 0000 

4.8904 

4.5000 

3. 8889 

4. 0000 

3. 3764 

3.5000 

2.9319 

3. 0000 

2.4945 

8. 5000 

2. 063? 

8. 0000 

1.6393 

1 .5000 

1.8209 

1 .  ft  ft  ft  ft 

0.  8093 

0 . 5  ft  ft  ft 

ft  .401  3 

0 . 

ft. 

7 

einitial 

0. 7500 

7. 0000 

0.6875 

7. 0000 

0.6250 

7. 0000 

0.5685 

7. 0000 

0.5000 

7. 0000 

0.4375 

7. 0000 

0.3750 

7.0000 

0. 3185 

7. 0000 

0.8500 

?.  0000 

ft. 1875 

?. 0000 

0. 1250 

7. 0000 

ft. 0625 

7.  00 00 

0. 

7. ft  ft  ft  0 

VOIP  PATIOS  ♦♦♦♦♦ 


E 

EFINAL 

7.  0000 

7. 0000 

6 . 66 1  0 

6.5168 

6.4651 

6. 1820 

6. 3076 

5.931 1 

6. 1735 

5.7405 

6. 0540 

5.5889 

5.9439 

5.4565 

5.8405 

5.3419 

5. 7488 

5.2369 

5.6478 

5. 1447 

5.5561 

5. 0594 

5.4661 

4 . 982  0 

5. 377? 

4. 9826 

♦♦♦♦♦  stresses  ♦♦♦♦♦ 


♦  ♦♦♦♦  F'OFE  PRESSURES  ♦♦♦♦♦ 


XI 


C- 

251  0 

4 

T  oc‘3 

4 

2904 

3889 

3 

3764 

d 

9319 

4945 

c 

06  3 ’3 

1 

6?9:3 

t 

8809 

8  03  3 

f 

4013 

0. 


total 

58.  9769 
98.  868 9 
138.5668 
168. 1  90S 
803.8485 
.???.  9 1?  5 
871.9313 
305.8344 
938.9444 
37 1 . 8787 
404.4505 

468. 5374 


EFFECTIVE 
-0. 00 CO 
4.4541 
7.6939 
1 0. S3: O 
13. 86 1 3 
16.8505 
80. 052* 
33. 4638 
87.8259 
31.1940 
35.5784 
40.4323 
45. 4778 


TOTAL 

58. 9769 
91 . 8089 
184.8749 
157.3578 
189.  397c' 
220.9630 
251 . 8790 
238. 1712 
311.7185 
340.6847 
368.8781 
396 .  '8359 
483. 0596 


STATIC 
58. 9769 
39 . 4379 
118.9138 
147.7158 
175.9495 
808.6885 
830.9818 
857.3594 
884.3444 
310.4537 

336.8005 

361.593? 
386 ■ 6: 7 4 


EXCESS 

n. 

8. 3709 

5 . 956 1 

9.6480 

1 3. 438? 
17.8745 
80.897c 
24.3118 
87 . 3741 
30.2310 
32 . 6776 
34 . 6487 
36 . 4cc'8 


TIME  =  ft. 73 ft OE  02 


DEGREE  DF  CONSOL  I BAT I ON  = 


settlement  » 


0.7490  FINAL  SETTLEMENT  = 


0. 708691 
1 . 0569 


BUTTON  BOUNDARY  GRADIENT  = 
MATER  LEVEL  ABOVE  BOTTOM  = 


C4 


i 


t 


4.  This  page  and  the  next  contain  the  results  after  8  years  of 
consolidation.  A  total  of  14.0  ft  of  dredged  fill  has  been  deposited 


♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•: URRENT  COND I T  ION  ~ 


IN  DREDGED 


♦♦♦♦♦  COORDINATES  ♦♦♦♦♦  ♦♦♦♦♦  VOID  RAT TO?  ♦♦♦♦♦ 

H  H'l  Z  E  INITIAL  E  PRIMAL 


14. no o o 

11.3111 

1.7500 

13. 5  noil 

10.3550 

1 > 6875 

13. noon 

10.3570 

1 . 6550 

18.5000 

9.3011 

1 . 5685 

12. nnoc 

9.4545 

1.5000 

11.5000 

9. 0155 

1.4375 

1 1 . 000 0 

P.5830 

1 > 3750 

10.5000 

8.  1559 

1.3185 

1  0.  000 ii 

7 .  73  ?*: 

1.8500 

9 .  5  0  0  0 

7.3158 

1.1875 

•?.  0000 

6.9050 

1. 1850 

3.5000 

6.4950 

1 . 0685 

3. 0000 

6. 0357 

1. 0000 

7*  •  5  0  f*  0 

5 . 6859 

0. 3375 

7. 0000 

5.8834 

0.8750 

6.50O0 

4.8873 

0 . 8 ! 85 

6 .  0  0  o  0 

4.4943 

0.7500 

5.50  0 0 

4. 1045 

0.6875 

5.  000  0 

3.7178 

0.6850 

4.50  0  0 

3.3340 

0.5685 

4  .  0  O'  0  0 

£.9531 

0.50 0  o 

3.5000 

£.5750 

0.4375 

3. 0000 

8. 1996 

0.3750 

£.5000 

1 . 8868 

0.3185 

E.  00 of. 

1 . 4566 

0.8500 

1 . 5000 

1 . 0889 

0. 1875 

1.00  0 o 

0. 7836 

0. 1850 

0 . 5  0  0  0 

0  ■  3  c*  06 

0.  06»?5 

0. 

0. 

0. 

7. 0000 

7. 0000 

7. 0000 

7. 0000 

6.6041 

6.5162 

7. 0000 

6 . 3886 

6. 1880 

7. 0000 

6.8148 

5.9311 

7. 0000 

6. 0808 

5.7405 

7. 00 00 

5.9696 

5.5889 

7. 0000 

5.8756 

5.4565 

7. 0000 

5 . 7934 

5 . 34 1 9 

7. 0000 

5.7196 

5 . 836  9 

7. 0000 

5.6580 

5. 1447 

7. 0000 

5.5898 

5. 0594 

7. 0000 

5.5898 

4.9880 

7. 0000 

5.4789 

4.9886 

7. 0000 

5.4179 

4.8633 

7 .  0  0  0  0 

5.3646 

4.80  39 

7. 0000 

5.3187 

4.7560 

7 .  00  0  0 

5.8680 

4.7090 

7. 0000 

5.8183 

4.661 9 

7. 0000 

5. 1640 

4.6148 

7. 0000 

5.1171 

4.5754 

7. 0000 

5. 0718 

4.5395 

7. 0000 

5.  0278 

4.5036 

7. 0000 

4.9851 

4 . 4676 

7. 0000 

4.9436 

4.4317 

7. 0000 

4.9033 

4 . 3966 

7. 0000 

4.3641 

4. 3676 

7. 0000 

4.8259 

4.3335 

7. 0000 

4.7886 

4. 3095 

7. 0000 

4.7519 

4.2304 

C5 


Var  i  lit  i  - 


ro  ro  to  *j«j  co 


♦♦♦♦♦  STRESSES  ♦♦♦♦♦ 


XI 

TOTAL 

EFFECTIVE 

t 1.3111 

75. 3995 

-0. O0O0 

10.3850 

112.5595 

5.3368 

10.3570 

143.5838 

9. 1334 

3.801 1 

183. 663? 

18.8338 

8.4545 

213.5535 

16. 1951 

3.  0155 

£58. 7694 

19.8080 

3.5330 

£36.5331 

88.3060 

3. 1553 

320. 0644 

85. 0774 

7 . 7336 

353. 2375 

88. 1767 

7.3158 

336. 1354 

31.0139 

6.  8080 

418. 7796 

33.7318 

6.4880 

451.1359 

36.9892 

6. 0357 

433. 3655 

40. 0640 

5. 6689 

515.3870 

43.  0.321 

5.8834 

547. 0774 

46.3021 

4 . 8373 

578. 6886 

49.6774 

4.4943 

609.9673 

58.9724 

4. 1045 

641. 1173 

56. 1973 

3.  7173 

678. 0756 

59.8386 

3.3340 

708.8433 

63.6297 

£.8531 

73.3.  4418 

67.8576 

S.  5750 

763.8599 

70.7771 

8.  1896 

794. 1 096 

74.7173 

1 . 3868 

384. 1951 

79.4359 

1 . 4566 

354. 1818 

34. 1177 

1 . 0339 

839. 3984 

88 . 6852 

0. 7836 

913.5187 

93. 0187 

0 .  36  Or- 

948 . 9*57 

97.6540 

tf . 

978.3146 

102.9677 

TIME 

*  0. 298 OE 

04  DEGREE 

SETTLEMENT  =  £.6989 

BOTTOM  BOUNDARY  GRADIENT  = 
WATER  LEVEL  ABOVE  BOTTOM  = 


♦♦♦♦♦  PORE  PRESSURES  ♦♦♦*♦ 


TOTAL 

STATIC 

EXCESS 

75.3995 

75.3995 

0. 

107.2238 

105.7345 

1 . 4338 

139.4047 

134.9332 

4 . 4666 

170.9750 

163.3382 

7.5363 

802.3584 

191.2535 

11.1 049 

833. 5673 

218.6444 

14.9230 

£64.8881 

245.6381 

13.6440 

894.9371 

272.8394 

88.6977 

385. 0607 

293.6375 

26.4233 

355. 1215 

324.7104 

30.4111 

384.9984 

350.5296 

34. 4683 

414. 1968 

376.1109 

33. 0353 

443.3015 

401.4655 

41.3360 

478.8950 

486 . 6  02  0 

45.6930 

500.7754 

451.5874 

49.2479 

528.9451 

476 . 2476 

52.6976 

556.9954 

500.7673 

56.8876 

534.9196 

525. 0983 

59.3272 

612. 1931 

549.8856 

6.2 . 96  74 

639.2135 

573. 1733 

66. 0453 

666. 1 836 

596.9418 

69. 848 

69.3.  0823 

680.5349 

78.5479 

719.3922 

643.9596 

75.4387 

744. 7092 

667.8801 

77.4891 

770. 0035 

690.3218 

79. 6823 

795.8672 

713.8674 

31.9998 

820.4940 

736. 0687 

34.4313 

345.3317 

753.7103 

86.6209 

369. 3469 

731.8146 

S3. 1323 

OF  CONSOLIDATION  *  0.768640 

FINAL  SETTLEMENT  =  3.4838 

-S3. 3337 
IS. 5135 


( 
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5.  After  14  years,  conditions  In  the  dredged  fill  layer  are  as 


shown  below  and  on  the  next  page. 


♦♦♦♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦CUPPEhT  COND I T I QMS 


IN  DREDGED 


COORBINRIES  *♦♦♦♦ 


A 

XI 

2 

14.  II (Ml It 

1 0. 3073 

1.7500 

13.5  ('  ft  0 

10. 3E 52 

1 . 6975 

13.  i'i  (i  (i  0 

9. 86~4 

1.6250 

1?. 50 MO 

9.  4231 

1 . 5625 

1  £ .  00 (ill 

3.3916 

1.5000 

11.5000 

•5.  5694 

1 . 4375 

1  1 . 00  0  0 

3. 154° 

1. 3750 

10.5000 

7. 7493 

1.3125 

to. 0000 

7.3445 

1.2500 

3.5000 

6.9474 

1. 1S75 

'1.  00  0  0 

6 . 5549 

1. 1250 

8 .  5  0  0  0 

6. 1666 

1 . 0625 

S. 0000 

5.7323 

1.0000 

7.5000 

5.4016 

0.9375 

7. 0000 

5. 0242 

0.  8750 

6.5000 

4.6499 

0.2125 

6. 0000 

4.2784 

0. 7500 

5. 5000 

3.9  096 

0.6375 

5. 0000 

3.5433 

0.6250 

4.500  0 

3. 1795 

0.5625 

4 .  0  0  0  0 

£,•31 79 

0.5000 

3.5000 

2. 4586 

0.4375 

3. 0000 

2. 1014 

0.3750 

5. 5000 

1 . 7462 

0.3125 

3. 0000 

1 . 393 1 

0.2500 

1 .5000 

1 . 0420 

0. 1375 

i. oc on 

0.6923 

0. 1250 

0.50  0  0 

0. 3455 

0. 0625 

0. 

0. 

0. 

♦♦♦♦* 

VOID  PftTIO: 

EINITIftL 

E 

EFINftL 

7. 0000 

7. 0000 

7. 0000 

7. oono 

6.4902 

6.5162 

7. 0000 

6.2046 

6. 1820 

7. 000H 

5.9903 

5.9311 

7. 0000 

5.3239 

5.7405 

7. 0000 

5. 6904 

5.5829 

7. 0000 

5.5734 

5.4565 

7. 0000 

5.4312 

5.3419 

7. 0000 

5.3945 

5.2369 

7. 0000 

5. 3161 

5. 1447 

7. 0000 

5.2447 

5.  0594 

7. 0000 

5. 1791 

4.  988  0 

7. 0000 

5. 1 192 

4.  9226 

7. 0000 

5. 0642 

4 . 863  3 

7. 0000 

5. 0134 

4.30 

7. 0000 

4 . 966 1 

4. 756m 

7. 0000 

4.9213 

4.7090 

7. 0000 

4.3802 

4 . 66 1 9 

7.  ft  0(i0 

4.3407 

4.6143 

7. 0000 

4.3031 

4.575.4 

7. 0000 

4.7671 

4.5395 

7. 0000 

4. 7322 

4.5035 

7. 0000 

4.6985 

4.4676 

7. 0000 

4.6657 

4.4317 

7. 0000 

4.6337 

4. 3966 

7. 0O00 

4. 6026 

4.3676 

7. 0000 

4.5721 

4.3385 

7. 0000 

4.5423 

4.3095 

7. 0000 

4.5131 

4.2804 

C7 


STRESSES  ♦♦♦♦♦ 


XI 

TOTAL 

EFFECTIVE 

1  0.  307  3 

180. 0790 

-0. 0000 

1  0. 5252 

158. 3*18 

7.2659 

9.3884 

1*2.4441 

13. 099* 

9.4831 

888.  9297 

13.5267 

8.9918 

880.6797 

24. 0093 

8.58*4 

£93.8490 

29.404? 

8.  154? 

328.542? 

34.3643 

7.7488 

358.83  0 1 

39. 6160 

7. 3445 

390. 75*7 

44 . 35*6 

8 .  '*•*  7  4 

422 • 3678 

49.4526 

8.554* 

453. 624 3 

54. 096? 

8. I860 

484.7339 

58. 6691 

5.7383 

515.5388 

63.464* 

5.4018 

548.  1 2  00 

67.3614 

5. 0848 

578. 4358 

71. *835 

4.649* 

806.47*1 

76.9004 

4.2784 

31.9397 

3. 9096 

ir.ir-tr* . 

66. 730? 

3.5433 

698. 2088 

91.3134 

3. 1795 

725.732? 

*5.6393 

2.3179 

755.11*8 

100.773* 

2.4588 

784. 3675 

105.3860 

8.1014 

313.4981 

110.7213 

1.7482 

84£ . 4666 

115.4777 

1.3931 

871 . 3854 

120. 1077 

1 . 0480 

900. 061 0 

1 £4 . 6248 

0.8928 

988. 6766 

130.2925 

0. 3455 

957. 1746 

135.9580 

0. 

985.5574 

141.515) 

PORE  PRESSURES 


total 

STATIC 

EXCESS 

120. 0790 

120. 0790 

0. 

149.7253 

150. 1662 

-0. 440* 

179.3443 

178.7941 

0.5501 

203.4030 

206.4547 

1 . 9483 

236. 6699 

833. 3797 

3.8908 

264.4443 

259.7240 

4.720? 

898. 1784 

285.592? 

6. 585? 

319.2141 

311. 0551 

8. 1590 

346.4000 

336. 159? 

10.2404 

378.9152 

360.9428 

1 1 . 9724 

399.5875 

335.4343 

14. 1532 

426. 0647 

4  09 . 6589 

16.405* 

458. 0739 

433.6388 

18.4351 

478.2586 

457.3951 

20. 3636 

504. 566? 

480. *452 

23.6215 

529.7787 

504. 3041 

85.4748 

554.6949 

527.4846 

27.210? 

579.7481 

550.4978 

8*. £443 

604. 8844 

573. 3528 

31.531c 

630. 0*35 

596. 0577 

34. 0357 

654. 3453 

613. 61*2 

35. 786) 

673.5415 

641. 0425 

37.49*0 

708.7608 

663.3321 

39. 4832 

726.9398 

685.4918 

41.4*73 

751.2173 

707.5254 

43.6983 

775.4369 

729.4360 

46. 0008 

7*8.3341 

751.2265 

47. 1575 

821 .2166 

778.8996 

43.3170 

844.  0-*84 

794.4574 

49.5849 

TIME  =  0.511  OF  04  DEGREE  OR  CONSOLIDATION  *  0.91 £65* 

SETTLEMENT  =  3.192?  FINAL  SETTLEMENT  =  3.4982 

BOTTOM  BOUNDARY  GRADIENT  =  -£2.4485 

MATER  LEVEL  ABOVE  BOTTOM  =  12.731? 


6.  The  input  data  file  for  the  soft  compressible  layer  example 
is  given  below. 


J  00 

11  1 

C 

101 

1  3 

.=•0  0 

£.80 

20.0  2 

1.0  36  75 

=  01 

3.  00 

0.  6 

1.81 0E-03 

.=08 

£.95 

4.2 

1. 112E-03 

P  03 

8.90 

8.8 

1. 030E-03 

804 

£.85 

14.  0 

9. 494E-04 

205 

£.  80 

19.6 

3.35412-04 

flH 

£ .  75 

85.4 

8. 234E-04 

80? 

£.70 

38.  0 

7. 616E-04 

.=  08 

8.  65 

39.  0 

7. 0O0E-04 

.=*09 

£.60 

48 .  0 

6. 392E-04 

21  0 

£ .  55 

53.  0 

5. 738E-04 

=  11 

2.50 

70.0 

5. 287E-04 

212 

2.45 

86.  0 

4. 6SUE-04 

813 

2.  40 

1  04.  0 

4. 234E-04 

=  14 

£.35 

188.  0 

3. 330E-04 

315 

8.30 

154.  0 

3. 456E-04 

=  16 

2.85 

180.  0 

3. 096E-04 

.=  1? 

£.20 

£  38.  0 

£. 736E-04 

;=  1 8 

2.  15 

238 .  C 

8.448E-04 

819 

8.10 

344.  0 

2 . 160E-04 

=20 

8.  05 

420.  0 

1 . 944F-04 

?£  1 

c.  00 

510.0 

1.714E-04 

OOO 

>.  C. 

1 . 95 

640.  0 

1 . 512E-04 

823 

1.90 

780.  0 

1. 325E-04 

=24 

1 . 35 

950.  0 

1. 170E-04 

=£5 

1.80 

1160. 0 

1. 034E-04 

=89 

1.75 

1400. 0 

9. 0 ODE -05 

227 

1.70 

1700.  0 

7. 720E-05 

828 

1.65 

8040.  0 

6. 624E-05 

829 

1.60 

2540. 0 

5.332E-05 

830 

1.55 

31 00. 0 

5. 1 18E-05 

831 

1.50 

3750.  0 

4. 392E-05 

238 

1.45 

4600.  0 

3. 773E-05 

833 

1 . 40 

5540.  0 

3. 197F-05 

834 

1 . 35 

6800.  0 

£. 736E-05 

p.5=r, 

1.30 

3400.  0 

£. 333E-05 

836 

1.25 

1  0400.  0 

1 . 997E-05 

30  0 

0  0 

0  1  0 

62.4 

901 

0  0 

0 

400 

0.6 

i.oe-03 

1.0 

401 

0  10 

1 .  o  ? 

408 

365 

500.  0 

21.0 

403 

730 

500.  0 

£1.0 

404 

1  095 

0.  0 

81.0 

500.  0 
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7.  Conditions  in  the  compressible  layer  after  3  years  are  shown 
below.  The  total  layer  depth  differs  from  the  input  value  slightly  due 
to  the  iterative  method  of  calculating  the  material  coordinate  and  the 
fact  that  Lagrangian  coordinates  are  reset  to  match  the  material 
coordinate. 


•♦♦♦♦♦♦♦♦♦♦♦•♦CURRENT  CONDITIONS-  IN  COMPRESSIBLE  FOUNDAT I □!<♦♦♦♦♦♦♦♦♦♦♦♦♦ 


♦♦♦♦♦ 

COORDINATES 

♦♦♦♦♦ 

♦  ♦♦♦♦ 

VOID  PATIOS 

♦  ♦♦♦♦ 

A 

VI 

■7 

EINITIAL 

E 

EFINAL 

19. 9825 

13. 5992 

6.4379 

2. 4344 

1.7203 

1.7208 

17. 7837 

16. 7999 

5.3391 

2.3113 

1 . 8267 

1.7087 

1?. 6638 

14.9337 

5. 1903 

2.2134 

1.9213 

1 . 6969 

1 3. GOG? 

13. 01 64 

4.5415 

2. 1450 

1.9817 

1 . 6862 

1 1 .  5864 

1 1 . 074? 

3.S92S 

2. 0958 

1.9970 

1 . 6755 

•a.  s??  3 

9. 1359 

3.2440 

2. 0392 

1 . 9744 

1 . 6648 

7. £41? 

7.2803 

3. 5952 

1.9991 

1 . 9274 

1 . 654 1 

5. 7047 

5.-3413 

1.9464 

1.9711 

1 . 3627 

1 . 6455 

3. 736 1 

3.5083 

1 . 2976 

1 . 9436 

1 . 7358 

1 . 6382 

1.3348 

1 . 7879 

0. 6433 

1.9176 

1.7027 

1.6309 

0. 

0. 

0. 

1.3930 

1 . 6253 

1 . 6236 

♦♦♦♦♦  STRE 

SSES  ♦♦♦♦♦ 

♦♦♦♦♦ 

POPE  PRESSURES 

♦♦♦♦♦ 

XI 

TOTAL 

EFFECTIVE 

TOTAL 

STATIC 

EXCESS 

18. 

5'-?2 

1724. SI  0 £ 

1575. 0000 

149.3102 

149.8102 

-0. 0000 

16. 

7C199 

1909. 9566 

1047.8234 

862. 1332 

262. 0843 

600. 0439 

14. 

9337 

2099.2790 

720.3325 

1378.9465 

378 . 5344 

1 000. 41 21 

13. 

0 1 64 

2291.7961 

557.5311 

1734.2150 

493. 1792 

1236. 0358 

11. 

074? 

2485. 8404 

517.9196 

1967.9208 

619.3512 

1348.5696 

9, 

1359 

2679.6324 

576. 6408 

2103. 0416 

74  0. 32  09 

1362.7207 

7. 

2203 

2872. 0377 

703.3944 

2168.6934 

859.3540 

1308. 3394 

er 

341  3 

3062.2064 

906.7564 

2155.4501 

977. 1 003 

1 178.3493 

3. 

503:3 

3249.4534 

1228.2871 

2021. 1713 

1091.4300 

929. 6913 

1. 

72 79 

3433. 4-306 

1 683.  6487 

1749.7819 

1202.5799 

547.2020 

0. 

3614. 1230 

2292. 1446 

1331.9784 

1310.4000 

21.5784 

TIME 

=  0.1095E 

04  DEGREE 

OF  CONSOLIDATION  =  0. 

f 

516665 

SETTLEMENT  = 

1 . 3839 

FINAL  SETTLEMENT  =  2 

.6774 

BOTTOM  BOUNDARY  GRADIENT  =  90'?.  1478 

MATER  LEVEL  ABOVE  BOTTOM  =  81.0000 
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8.  Compressible  layer  conditions  after  6  years  are  shown  below. 


♦♦♦♦♦♦♦♦♦♦♦♦♦♦CURRFNT  CONDITIONS  IN  COMPRESSIBLE  FOUNDATION 


COORDINATES 

♦  ♦♦♦♦ 

♦  ♦♦♦♦ 

VOID  RATIOS 

A 

XI 

2 

EINITIAL 

E 

EFINAL 

18.9825 

17.9417 

6.4879 

2. 4844 

1.7203 

1 . 7208 

17.7887 

18.  16-32 

5. 839 t 

2.- 31 19 

1 . 7615 

1 . 7087 

15. 8892 

14.3596 

5. 1903 

2.2134 

1.7970 

1.6970 

1 9 . 8085 

12.5362 

4.5415 

2. 1450 

1.3217 

1 . 6862 

1 1 . 5284 

10.7016 

3.3928 

2. 0852 

1.8308 

1 . 6755 

9.5993 

8. 8668 

3.2440 

2. 0392 

1.3221 

1 . 6648 

7.8415 

7.  04.31 

2.5952 

1.9991 

1 . 7972 

1.6541 

5.7047 

5.2398 

1.9464 

1.9711 

1.7601 

1 . 6455 

3.7881 

3 • 463  3 

1 . 2976 

1 . 9436 

1.7154 

1 . 6  382 

1 . 8848 

1.7163 

0. 6488 

1.9176 

1 . 6686 

1.6309 

0. 

0. 

0. 

1.3930 

1 . 6248 

1 . 6236 

♦♦♦♦♦  STRESSES  ♦♦♦♦♦  *♦♦♦♦  PORE  PRESSURES 


XI 

TOTAL 

EFFECTIVE 

TOTAL 

STATIC 

EXCE 

17.9417 

1765.8365 

1575. 0201 

190.8165 

190.8365 

-0. 

16-.  16-32 

1949. 6899 

1 344.8958 

6  04.  794  0 

301.31 78 

302. 

14. 3596 

2135. 1 075 

1174.5330 

960.5745 

414.3622 

546. 

12.5362 

2321.7578 

1 068. 8232 

1252.9346 

528. 1404 

724. 

10.7016 

2509.1119 

1 030. 6543 

1478.4571 

642 . 6224 

g  js. 

8.8669 

2696.4732 

1067. 0027 

162®. 4705 

757.1116 

872. 

7 .  04  ?  1 

2883. 1429 

1173.-3985 

1709.7444 

870.9092 

838. 

•  L  7  ” 

3  '•■66.5421 

1351.6883 

1716.8537 

983. 436: 

7  33 . 

3. 4633 

3252 . 267 0 

1607.7970 

1644.4700 

1094.2380 

55  0. 

1.716® 

3434. 1201 

1913.7976 

1520. 3225 

1203.2692 

317. 

0. 

3614. 1232 

2292.0714 

1322. 0513 

1310.4  0  0 0 

1  1  . 

TIME 

=  0.2190E 

04  DEGREE 

□F  CONSOLIDATION  =  0. 

1 

SETTLEMENT  »  2.0408  FINAL  SETTLEMENT  =  2.8774 

BOTTOM  BOUNDARY  GRADIENT  =  491.8586 

WATER  LEVEL  ABOVE  BOTTOM  =  21.0000 


i_n  t.n  co  *-*  '.o  lh  v.i  '£}  n  ro  •  )  ro  o.<  w  _n  o>  “-j  co  •£■  o  ro 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Cargill,  Kenneth  W. 

Consolidation  of  soft  layers  by  finite  strain  analysis  / 
by  Kenneth  W.  Cargill  (Geotechnical  Laboratory,  U.S. 

Army  Engineer  Waterways  Experiment  Station).  --  Vicksburg, 
Miss.  :  The  Station  ;  Springfield,  Va.  :  available  from 
NTIS,  1982. 

113  p.  in  various  pagings  ;  ill.  ;  27  cm.  —  (Miscellaneous 
paper  ;  GL-82-3) 

Cover  title. 

"March  1982." 

Final  report. 

"Prepared  for  Office,  Chief  of  Engineers,  U.S.  Army 
under  CWIS  Work  Unit  No.  31173,  Task  34." 

Bibliography:  p.  63-64. 

1.  Computer  programs.  2.  CSLFS  (Computer  program). 

3.  Difference  equations.  Nonlinear.  4.  Soil  consolidation. 

I.  United  States.  Army.  Corps  of  Engineers.  Office  of  the 
Chief  of  Engineers.  II.  U.S.  Army  Engineer  Waterways 


Cargill,  Kenneth  W. 

Consolidation  of  soft  layers  by  finite  strain  :  ...  1982. 

(Card  2) 

Experiment  Station.  Geotechnical  Laboratory.  III.  Title 
IV.  Series:  Miscellaneous  paper  (U.S.  Army  Engineer 
Waterways  Experiment  Station)  ;  GL-82-3. 
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